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Synopsis 


Joints are inevitable in structural applications. Adhesively bonded joints are 
being preferred to mechanical fasteners because of distinct advantages they offer, 
such as uniform stress distribution, stiff connection, smooth surface contour, 
sealing against corrosion, excellent fatigue properties, no fretting problems, 
relatively light weight, damage tolerance and stealthiness. 

Use of composite materials in structural applications is also growing because 
of their excellent properties like high strength to weight ratio, stiflbiess to 
weight ratio, dimensional stability, electncal insulation, smooth outer surface etc. 
Composite materials provide a scope to the designer to tailor the material and/or 
structure to suit the applications. In addition, many complicated structural 
parts can easily be moulded using composites. Adhesive joints are found to 
be more suitable with composite materials particularly in case of the fibre 
reinforced plastics. With the development of various structural adhesives, these 
joints are being employed in a wide variety of structural applications ranging 
from aerospace, automobile components to those found in microelectronic 
components. In the primary adhesively bonded structures technologr (PABST) 
programme, an expert group is working to promote the use of bonded structures 
in aircraft industry. Some of the important applications of adhesively bonded 
structural members are aircraft’s rudder floor panel, engine cowling, elevator 
and flq), ceramic tyle of re-entry vehicles, integrated fuel tank, composite 
propeller blade, helicopters, launch vehicle attachjqaent (LVA) ring, collapsible 
tube mast. For structural integrity, it is, therefore, important to ass^ their 
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strengths and durability a priori so that joints could be designed in a rational 
manner resulting in an improved strength and more reliable structure. 

From a critical review of literature it is observed that polymeric adhesive 
materials generally exhibit material non-linearities (elasto-plastic and elasto- 
visco-plastic behaviour). In view of the versatility of finite element method 
this seems to be the best mathematical tool to analyse the bonded joints 
in presence of non-lineanties . It is also a convenient tool for associated 
parametric studies. Furthermore, fracture mechanics approach results in better 
mathematical model to predict the joint strength. Most of these investigations 
reported are either purely theoretical or computational and very few have 
experimental data to support their results. Further, from a comparison of 
results of various researchers it is observed that there is a large amount of 
scatter in the results^*^. Thus, there is a need for a better technique with a 
suitable mathematical model representing crack propagation for the analyses of 
adhesively bonded joints. It is also necessary to conduct simple experiments to 
determine the joint strength and compare them with those obtained from the 
finite element analysis. This will be complementary to the results obtained by 
the finite element technique. 

The thesis includes the studies of adhesive joint based on (i) linear elastic 
analysis (ii) elasto-plastic analysis and (iii) elasto-visco-plastic analysis. The 
adherend and adhesive materials are considered to be linear and non-linear, 
respectively. In all the above cases detailed parametric study is carried out. 
Experiments have been conducted on specially designed specimens to determine 
the fracture parameters and the strength of the joints. 

In Chapter 1, various types of prevalent adhesively bonded joints are 
discussed. This Chapter provides an overview of the joints, composite materials, 
non-linear finite element analysis, fracture mechanics approach and experiments 
to det^mine the joint strength. Adams and Wake^ have presented an excellent 
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review of the analyses, experiments and designs of these joints carried out by 
several researchers. Review of relevant literature carried out in this Chapter, 
includes (i) linear elastic analysis, (ii) material (elasto-plastic and elasto-visco- 
plastic) and geometrically non-linear analysis, (iii) fracture mechanics, (interface 
failure, fracture, delamination etc.), (iv) parametric study, (v) environmental 
effects, (vi) testing of bond quality and (vii) design of adhesively bonded joints 
(single lap, double lap, butt, scarf, tubular etc ). The Chapter also summerises 
the various findings so far, the lacunae and the motivation to carry out the 
present research work. 

In Chapter 2, the linear elastic finite element analysis of the joints, the 
modelling of crack opening is presented. Automatic mesh generation is adapted 
for the analysis of the joint. Chebyshev’s function is employed to refine the 
mesh at the critical zones. Stress distribution in the adhesive across the overlap 
of single lap joint^ is examined in detail. Stresses and strains distributions 
across the overlap and through the adhesive thickness of stepped lap joint® 
are also looked into. These results shows good agreement with corresponding 
results from the literature^*®. Crack is modelled with paired nodes. Employing 
Irwin’s crack closure technique strain energy release rates are computed. A 
single edge notch (SEN) specimen is considered to determine the minimum 
refinement required to obtain reasonably accurate (close to exact value) strain 
energy release rate. Further, a cracked lap shear (CLS) specimen® is considered 
to validate the finite element code computing strain energy release rates. 

Parametric study of the double lap joint is also presented in this Chapter. 
Seven possible fracture criteria are considered and the corresponding predicted 
failure loads are compared. Strengths predicted with the criterion {Gj + Gjj > 
Gjc + Gjjc) are found to be closer to those obtained from the experiments 
than with the rest of the criteria. Hence, this criterion is considered in the rest 
of the analyses. Effect of use of the same and different adherend materials 
(sti&ess) on joint strength are examined. Effect of equal and unequal adherends 
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(thickness) on joint strength is also investigated Response of the joint of 
metal to metal, composite to metal and composite to composite is discussed 
Effects of various parameters on the strength of composite to composite joint 
have been looked mto These parameters include stacking sequence, boundary 
conditions, crack location, adhesive material properties, overlap length and 
adhesive thickness on joint strength are studied Some of the conclusions drav/n 
on the basis of this study are 


• Use of Chebyshev’s function is found to be effective for automatic mesh 
generation with vanous refinements at critical zones 

• In case of SEN specimen, with 2688 elements and proper mesh refinement 
near the crack-tip, error estimated in the value of fracture parameter is 
less than 1 6 percent 

• Result based on two dimensional finite element analysis with fracture 
mechanics approach agreed well with the expenmental result for CLS 
specimen and the difference is less than 2 7 percent 

• Introduction of 0° ply m the immediate neighbourhood of adhesive results 
m a stronger joint 


Elasto-plastic analysis of the adhesively bonded joint is presented m 
Chapter 3 The finite element model, flow chart and the algorithm for elasto- 
plastic analysis are also included Finite element code is validated on the 
cantilever beam problem^ Spread of plastic zones based on elasto-plastic 
analysis for vanous boundary conditions and crack location are presented 
Effe^ of stacking sequence, adhesive material properties, overlap length and 
adhesive thickness on the joint strength are investigated 
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• The elastic analysis predicts a lower failure load as compared to the elasto- 
plastic analysis This might be due to the peak stresses of elasto-plastic 
analysis are reduced 

Visco-plastic analysis of the adhesive joint is presented in Chapter 4 Finite 
element model, flow chart and the algonthm employed for the analysis are 
included Finite element code developed is validated on the perforated plate 
problem® The visco-plastic stress analysis results show good agreement with 
those available in literature for a single lap joint and a stepped lap joint Effects 
of stacfang sequence, fracture toughness of the adhesive, overlap length and 
thickness of the adhesive on the joint strength are examined 

• The elastic analysis predicts a lower failure load as compared to the visco- 
plastic analysis This might be due to the peak stresses of visco-plastic 
analysis are reduced 

The evaluation of fracture properties requires the knowledge of critical 
strain energy release rates (fracture toughnesses) in mode I and II In Chapter 5, 
some simple expenments are performed to evaluate these values are discussed 
Joint strengths are measured from these expenments and are observed to be 
close to those obtamed from the linear elastic finite element analysis Vanous 
possible fracture critena are considered Among all possible fracture cnteria, 
{Gj -t- Gii > Gjc + Gjic) IS found to be most appropriate 

The thesis concludes with Chapter 6 Results and discussion of the present 
work are presented in this Chapter Some of the general conclusions drawn on 
the basis of the present work are 

• The finite element approach with paired nodes modelhng the crack provides 
a simple computational tool to predict joint failure Howevct, flhc gnds 
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near the crack-tip are necessary for accurate estimation of strain energy 
release rates 

• High Young’s modulus of the adhesive, large overlap length and small 
adhesive thickness individually results in stronger bonded joint 

• A proper choice of layup sequence results in an efficient bonded joint 

• The computed bond strengths of double lap joint based on finite element 
analysis with fracture mechanics approach and those determined from the 
expenment show good comparison 

• The fimte element code developed predicts linear and non-linear behaviour 
of the adhesive lap joints fairly accurate^’® 

• The elastic analysis predicts a lower failure load as compared to the elasto- 
plastic and elasto-visco-plastic analyses This might be due to the peak 
stresses of non-lmear analysis are reduced 

• The results obtained from parametnc study are similar for linear elastic, 
elasto-plastic and elasto-visco-plastic analyses 


Limitations of the present analyses are also indicated Finally few suggestions 
are made for further extension of the present work 
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Chapter 1 
Introduction 


1.1 Motivation 


Joints are inevitable in structural applications Though, they constitute a small 
fraction of the total weight of a structure, they are important as load transfer 
and stress difihision takes place through these joints Therefore, a good design 
could result in better utilisation of bulk matenals Adhesively bonded joints are 
being preferred to mechanical fasteners because of distinct advantages such as 


• uniform stress distribution 

• stiff connection, 

• smooth surface contour, 

• sealed against corrosion, 

• excellent fatigue properties, 

• no fretting problem, 

• relatively light weight. 


1 



2 


CHAPTER 1 INTRODUCTION 


• damage tolerance, 

• stealthiness etc 


An adhesively bonded joint and various mechanical fasteners, which could 
carry same amount of load (15 kN) are shown in Figure 1 1 Outer members of 
these timber joints, were 3| X 1 mch and the inner 3| X 2 inch These figures 
clearly shows that adhesive bonding saves sigmficant amount of matenals and 
m addition has same strength as the mechanical fasteners 

With the development of high performance structural adhesives, these joints 
are being employed in a wide variety of structural applications ranging fi:om 
aerospace, aircraft, automobile components to those found in microelectronic 
components In the aerospace industry extensive applications of adhesively 
bonded joints are found These include, honeycomb sandwich floor panel, 
rudder floor panel, engine cowling, elevator, flap, ceramic tyle of re-entry 
vehicle, integrated fuel tank, composite propeller blade, helicopter blade, launch 
vehicle attachment (LVA) ring, collapsible tube mast etc 

A group of researchers are working in the primary adhesively bonded 
structures technology (PABST) programine*^with the objective to develop a 
light transport aircraft, whose most parts will be of composite materials and 
adhesively bonded joints This is expected to result in less maintenance and 
operational costs and also the required runway length will be significantly short 
as compared to the present ones 


Non-aerospace applications of adhesive bonds includes, Rolls-Royce carbon 
fibre reinforced plastic (CFRP) cowl door, storage barrel, lamp post, tire 
reinforcement, abrasive wheel, brake lining, safety glass, William Honda grand 
pnx car shaft, decking, furnitures, magnets, yacht etc ' - . ^ 
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Composite materials arc now used extensively in a wide range of applications, 
because of specific properties like (i) high specific strength and stiffness, (ii) 
dimensional stability, (in) electrical and thermal insulation, (iv) no corrosion, (v) 
smooth outer surface, and (vi) material could be tailored to suit the application 
Optimal use of the bulk material and desired redistribution of material properties 
result in ’engineered composite structure’ 

Adhesively bonded joints are also found to be more suitable with fibre 
reinforced plastic (FRP) composites because of the facts that conventional 
welding, brazing methods are not applicable for these composites Moreover, 
the adhesives used in the bonding are similar m nature as the basic matrix used 
for the reinforcement of the composites 

There are large number of parameters, which affect the performance of 
the adhesively bonded joints of fibre reinforced plastics Following parameters 

• relative stiffiiesses of laminates (adherends) 

• relative stiffnesses of adhesive and laminate (composite adherend) 

• relative thicknesses of laminates (adherends) 

• relative thicknesses of adhesive and laminate (adherend) 

• overlap length 

• employment of various types of joints such as 

smgle lap, double lap, stepped lap, butt, scarf, tubular etc 

• displacement boundary conditions 

• stacking sequence of layups in the laminates (adherends) 

• surface preparation, use of pnmers, mechanical properties of the adhesive 

matenal, and shape of the adherends etc 
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are found to have strong influence on the performance of these joints 
Adhesives used in these joints are usually of polymeric materials, which exhibit 
non-linear material behaviour such as elasto-plastic and visco-plastic behaviour 
under high stress and/or high temperature Moreover, in adhesive bondings, 
mainly two types of failures are observed viz (i) cohesive type damage 
occurring in adhesive matenal, and (ii) adhesive type damage occurring at 
the interface of adhesive and adherend 

The motivation for present study is, therefore, to study the response and 
fracture of the adhesively bonded jomts of composite matenals, in the presence 
of material non-lineanties It has been observed that there is large amount 
of scatter in the expenmental data Similar observations are also found in the 
theoretical results Hence, in order to confirm that the theoretical results are 
realistic, experiments need to be conducted on well designed specimens 


1.2 Review of literature 


1.2.1 Elastic analysis 


I Lap joints Figures 1 2a and 1 2b 

Shear lag analysis of adhesively bonded joints was carried out by Volkerson 
(1938) He observed that absolute maximum of shear stress was independent of 
overlap length However, bending moment caused by non-colinear forces was 
not accounted Further, adherends bending allowing the joint to rotate gave rise 
to geometrical non-linearity was also neglected Goland and Reissner (1944) 
determined stress distribution m the adhesively bonded joint under tbnsile load 
Th^ also obtained solutions for thin W thibk demertt layers In this analysis. 
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the effect of flexibility of the joint was considered only for thick cement layer 
However, shear deformation and tensile stresses across the adherends were 
neglected Further, the joint was considered to be monolithic (zero glue line 
thickness) 

Chen and Cheng (1983) obtained a closed-form solution, which was based 
on two-dimensional elasticity theory in conjunction with the variational principle 
of complementary energy and was adaptable for various possible adhesive layer 
flexibility Using an impact test. Hams and Adams (1985) measured strength 
and energy absorption of bonded single lap joints with four different epoxy 
adhesives and three different aluminium alloy adherends Static joint strength 
was found to be not significantly affected by loading rate Energy absorption 
observed to be large when joint strength was significant for plastic deformation 
to occur in the adherends pnor to failure From the crush tests of bonded 
and spot welded q^linders, they found that bonded joints were supenor in 
performance as compared to spot welded ones Adams and Harris (1987) 
investigated the effect of local geometry at the edges of the overlap m single 
lap joints Three types of joints viz (i) joints with square edged adhesive 
layer, (ii) joints with fillet of adhesive and (iii) joint with an adhesive fillet 
and tapered adherend for jomt strength were analysed They found that finite 
element analysis was capable of predicting the significant strength increase, 
which was achieved m single lap joints by filleting the adhesive at the edges of 
the overlap and rounding the ends of the adherends 

Pindera and Wang (1990) obtamed the three dimensional stress distnbution 
in a symmetnc adhesively bonded joint using techniques of analytical and optical 
isodynes Cheng, et ol (1991) studied the behaviour of the adhesively bonded 
single lap joints with non-identical adherends under the action of tensile forces 
Their results showed that the intensities of the stresses were ^ways much greater 
in a small zone at both ends of overlap region of the joint ARALL laramates 
are q family of structural composite materials, wfiidh consist of^thm alunp^nium 
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alloy sheets alternating with aramid fibre/epoj^ prepreg These laminates are 
developed for layers developed for fatigue applications Using finite element 
modelling and shear lag analysis failure of the adhesive m ARALL was predicted 
by Long (1991) He inferred that the failure load was a function of joint 
configuration and overlap length 

Dorn and Weiping (1993) studied the influences of spew fillets, different 
adhesives and metal adherends on the stress and strain states in adhesively 
bonded plastie/metal single lap joints the strength and failure properties of 
these structures were undertaken by finite element technique (ADINA software) 
and experimental approaches They observed that the introduction of spew 
fillets lead to significant reduction m peak values of adhesive shear and peel 
stresses Further, the use of aluminium adherend instead of steel in plastic/metal 
bonded joint could decrease the adhesive stiffness imbalance and hence improve 
the symmetry of adhesive stress distnbution They also found that the failure 
occurred in the plastic adherend in the spew end region and suggested non-lmear 
analysis for the accurate prediction of joint strength 


I Butt joint Figure 1 2c 


Using linear elastic finite element method, adhesive butt joints were investigated 
by Alwax and Nagaraja (1976) They concluded that there were significant 
variation m stress pattern across the adhesive thickness particularly near the 
adhesive-air interface Further, they observed that, there were surface layers 
at the interface, whose charactenstics were different from the bulk adhesive 
material However, this work did not include the influences of these surface 
layers and the viscous behaviour of these layers Sawa, ct al (1989b) extended 
them earifer work to include three dmiensional stress analysis of the adhesive 
iHrtt joints Ibimna. et^oll (19^#a> fh^^hree «dnal stress analysis 



1 2 Review of literature 


7 


of adhesive butt joints subjected to tensile loads The effect of debonded area 
and a spew fillet on the principal stress distributions were numerically calculated 
It was shown that, the joint strength decreased considerably when the joint had 
a debonded area and got increased when the joint had a spew fillet 

Temma, et al (1990b) carried out two dimensional stress analysis of 
adhesive butt joints, in which two dissimilar finite strips were joined by an 
adhesive subjected to cleavage loads Effects of elastic moduli ratio of adherend 
and adhesive, Poisson’s ratio, thickness of the adhesive and the cleavage load 
on the stress distribution at the interface were studied They measured strains 
produced at the adherends with strain gauges and stresses by photoelastic 
techmques These results were comparable with those obtained from numencal 
methods Temma, et al (1991) examined the stress distnbution of an adhesive 
butt joint having a circular hole in the adhesive when the joint was subjected 
to an external bending moment and also the influence of the hole size on the 
strength of the joint They also conducted photoelastic expenments to measure 
the stress distnbution in the adhesive 


I Scarf joint Figure 1 2e 

Stress distributions, m a plane scarf joint with arbitrary angle of scarf and 
vanous matenal (elastic) parameters associated with adherends and adhesives, 
were analysed by Chen and Cheng (1990) The analysis was based on the two 
dimensional elastic theory m conjunction with the principle of complementary 
energy Some typical features of the stress distnbution were shown Usmg 
elastic finite element analysis, Ikegami, et al (1990) investigated the strength 
of adhesively bonded scarf joints comprised of an epoxy adhesive and two 
adherends namely glass cloth laminated plastics and metal They also examined 
the effects ot scarf length and elastic imbalance on the deformation and strength 
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of the joint 

Marcolefas, et al (1991) suggested a method for the analysis of adhesive 
bonded metal to composite scarf joint, in which both the metal and the adhesive 
exhibited non-lmear behaviour A failure cntenon was proposed for this type 
of scarf joint and for a given joint this method calculated the bond strength 
as a function of overlap length Further, they observed that thermal mismatch 
of adherends and adhesive had a severe effect for intermediate overlaps only 
However, even with elasto-plastic adhesive behaviour reliability of this analysis 
was limi ted, because of complexities due to non-identical adherends and non- 
umformities in the adhesive shear stress distnbution 


I Tubular pint Figure 1 2f 

Using interface finite elements Habib, et al (1988) earned out stress analysis in 
axisymmetnc structures, which took into account the interface conditions due to 
geometneal continuity and interface equilibrium Further, stress distnbution m 
tubular lap joints under internal pressure and mfluence of the matenal properties 
of the adherends on the mterfacial stress distribution were investigated They 
observed significant vanation of the stresses across the adhesive thickness, 
specially near the two edges of the joint The stress variation was observed to 
be non-uniform across the thickness Further, the edge effect was significant 
up to a distance of 15 times the thickness of the adhesive from the edge 
Imanaka, et al (1989) conducted fatigue tests and also investigated analytically 
adhesive bonded joints in shafts. They found that the maximum normal and 
shear stresses are the key parameters for estimating fatigue strength under 
pulsating tensile and torsional loading conditions, respectively These maximum 
stresses under tensile or torsional load conditions increased with the increase 
Quilling (oldriafiy Idiigth Fui-ther, it was observed that shorter cddpling 
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lengths resulted in a higher fatigue life and proper tapenng in the coupling 
zone improved the fatigue strength 

Harj-Rabab, et al (1990) employed rectangular and axisymmetric interface 
finite elements for the analysis of adhesive bonded joints The stress distributions 
were compared with the results of closed-form solution and of displacement finite 
elements The influence of defects on the stress distribution was also studied 
Though this method provided accurate stress distribution at the interface and 
solved local effect problems, the drawbacks of these elements were (i) a large 
number of degrees of freedom per node were needed, (ii) difficulties to solve 
resulting equations and (iii) computation time was more Further, this analysis 
was limited to linear elastic behaviour of the constituents of the bonded joint 

Sawa, et al (1992) investigated the strength of the joints combining 
adhesives with bolts The stress distribution m the adhesives and a vanation of 
axial bolt force were analysed usmg the theory of elasticity for axisymmetncal 
structures They found that the strength of the combined joint was greater than 
that of conventional adhesive joints, and the strength increases with decrease 
m the bolt pitch diameter They also performed expenments to determme 
joint strength, maximum stress m blots for vanous loading conditions Sawa 
and Kobayashi (1988) descnbed the strength charactenstics of a combined joint 
consisting of an adhesive with bolts Hollow <ylinders were considered as 
adherends and stress distnbution at the boundary of adhesive layer, the force 
ratio and the strength of the joints were examined They found that the strength 
of the combination joint was greater than that of an purely adhesive jomt 

Using finite element technique Davies and Khalil (1990) analysed double 
containment joint made up of steel and epoxy, and based on overaU stiffiiess of 
the corner joint designs were recommended With modest increase m elastic 
modulus of adhesive^ they observed substantial increase m joint stiflEoess Chen 
and Cheng (1992) analysed the stress distribution m adhesive-bonded tubular 
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lap joints subjected to torsion The analysis was based on the elasticity theory in 
conjunction with the variational principle of complementary energy Adherends 
with different thicknesses and different adhesive materials were considered in 
the analysis 


1.2.2 Non-lmear (material) analysis 
I Elasto-plastic analysis 

By means of elastic-plastic analj^ical technique, adhesively bonded double 
lap jomts were analysed by Hart-Simth (1973) On the basis of the study, 
he concluded that (i) there was a definite limit to the bond strength that 
could develop between specified adherends and a particular adhesive (ii) the 
maximum potential bond shear strength was determined by the adhesive strain 
energy i c shear stress per umt bond area (m) adhesive peel stresses were more 
critical at the ends of the overlap than those associated with shear stresses He 
also studied the effect of thickness on the strength of double lap joints Further, 
Hart-Smith (1981) analysed the adhesive-bonded joints using the continuum 
mechanics approach and concluded that (i) adhesive shear stress distribution 
was inherently non-uniform for single lap, double lap, stepped (Figure 1 2d) 
and scarf (Figure 1 2e) bonded joints (ii) the inclusion of adhesive non-linear 
(el^to-plastic) behaviour was vital if the bonded joints analyses were to predict 
failure (iii) very little portion of a typical bonded joint was actually loaded 
beyond the adhesive elastic compatibility, and plastic behaviour was confined 
to the immediate vicinity of the ends of the bonded overlap or to the ends of 
the steps in a stepped lap joint (iv) a precise representation of the adhesive 
non-Jmear bdjaviour was unnecessary, any model which was mathematically 
trieteble could (v) adh^ive peel stresses, which occurred in conjunction 
slr^^es could' be tdlfevltded by suitable tapermg of the adherend 
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or local thickening of the adhesive (vi) the strength of bonded joints of the 
simplest configuration was limited, so the use of adhesive bonding in thicker 
structure demands more complex geometries of the joints and (vii) bonded 
joints should always be designed in such a way that the adhesive was never the 
weak link 

Ramamurthy and Rao (1984) analysed lap joints with elasto-plastic adhesive 
behaviour and uniform adherends This analysis also accounted for arbitrary 
adherend shapings They concluded that strength of a bonded joint depended 
on adhesive behaviour and the strength of adherend materials However, this 
work was limited to isotropic matenals and did not consider the effects of 
other important parameters like adhesive thickness and overlap length etc on 
bond strength Czarnocki and Piekarski (1986) observed that the failure of an 
adhesive layer in a symmetric lap joint originated in the central plane of a joint 
at the front edge of an adhesive layer and was caused by the biaxial state of 
normal stresses Further, the most severe state of stresses existmg m the central 
plane of a joint was not effected by joint width (within the range investigated) 
The load capacity of a joint per unit width was found to be constant They also 
found that use of weaker but sufficiently flexible adhesives resulted m a higher 
load capacity than that obtained with a strong but rigid adhesive, if the length 
of the jomt was not restricted Further, application of a more flexible adhesive 
could significantly reduce stress concentration in adherends However, matenal 
non-hneanty was not accounted in the analysis 

Using non-hnear (elasto-plastic) finite element analysis of both clad and 
unclad configurations Crocombe and Evan (1988) investigated the effect of 
the cladding layer on the strength of a bonded joints They found that for 
the unclad joints yielding occurred first at the tensile zone and spreads frcmi 
the overlap region Analysis of clad jomts showed a similar yielding in the 
overlap region but the adhesive fillet were stresses to a much higher level by 
the excessive deformation of the cladding layer Adams (1989) inferred that 
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finite element technique was one of the best way of treating the non-Iinear 
(joint rotation) mechanics, material properties, adhesive fillet and tapenng of 
the corner in bonded joints He showed that the strength of the joint could be 
increased up to 500 percent by using composite adherends effectively Bassani 
and Qu (1989) calculated the crack-tip stress fields on the interface between 
two semi-infinite anisotropic elastic solids They estimated small scale yielding 
plastic zones surrounding the crack-tip They also found that crack-tip stress 
had the standard square-root singularity and the corresponding stress intensity 
factors were separable m the three modes of loading 


Bigwood and Crocombe (1989a) earned out elastic and elasto-plastic 
analysis of general jomt configurations and suggested a two parameter formula 
for preliminary design Lee, et al (1989) developed computer programs for 
analysis of adhesive joints in co-anal tubular structure This program took 


mto account plasticity in adhesive and predicted the stress m a lap shear joint 
subjected to any combination of tensile, shear, bending and thermal loads 
They conducted tests to determine adhesive properties and the shear strength, 
stiffness and durability of structural joints Groth and Brottare (1989) earned 
out a fimte element analysis for determination of the joint stiffness as function 
of the nominal strain Adhesive matenal was considered to exhibit elasto-plastic 
behaviour They observed that apparent stiffness of the adhesive layer kept 
on decreasing due to plasticity at the interface and saturates below a nominal 


strain level of 05 percent However, this work was limited to thick adhesives 
Additionally, Bigwood and Crocombe (1990) earned out the elasto-plastic finite 
difference analysis of adhesively bonded joint considering as a plane strain 
problem and compared these results with those obtained from elasto-plastic 
finite element analysis They assumed that the material behaviour of adherends 
and adhesive was linearly elastic and elasto-plastic, respectively The joint 
modelled a.s an adherend-adhesive sandwich allowing the application of 

Z"TT '' the 

adherends > They- also studied the dfect of adhesive thickness oa the strength 
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of lap joint Schemes for estimation of J integral for predicting elastic-plastic 
fracture was explained by Brust, ct ol (1990) Edlund and Warbrmg (1990) 
carried out two-dimensional, elastic-plastic finite element analysis of a shear 
loaded adhesive joint Crocombe and Bigwood (1992) outlined an approach for 
including the adherend non-lmearity into the adhesively bonded joint analysis A 
bi-linear elasto-plastic model incorporating strain hardening was used to model 
the adherend stress-strain behaviour However, this approach was complicated 
due to the complexities of the governing equations, involved in the analysis 

Thomsen (1992) considered the problem of an adhesively bonded lap joint 
between two dissimilar orthotropic circular cylmdncal laminated shells, which 
includes elasto-plastic behaviour of the adhesive He observed that the overlap 
length and stiffnesses of adhesive and adherends significantly influence the 
the stress distribution in the joint He also showed that with elasto-plastic 
analysis considerable reduction in peak stresses were achieved as compared 
with purely elastic analysis Adams and Mallick (1993) presented a method to 
predict the strengths of elasto-plastic adhesives at room temperature For most 
cases, predictions were within 5 percent variation from the experimental values 
However, these predictions were not accurate for the cases of brittle adhesives 
at low temperatures (— 55°C) and high overlap lengths 


I V5sco-eJasfrc and visco-plashc analyses 

Creep behaviour of structural adhesives used m lap joints with steel adherends 
was studied by Allen and Shanahan (1975) Adhesives employed were (i) a 
phenolic resin/polyvinyl formal composite Redux 775 and (u) a modified eposy 
novolac Redux BSL 906 It was shown that the separation of adherends due 
to adhesive shear creep should be associated with the maximum adhesive shear 
stresses found at ,the eqds pf the overlap. Due tp^tlie crpep m adhestve bonded 
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lap joints, Allen and Shanahan (1976) observed reduction in stress concentrations 
at the ends of the overlap as compared with elastic solution They also found 
that creep curves obtained from the lap joints did not correspond directly 
to those obtained from experiments employing bulk adhesive material under 
similar loading conditions Further, they concluded that if creep was allowed, 
the long joint had a secondary safety factor in that stresses causing creep were 
reduced rapidly However, there was considerable amount of scatter reported 
in the expenmental results 


Devitt, et al (1980) suggested an approach for the determination of 
delamination fracture toughness in the opening mode, which was developed 
usmg a split beam and a non-linear visco-elastic analysis They determined 
the visco-elastic behaviour of the resin from creep/recovery tests They also 
showed that the measured range of energy release rates for a range of crack 
growth rates, was consistent with predictions from an idealised visco-elastic 
crack growth theory Using finite element program (SAAS III) Sen and Jones 
(1980a) determined the stresses and strains across the adhesive thickness in a 
double lap bonded joint, subjected to a quasi-static load Visco-elastic stress 
analysis was performed with the finite element program and Schapeiy’s direct 
method of transform inversion Nagaraja and Alwar (1980) carried out finite 
element analysis of adhesively bonded plate lap joint considenng the visco-elastic 
behaviour of the adhesive matenal They found that the decrease of stresses 
with time was of the order of 18 percent for normal stresses and 45 percent 
for shear stress Further, the vanation of stresses across the adhesive thickness 
was confined only to the regions near the edges This method was also found 
applicable to study the response of solid propellants in rockets 


Ratwani, et al (1982) disdlfesed the importance of consideration of visco- 
ete of adhta.ve for strength prediction of the bonded joint They 

fee jnaximum* shear stn^ m the adhesive layer of a stepped 
<^adrl6^|ladeaflyh^^the^#irie dependent behaviour was 
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taken into consideration and the shear distnbution becomes fairly uniform after 
about 2 hours Henriksen (1984) developed an algorithm for the analysis of 
non-linear visco-elastic adhesive matenals This formulation also took into 
account high temperature, cyclic creep, finite and infinitesimal strains in adhesive 
joints Based on J integral, Schapeiy (1984) developed a method of quasi- 
static deformation and carried out fracture analysis for a non-lmear visco-elastic 
matenal Their analysis, included crack growth, crack closing and crack healing 
and, was applicable to problems involving ablation, interfacial contact and 
interfacial separation Groth (1986) presented a general method to calculate the 
stresses in bonded joints using the substructunng technique, where large parts 
of structure were linearly elastic and remaining part exhibited non-linearity 
He also pomted out that visco-elasticity of the adhesive was very important 
because large redistributions of the stresses occurred during the penod when 
the joint was loaded He also concluded that normal stresses changed rapidly 
at the terminus and more number of elements were required there in order to 
evaluate the stresses accurately Yadagin, et al (1987) studied the visco-elastic 
analysis of adhesively bonded joints using finite element method with curved 
isoparametric elements Hereditary integrals used to represent the stress-strain 
relations resulted m a stiff model while the use of relaxation modulus resulted 
m a flexible model They found that the average of these two solutions resulted 
m a better approximation 

Roy and Reddy (1988a) developed a general purpose finite element analysis 
software called NOVA to analyse bonded joint, which takes mto account 
the visco-elastic adhesive behaviour and the geometrical non-hneanty The 
effect of temperature and stress levels on the visco-elastic response was taken 
mto account by a non-lmear shift factor definition Penetration sorption was 
accounted by non-lmear Fickean diffusion modfi m which diffusion coefficient 
was dependent on the temperature, penetrant concentration and the dilatational 
stram The adhesive layer was modelled usmg Schapeiy’s non-lmear smgle 
integral constitutive law for uniaxial and multiaxial state of stress by Roy and 
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Reddy (1988b) The effects of temperature, stress level and penetrant absorption 
on the visco-elastic response were studied However, the effect of change m 
Poisson’s ratio with time in the polymer was not taken into account, which had 
a significant bearing on the final response of the joint 

Needleman (1990) discussed vanous issues associated with the analysis of 
interfacial failure phenomena and inferred that visco-plastic material behaviour 
was responsible for limiting the achievable stress level He also indicated that 
a key parameter controlling decohesion was the ratio of the strength of the 
interface to the strength of the ductile solid Chalkley and Chiu (1993) carried 
out experiments to study the shear stress-strain behaviour of FM 73 epojg^ 
They found that FM 73 exhibited a more prominent visco-plastic behaviour 
than visco-elastic behaviour Based on the results of thick-adherend and short- 
overlap specimens, it was concluded that due to the time-dependent nature of 
the adhesive, true shear stress/stram testing could be accomphshed only under 
stram control But, m practice thin adherends and long overlap lengths were 
employed to achieve higher strength 

Fujii (1993) presented a numencal model using a one-dimensional approxim- 
ation for analysing the dynamic response of adhesively bonded sandwich beams 
The visco-elastic characteristics of the adhesive damping layers were represented 
by generalised Maxwell model To determine the dynamic response, they 
proposed a new method which combines the transient analysis based on finite 
element method with the expenmental modal analysis This approach was 
validated with an impact test on a sandwich beam bonded with chloroprene 
rubber adhesive damping layer Hiregoudar (1993) carried out finite element 
analysis of single lap joint This analysis included the geometrical and material 
(visco-plastic) non-linearities However the analysis was limited to ideal bonded 
|oint, having no defects at All Su and Mackie (1993) earned out two-dimensional 
creep analjrsis of structural adhesive joints They observed significant amount of 
#Cous ddfbrmation even at room tmperature. ? At elevated temperature (45®C) 
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the Young’s modulus was observed to be reduced substantially as compared 
with that at room temperature Both normal stresses and shear stress values 
got reduced because of the creep in the adhesive joints The reduction of shear 
stress was rather limited Further, with creep there was an increase in strains, 
and maximum increase was observed in shear strain 


1.2.3 Non-linear (geometrical) analysis 


Results of stress analysis of adhesively bonded joint and calculation of strain 
energy release rates for cracked lap shear (CLS) specimens at four different 
debond lengths were reported by Johnson (1987) The various numencal 
techniques considered were two- and three- dimensional finite element analyses, 
beam theory, plate theory, and a combination of beam theory The results were 
clustered into two groups, geometncally hnear and geometncally non-linear 
analyses He suggested that care was needed to ensure proper modelling of 
the specimen in order to get correct result eg adhesive bondline must be 
modelled with proper grid density and full bondhne and not just near crack-tip 
region, needed to be mcluded in the model Further, he pointed out that, 
there was a small amount of contnbution from out of plane shear deformation 
(mode III) of adherends Reddy and Roy (1988) developed a two-dimensional, 
finite element code using an updated Lagrangian formulation to account for 
the geometric non-lmearity for the analysis of adhesively bonded joints. The 
effects of boundary conditions and mesh on the stress distnbutions in lap joints 
were also investigated It was observed that the ^e of sboundary conditions 
assumed cbuld siipaSeafitly alter die stress distnbution iii ^ joint 
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1.2.4 Fracture mechanics approach 


Using beam theory and J integral, Femlund and Spelt (1991a) developed an 
numerical solution to calculate the strain energy release rate of adhesive layer, 
which included both geometrical and material non-lmearities They observed 
that (i) double cantilever beam showed a predominant mode I fracture, (ii) end 
notch fracture (ENF) specimen produces almost pure mode II, and (iii) the 
cracked lap shear (CLS) specimen joint showed a mixed mode fracture Using 
modified crack closure integral and constant strain finite element stress analysis 
(AXISOL software), Rybicki and Kanmnen (1977) evaluated stress intensity 
factors and applied to three configurations, (i) a double cantilever beam test 
specimen, (ii) a finite width plate containing an embedded crack and (ni) a 
bolt fastened double lap joint containing two radial cracks They observed that 
the results obtained with the use of higher order elements were less accurate 
as compared to those obtained with constant strain finite elements, which uses 
nodal forces and nodal displacements However, their work was limited to 
linear elastic analysis and isotropic materials 


Using the finite element computer code NASTRAN, Ratwani and Kan 
(1986) earned out the analysis of delammation propagation m the composite 
layer of metal to composite stepped lap joints They found that the delammations 
did not grow to be catastrophic in the composite layer for some resin material 
■:^e presence of both debondmg and step-end gaps in the adhesive layer 
r^resented a critical condition from the point of view of delammation initiation 
and propagation Janardhana, et al (1986) employed finite element model to 
^udy ^e effwt of debond length on various lap joints They also examined 
n, f m ‘l^Jammation damage in composite to metal joints and conducted 

graphrte/epoxy and aluminium 
^98$) Wrwd ^ut ei^erimental and analytical investigations 
loMel to Ihafdctense 
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debond growth mechanisms under mode I and mixed-mode (I and II) loadings 
Failure occurred in the form of debond growth in all the specimens tested Total 
strain energy release rate appeared to be the driving parameter for cohesive 
debond growth under static and fatigue loadings Also the values of total strain 
energy release rates measured for debond growth were an order of magnitude 
below the static fracture toughness (critical strain energy release rate) 

Groth (1988a) predicted fracture load of the single lap joint with failure 
criterion, which was based on matenal-mduced singularities appeared at the 
terminus of the joint These results were compared with those obtained 
from the experiments Pilarski and Rose (1988) used an ultrasonic oblique- 
incidence technique to detect the interface quality in alumimum to aluminium 
adhesively bonded structure However, these measurements were sensitive to 
frequency range used, angle of incidence, attenuation and reflection coefiicient 
Aivazzadeh, et al (1988) presented several tnangular and rectangular mixed 
finite elements for the interface stress analysis m adhesive butt joints For two- 
dimensional interface problems such as sandwich structures and laminates they 
employed these elements and evaluated stress distnbutions in these structures 
However, these elements had the drawbacks of mixed elements and also it was 
difficult to extend this to non-lmear analysis 

Aivazzadeh and Verchery (1986) presented rectangular interface finite 
elements for the stress analysis of adhesive joints The results obtained 
from these elements showed that, stress distribution could be evaluated more 
accurately than those obtained from conventional displacement elements However, 
this work was limited to two-dimensional linear elastic and rectilinear interfaces 
Hatton, et al (1989) suggested a method which uses two stress singulanty 
parameters for evaluating adhesive strength and applied to several kmds of 
moulded models, composed of epoxy resin and metal (Fe-Ni) alloy sheets, 
and plastic encapsulated large scale integration (LSI) models The initiation 
and extension of delamination were estimated theoretically Occurrence of 
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delammation and its propagation in the molded models were observed using 
scanning acoustic tomography Crompton and Clark (1989) examined the near 
crack tip stresses with finite element analysis and stress intensity calibration 
equations and found that for a 0 5 mm bond line thickness the stress singulanty 
deviates from the conventional square root singularity 


Groth and Brottare (1988) proposed a fracture parameter namely singular 
intensity factor equivalent of stress intensity factor (SIF) for the failure analysis 
of butt joint with various adhesive materials Sethuraman and Maiti (1989) 
developed an element with a constant shear strain field for the analysis of 
doubly bonded crack-stiffened panels serving as crack arresters They found 
that the stiffener accrued more benefits m the case of an edge crack than in the 
case of a centre crack, when both were subjected to the same mode I loading 
Further, increasing the shear modulus or reducing the thickness had a beneficial 
influence on the stress intensity factor of the stiffened structure Moreover, 
the stiffener was observed to be most effective when it was located close to the 
crack tip ^n on the side opposite to the ligament Crocombe, et al (1990) 
outlined various approaches applied in analyses of adhesively bonded joints to 
predict the strength of these joints The applications of the singularities at joint 
bimatenal interfaces were discussed and failure criteria were proposed 


Kairouz and Matthews (1990) postulated various mechanisms of failure in 
bonded CFRP single lap joints These mechanisms were found to be dependent 
on the layup and stacking sequence of the adherends From the joints tested, 
strength of the joint was found not to be strongly dependent on stacking 
sequence, but failure mechanism was influenced by the direction of surface 
fibres and was determined by stresses induced by bending They concluded that 
for adherends with 0° ply/phes on Ae surface failure initiated as tensile fracture 
and was followed by delammation within 0° ply Further, they 

“ " ' proceeded inwards 

“ by rntralkmmar and mterlaminar 




^ I * 


^ i 





1 2 Review of literature 


21 


shear cracking A numerical study of the coin-tap test was carried out by 
Mackie and Vardy (1990) Using a two-dimensional finite difference model 
of elastic waves in an adhesively bonded structures, they found a relationship 
between defect size and the frequency of vibration of the material However, 
this technique was effective for large civil structures only 

Mecklenburg, et ol (1990) applied the energy separation method for 
characterising the fracture resistance of adhesives and compared it with elastic 
and plastic fracture parameters such as strain energy release rate and J integrals 
of bonded half compact tension test specimens Suo (1990) discussed several 
mechanics problems of brittle adhesive joint in fracture testing of nominal 
toughness, configuration stabihty of mode-I in-layer fracture, tunnel crack 
and parallel debond due to residual stresses He concluded that parallel 
debonding was likely to take place for adhesives in residual compression than m 
tension, provided other conditions were the same Khalil and Bayoumi (1991) 
determined the fracture energy and fracture toughness of an epoxy adhesive 
matenal at different loading speeds with cleavage test specimens They observed 
that with increased loading speed the effect of loading rate on the fracture 
energy and toughness decreased Correlations between the fracture energy and 
fracture toughness were obtained for both aluminium and brass adherends 
Liechti and Chai (1991) developed a biaxial loading device, which produced 
steady crack propagation under all loading directions Energy release rates 
were extracted from linear elastic finite element solutions that matched with 
the measured normal crack opening displacements (NCOD) Zenkert (1991) 
mvestigated the debondings m foam core sandwich beams with cracks m the 
interface between the face and core Using fracture toughness criterion, he 
predicted the fracture loads for beams with simulated debondmgs subjected to 
four pomt bendmg' He also inferred that interface debondings had significant 
influence on the load bearing capacity of foam core sandwich beams . 


Ishtkaiya, et (1991) mvesti^atdd the mixed mode fractqre criterion of 
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adhesive joint and earned out static fracture tests on acryl and epojq^ adherends 
and epoxy resin adhesives From the experimental results, they concluded that 
total strain energy release rate was a dominant parameter, which characterise 
the fatigue crack propagation in adhesive joints Further, they inferred that the 
critena for fatigue crack growth and static fracture toughness were different 
Using beam theory, Fernlund and Spelt (1991b) presented a method to predict 
failure load of structural adhesive joints in design application However, the 
calculation of the strain energy release rate and the analysis was limited to the 
Imear relation of load and displacement Further, Fernlund and Spelt (1991c) 
conducted experiments for the study of crack propagation in cracked lap shear 
(CLS) specimen subjected to a range of loading conditions They concluded 
that the critical strain energy release rate was not a material property but 
increased with decrease in thickness of adhesive layer and with decrease in 
GilGji ratio Further, they found that loads required to propagate a crack in 
a cracked joint for a specific adhesive system, provided a conservative estimate 
of the failure load for the corresponding uncracked joint 

Miskioglu, et ol (1991) proposed a linear elastic model for a bimatenal 
interface crack to determine the direction of crack propagation during surface- 
layer removal by scraping Photoelashcity method was used to determine the 
loading at which a scraper would need to be delivered to the surface layer 
for continued mterfacial crack propagation and the results were compared with 
those obtained from finite element analysis In both the analyses fnctionless 
contact between wedge and layer was assumed The data-collection region 
was r^tricted to a zone very close to the crack-tip There was discrepancies 
between the measured and calculated fringe orders far away from the crack-tip 
Similar discrepancies were observed in the very close proximity of crack-tip, 
>i^%ere the solutioqs were influenced by material and geometrical non-hnearities 
Stress and displacement fields for an mterfacial crack opening in an anisotropic 
bimatenal were represented by two real matrices respectively by Qu and Li 
(|9il)^ Thiy fi^nd teat«tet^©*c%«ild4m.cha*aeterise(l hy ihe'tractions along 
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the interface for dislocation and crack problems However, exact number of 
non-dimensional parameters needed to describe general anisotropic bimatenal 
system was not considered 

Wassell, et ol (1991) developed a test method, which measures the crack 
growth rate for an adhesive material in a bondline that contained a cohesive 
crack They used finite element technique to verify the experimental results 
(material properties) The test was found to be sensitive to crack propagation 
rates of the order of 10'^ mm/cycle However, this work was limited to Imear 
elastic adhesive material Bishopp (1992) carried out expenments with equal 
adherends double cantilever beam (DCB) specimen, which generated mixed 
mode fracture data called as fatigue envelope He observed that the mode 
ratio and critical energy release rate (Gc) were independent of crack length 
and fracture toughness (Gic) was independent of adherend thickness He also 
observed that critical energy release rate Gc was increased with increasing 
nominal phase angle of loading and Gjjc was approximately three times higher 
than Gic Further, he found evidence of effects of crack-face fnction in mode 
II However, these results were based on brittle nature of structural epoxy 
adhesive Davies, ct al (1992) attempted to establish the reproductabihty of 
values for mode I and mode II tests on two materials m different laboratones 
Large amount of scatter was observed for the results obtained from end-loaded 
split (ELS) specimen They determined Gjc and Gjjc values from DCB and 
end notched flexure (ENF) specimen, respectively They also observed that the 
values of Gw and Guc increased slightly with increase m specimen thickness 

Dung (1992) introduced a numencal model of void growth and analysed 
the behaviour of a smgle (yhndncal or sphencal or ellipsoidal void m a plastic 
matenal containing a periodic array of voids The effects of the triaxiahty on the 
void growth was accounted with the microscopic and macroscopic factors such 
as the stress components, the hardening exponent, the volume fraction, shape 
and spacing of the voids and the constitutive softenmg This model was apphed 
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to a three-dimensional problem with the assumption of constant distnbution 
of the stress and strain fields in solid throughout the process of deformation 
He found that the predicted results were satisfactory for volume fraction and 
yield stress However, estimated void aspect ratio did not show the flattening 
behaviour of the void In the analysis of adhesive joints Krenk (1992) modelled 
the adhesive layer as a continuous distnbution of tension and shear spnngs He 
concluded that the stress concentration and the energy release rate depended 
on the adhesive thickness and stiffness This study revealed that large adhesive 
stiffness and small layer thickness resulted in a stiff bond with large localised 
stress concentrations Further, a double cantilever beam specimen under mode 
I and mode II loadings and a symmetric lap joint under mixed mode loading 
were analysed 


Edde and Verreman (1992) carried out stress analysis using beam theory 
approach for a clamped crack lap shear (CLS) adhesively bonded joint They 
used a crack-closure technique to obtain the total energy release rate in various 
modes for a crack propagating along the bond line However, the results 
showed good agreement with finite element solution only for small crack lengths 
The fatigue growth of planar defects within adhesively bonded T-peel joints 
was modefled by Ghlchrist and Smith (1993) In accordance with linear elastic 
fracture mechanics, crack growth was characterised by stress intensity factors 
and strain energy release rates It was assumed that defects within the adhesive 
were likely to be imtiated cohesively within the adhesive fillet region, propagated 
to the surface, formed a through-width defect and rapidly grew through the 
remaining bond They found that the static and fatigue strength of joints were 

maximum, when the adhesive filled the fillet region This analysis was limited 
only to cohesive type of failure 


and Zang (195.3) charactemed stresses single lap joints m the 
^t^nce of a void by Ovo non-dimensional parameters The influence of the 
locifion ankt see 'of tile voitf ofi ^efc stress Was also studied for identical and 
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non-identical adherends Whitcomb and Woo (1993) studied the debond growth 
in adhesively bonded tubular joints using geometrically non-linear finite element 
analyses They discussed behaviour of tubular and plate joints under tensile and 
flexural loads They also observed that the strain energy release rates converged 
rapidly with mesh refinement^ Richardson, et ol (1993) assessed the validity of 
modelling adhesive joints as two-dimensional problems Comparison of both the 
adhesive stress distnbution and the energy release rates associated with the crack 
growth from two- and three-dimensional analysis were made They suggested 
to cariy out a three-dimensional analysis to obtain modified load distnbution 
and apply these distributions to refined two-dimensional model, which would 
take care of non-uniform load transfer Refinement of three-dimensional mesh 
would influence the accuracy of the calculated load transfers and thus the final 
results This important aspect was not considered in their analysis 

Usmg a high-magnification video camera Chai (1993) evaluated damage 
at the tip of cracks m adhesive bonds deformmg m shear He concluded that 
ultimate shear strain was a key material property for controlling crack growth 
in brittle and ductile epoxy resin Qu and Bassani (1993) studied the aspects 
of fracture mechanics of interface cracks in anisotropic bi-matenals A general 
discussion of the oscillatory nature of the crack-tip fields, crack-face contact, 
effects of specimen size, interfacial toughness and matenal properties were 
presented They concluded that there was a coupling of all three crack-tip 
fracture modes in the normal anisotropic bimatenals, whereas in isotropic case, 
mode III was decoupled from mode I and II 


1.2.5 Parametric study 

Influence of bond thicbiesa on the stress distnbijtion in single lap adhesive 
joints ^was mvestigated by Ojalvo and Eidmoff (1978) They^ found that its^ 
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effect was more pronounced for the joints with short overlaps, thick adherends 
and stiff adhesives Sen and Jones (1980b) carried out a parametric study 
and design of a double lap joint, which was bonded with a visco-elastic 
adhesive They studied the qualitative influences of the vanous geometrical, 
matenal parameters and creep strains on the possible failure modes of the joint 
They found that the most important parameters were (i) ratio of adherend 
modulus to equivalent property of visco-elastic adhesive and (ii) thickness of 
the adhesive Using the double cantilever beam (DCB) adhesive joint specimen 
and scanning electron microscope (SEM), Chai (1986) studied the effect of 
bond thickness on the fracture behaviour of ductile adhesive Depending on 
the thickness of the adhesive continuous crack growth, unstable crack growth 
and a combination of both were observed Qu and Bassani (1989) studied the 
two-dimensional problem of an interfacial crack between two anisotropic elastic 
solids and derived necessary and sufiicient conditions such that the crack-tip 
fields were non-oscillatory However, this study was limited to perfect continuity 
of displacements across the interface 


Using the theory of elasUcity, Sawa, et at (1989a) analysed stress 
distnbutons of adhesive butt joints of dissimilar solid (ylindeis subjected to 
eirtemal tensile load The effects of adherend stiffness, adhesive thickness, ratio 
of adherend and adhesive Poisson’s ratios, normal aspect ratio and shear stress 
distnbuuon on the joint strength were studied by numencal calculations They 
found that the stress smgulanty became large at the edges of the interface where 
the adherend stiffiiess was large and fracture occurred at the mterface MaU 
»d R^amurthy (1989) earned out experimental investigabon on composite 
bonded joints to study the effect of bond thickness on debond growth rate 
Mder cyclic loading and critical strain energy release rate under static loading 
^ey observed that fracture and fatigue strengths were dependent on bondline 
thickness Using two stress smgulanty parameters Hatton (1991) investigated 
the shear strength of single lap bonded joints He studied the influence of 
adhesive layer thickness, overlap length and bonding edge angle on joint strength 
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He found that the shear strength of single lap joints got decreased when the 
overlap length was more than the substrate (adherend) thickness Further, he 
observed that fracture occurred at the upper bonding edge, when the bonding 
edge angle was less than 7r/3 and it occurred at the lower bonding edge, when 
the bonding edge angle was greater than 7r/3 It was suggested that shear 
strength could be improved by decreasing the bonding edge angle to less than 

7r/6 


Kim, ct al (1992) employed small fatigue test specimens and investigated 
the effects of adhesive thickness and adherend surface roughness on the fatigue 
strength of the tubular single lap bonded joint These tests were performed on 
single lap joint, single lap joint with scarf, double lap joint and double lap joint 
with scarf Chai (1992) studied mixed mode fracture and interaction between 
fracture modes (opening and shearing modes) in adhesive joints with vanous 
adhesive thicknesses under combined loadmgs He found that fracture energy 
was largely controlled by plastic deformation m the adhesive matenal Munz 
and Yang (1992) examined stress smgulanty occurring at the free edge of the 
interface of dissimilar materials under mechanical and thermal loadmgs They 
inferred that size of the component had a significant effect on the stress at the 
free-edge of the interface, especially for thermal loading 

Choi and Thangjitham (1993) presented a numerical method for analysing 
the steady state anisotropic thermo-elasticity problem of fibre remforced composite 
lammate containing a crack embedded within a given matrix nch interlaminar 
region It was shown that the interlayer model provided a resulting crack-tip 
stress field with the standard order of square root smgulanty Owing to the 
unique structural charactenstics of lammate stacking sequence, couplmg effects 
could anse from matenal properties mismatch between adjacent fibre layers, 
leading to the simultaneous presence of all fracture modes This aspect was not 
considered and further his work was limited to hnear elastic fracture mechanics 
(LEFM) 
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1.2.6 Environmental effects 

Parker (1990) exposed carbon fibre reinforced epoxy composite joints bonded 
with epojqr adhesives to 50° C and 96 percent relative humidity for 3 years 
and tested them From the experiments, he concluded that strength reduction 
was dependent on exposure time, test temperature and changes m the adhesive 
properties rather than m adherends Depending on the adhesives used, he 
observed that the joints failed either in the composite adherend or m adhesive 
or both Rao and Crocker (1990) descnbed a numencal model to predict the 
natural frequencies and the mode shapes of a bonded lap jomt for free vibration 
They compared numerical and expenmental results for natural frequencies of 
lap bonded joint in beam system They found that the natural frequencies of 
the system got mcreased with an increase m the overlap length Further, the 
damping ratio due to the presence of bonded joint in the system did not have 
significant effect on the system response Using numencal shape optimisation 
techmques, Groth and Nordlund (1991) performed shape optimisation of bonded 
joints By changmg the profiling of the adherends, they obtained strong and 
light joints Types of jomts studied were the single lap, double lap and double 
strap A console bonded to a ngid wall was also examined However, this 
analysis needed a suitable initial optimal geometry as the input, else the objective 
function could find a local (not necessarily global) minimum value 

Adams, et al (1992) found significant changes in stress state of lap joints 
caused by thermal effects The stress/strain properties of the polymeric adhesive 
also got changed enormously with temperature resulting in variation of strength 
of single lap shear joints They observed that strain rate capability was high but 
the load capability was low The reverse occurred at low temperatures Cheriy 
and Ye (1992) studied the behaviour of high temperature anaerobic (without 
otygen) adhesives, which were cured at room temperature in the absence of 
oxygen Experimental results from a commercial high temperature anaerobic 
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adhesive showed that the adhesive polymer became brittle as a result of fully 
cured state, and the adhesive near the edges of the bond did not get cured 
properly because of penetration of oxygen inhibiting it from curing The size 
of this uncured zone was found to be a key factor affecting the strength of 
the joints at room temperature However, this work completely neglected the 
influence of interfaces and the adherend properties Conley, et al (1992b) 
measured subcritical crack growth behaviour along polymer/glass interfaces for 
various polymer adhesives at different relative humidities They used the power 
law relationship between crack velocity and the strain energy release rate to 
discuss the delamination of thin films under residual tensile stress The shell 
epoxy was found to be most resistant to delamination as compared to other 
epoxies Probable reason for this retarded crack growth was mentioned to be 
high fracture resistance and voids at the interface 

Finger like crack growth behaviour in a urethane acrylate/glass interface 
was studied by Conley, et al (1992a) as a function of humidity and applied strain 
energy release rate C?, for the study A four point flexure specimen consisted 
of two glass plates bonded together with a finger like crack introduced at the 
adhesive/glass mterface They observed that finger spacing was independent 
of crack velocity Though crack growth at the interface was characterised by 
finger like perturbations along the crack front due to the viscous nature of 
this adhesive, calculation of G did not include viscous matenal behaviour of 
the adhesive Ritter and Conley (1992) earned out study of moisture assisted 
crack propagation at polymer/glass mterfaces They concluded that, above a 
threshold value of applied G, moisture assisted crack growth at the interface was 
dependent on the applied G by a power law expression, where the proportionality 
constant depended on adhesive tested and the relative humidity Polymer 
adhesive tested were epojty acrylate, urethane aciylate, Devcon epoxy and Shell 
epoxy They observed that high humidity enhanced the crack growth rate Su, 

ct al (1992) studied the effects of aging and environment on the fatigue life 

1 

of adhesively bonded joints They observed that those adhesives cured with 
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polyamine hardener showed unproved fatigue life and had high initial strength 
and Young’s modulus Further, the adhesive that had high moisture absorption, 
low initial Young’s modulus and low heat deflection temperature, lost it’s fatigue 
resistance completely Wingfield (1993) discussed the problems associated with 
the surface pretreatment of fibre reinforced composites of adhesive bonding 
To obtain a strong joint, they suggested removal of surface contamination such 
as silicones and flurocarbons However, the abrasive surface treatments like 
alumina glitblastmg and sihcon carbide abrasion would change the mechanical 
properties of the bulk composite 


1.2.7 Quality of adhesive bonding 


Various methods and specimens, for obtaining the mechanical characteristics of 
structural adhesives needed for the design and analysis of adhesively bonded 
structural joints, were described by Jeandrau (1986) Groth (1988b) designed 
a test specimen of double contoured cantilever beam, which resulted m a 
constant value of stress intensity factor in Mode I for a crack length up to 
90 mm of the total specimen length 100 mm However, this test specimen 
was vahd for mode I loading condition only The results obtained from 
five epo^ bases structural adhesives from single lap jomts were presented 
by Jeandrau (1991) Thomson and Thomson (1991) attempted to develop 
a quantitative accept/reject methodology based on fracture mechanics, which 
considered elements of matenal characterisation, stress analysis, failure modelling 
(including time and temperature dependent features) and fractography They 
also reviewed the earlier research related to improvement of techniques to 
measure the bond strengths They concluded that the correlations between 
ultrasonic data and strengths for adhesive and cohesive types of failures persisted 
in <Jompfex adhesive structures 
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Jiao and Rose (1991) introduced an interface layer model for evaluation 
of adhesive bonding quality However, this process was sensitive to frequency 
range, angle of incidence, wave type, type of specimen and signal processing 
procedures employed during the experiment Using Acousto-ultrasonic (AU) 
non-destructive technique Tiwan, et al (1991) examined the quality of adhesively 
bonded sheet metal used for automobiles AU is a hybrid of acoustic emission 
and ultrasonic techniques This technique assured quality of bonded joint and 
requires access to only one side of the component They observed that bond 
quality was related to the attenuation of the energy content of the received 
signal for mixed mode failure However, surface roughness and cleanliness 
did affect the AU signals and thus influenced the results Hanneman and 
Kinra (1992a) proposed a technique for ultrasonic non-destructive evaluation 
of adhesively bonded joints, which was applicable to perfectly elastic and linear 
visco-elastic materials, and also for isotropic and anisotropic materials 


Using piezoelectric transducers, transmissions of a plane longitudinal 
ultrasomc wave, through lap joint immersed m water, were measured by 
Hanneman and Kinra (1992b) to deduce the jomt properties such as adherend 
and adhesive thicknesses However, these results were sensitive to frequency 
range used, resonance mode, adhesive material properties and attenuation 
Parker (1993) studied the environmental durability of aluminium joints with 
different pretreatments He compared the durabilities of bonded joints made 
with aluminium pretreated by ceramic pickle and anodizing pretreatment 
Considering time of failure under stress, residual strength and mode of failure, 
he concluded that anodizing pretreatment resulted in a longer life of the bonded 
joint under normal laboratory environment 
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1.2 8 Design 


A method was developed by Fmdlater (1987), which enabled the effect of joint 
dimension to be compared with the load carrying capacity of a particular joint 
configuration However, these did not yield information about the levels of 
stress in a joint nor they were substitutes for numerical tool such as finite 
element analysis Bigwood and Crocombe (1989b) suggested a two parameter 
design formula which calculated the adhesive shear and tensile stresses in the 
overlap region However, the errors were introduced in case of adherends 
mismatch T his design formula yielded accurate stress values and was valid for 
balanced joint configurations Hashim, et al (1990) studied the design and 
assessment methodologies for adhesively bonded structural connections They 
concluded that the use of hot cured epoxy adhesives was both technically and 
practically viable for replacing fillet welded connections in lightly loaded steel 
structures Depending on the specific application and requirement a limited 
experimental and numerical technique used with a reliability assessment method 
might be adequate and might not require a higher factor of safety However, 
this work did not provide enough case studies to support the new idea of 
rehabihty assessment m joint analysis 


Adams (1992) discussed the importance of finite element analysis, reliable 
data (adhesive properties) for the joint analysis and the need to compare the 
results of finite element analysis with those obtained from experiments Lee 
and Selby (1992) devised a logical technique, which could be used to select 
most suitable adhesive for a purpose This was an implementation of expert 
^stem Therefore, this technique had the major hmitations such as need of an 
ei»:ellent data bank and regular updating of data 
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1.3 Objectives of the present work 

13 1 Status of the literature 

From a critical review of the available literature following observations regarding 
the status of understanding of the response of adhesive joints could be made 


• Vanous parameters like joint configuration (single lap, double lap, stepped 
lap, butt, scarf etc ), adhesive thickness, overlap length, quality of surface 
preparation, bonding quality, material properties of the adhesive and 
adherend, and shape of adherends play an important role on the performance 
of the adhesively bonded joints 

• Polymeric adhesives exhibit non-linear behaviour under high stress and/or 
extreme thermal conditions 

• There is significant variation m the results based on elastic, elasto-plastic 
and visco-plastic analyses Results of non-hnear responses of these joints 
are spread over a wide range 

• Fimte element technique appears to be one of the most versatile mathematical 
tool to handle the non-hnear behaviour and the parametric study of the 
adhesively bonded joints 

• Most of the analyses reported in literature are either completely theoretical 
or computational Very few of them have been corroborated by the 
expenmental studies 

13 2 Obiect ives of the present work 

It is^felt that, there is a need to systematically develop a finite element 
technique to study the behaviour of adhesively bonded joints With this view, 
in this thesis an attempt is made to 
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• develop an automatic mesh generation program, which could locally refine 
the finite element mesh near the critical zones such as sharp edges, crack- 
tips, interfaces and discontinuities, 

• be able to take into account various combination of adherend matenals 
like metal-metal, coraposite-metal and composite-composite, 

• incorporate fracture mechanics approach with possible failure criteria to 
predict the failure of the joint, 

• mclude non-linear behaviour (elasto-plastic and visco-plastic) of the polymeric 
adhesive, 

• study the effects of geometncal and matenal parameters on the performance 
and strength of these joints, 

• examine the growth of plastic zones under different displacement boundary 
conditions and as the crack propagates, 

• conduct some simple expenments to determine the strain energy release 
rates in pure modes I and II, respectively, 

• determine the bond strength of the joint, and 

• arnve at a simple theoretical fracture criterion, which fairly correlates well 
with the observed experimental results 


1.4 Layout of the thesis 


ftesenc work includes (a) elastic analysis (b) elasto-plastic analysis (c) visco- 
plastic analysis and (d) expenmental study of adhesively bonded joints 
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Chapter 1 serves to introduce the adhesively bonded joints Section 1 1 
presents the motivation to carry out the present research work This Section 
provides bnef introduction of adhesively bonded joints, composite material, 
non-hnear finite element analysis, fracture mechanics approach and experiments 
to determine fracture parameters and the joint strength Section 1 2 contains 
the review of the literature In this Section literature related to (i) linear elastic 
analysis, (ii) non-Iinear (material) analysis (iii) non-hnear (geometrical) analysis 

(iv) fracture mechanics approach (interface failure, fracture and delamination), 

(v) parametnc study, (vi) environmental effects, (vii) quality of adhesive bonding 
and (viii) design of adhesively bonded joints (single lap, double lap, butt, scarf 
and tubular joints) are discussed Matenal non-lineanty includes both elasto- 
plastic and visco-plastic behaviour Section 1 3 serves to focus on the objective 
of the present research work 

Chapter 2 deals with the linear elastic finite element analysis of the joint 
and the fracture mechanics approach Section 2 1 provides a brief introduction 
of the linear elastic analysis of the bonded joints Section 2 2 presents the 
crack closure technique and finite element formulation Procedure followed to 
compute the failure loads and the automatic finite element mesh generation for 
bonded joints and the local refinement of mesh at critical zones of the joints 
are discussed Section 2 3 deals with the validation of computer code generated 
for the (linear elastic) finite element analysis Stress and strain distributions in 
the adhesive across the overlap length of single lap and stepped lap joints, are 
investigated m subsections 231 and 23 2, respectively Also, these distributions 
are investigated across the thickness of the adhesive 

Section 24 intends to present the modelling of crack growth (quasi- 
static crack propagation) and validation of the finite element code computing 
strain energy release rates Crack is modelled with paired nodes Employing 
Irwm’s crack closure technique strain energy release rates are computed In 
subsection 2 41a single edge notch (SEN) specimen is considered tc^ deteraime 
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the minimum refinement required to obtain accurate (close to the close form 
solution) strain energy release rate Moreover, in subsection 2 4 2, a cracked 
lap shear (CLS) specimen is considered to validate the finite element code 
calculating strain energy release rate Both plane stress and plane strain models 
of CLS specimen are investigated Effect of mesh refinement on the computation 
of strain energy release rates is investigated in Section 2 5 Parametric study 
IS carried out in Section 2 6 These include effects of (i) crack location, 
(u) thickness and Young’s modulus of the adherend, (iii) stacking sequence 
in composite adherends, (iv) Young’s modulus of the adhesive, (v) overlap 
length of the adhesively bonded joint, and (vi) thickness of the adhesive on the 
strength of the adhesively bonded joints These parameters are considered one 
at a time Conclusions of linear elastic analysis of these joints are discussed in 
Section 2 7 

Elasto-plastic analysis is presented in Chapter 3 Section 3 1 provides 
a bnef mtroduction of the non-linear elasto-plastic analysis Finite element 
model, flow chart and the algorithm for elasto-plastic analysis are illustrated 
in Section 3 2 In this Section stress-strain relation, convergence of the elasto- 
plastic analysis and effect of strain hardening on bond strength are also discussed 
In Section 3 3, finite element code is validated with a beam problem from 
literature Section 3 4 shows the spread of plastic zone based on elasto-plastic 
analysis Growth of plastic zone for three different boundary conditions and 
for crack propagation are illustrated in subsections 3 41 and 3 4 2, respectively 
In Section 3 5 the effects of vanous parameters (considered one at a time) 
on the strength of composite to composite joints are discusses Subsections 
351, 35 2, 353 and 3 5 4 are dealt with influence of (i) stacking sequence in 
composite adherend, (ii) fracture toughness of adhesive (iii) overlap length and 
(Iv) adhesive thickness on joint strength, respectively Concluaons related to 
the elasto-plastic analysis are discussed in Section 3 6 

Visoo-plastic ariailysis is presented m Chapter 4 Section 4 1 provides a 
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brief introduction of the non-hnear visco-plastic analysis Strain rate, visco- 
elastic, visco-plastic, visco-elasto-plastic models of adhesive material, flow chart 
and the algorithm for visco-plastic analysis are illustrated in Section 4 2 In 
subsection 4 31, finite element code is validated with a perforated plate problem 
from literature Further, the visco-plastic stress analysis results are compared 
with those available from literature for single lap and stepped lap joints in 
subsections 43 2 and 4 3 3, respectively Section 4 4 discusses the effects of 
various parameters (considered one at a time) on the visco-plastic strength of 
composite to composite joints Subsections 441, 44 2, 443 and 444 are 
dealt with effects of (i) stacking sequence, (ii) fracture toughness of adhesive 
(ill) overlap length and (iv) adhesive thickness on joint strength, respectively 
Conclusions made from the visco-plastic analysis are discussed m Section 4 5 

Experimental work is presented in Chapter 5 A brief introduction of 
expenmental study is discussed in Section 5 1 Sections 5 2 and 5 3 describe 
the theories mvolved in determination of cntical strain energy release rates 
m pure modes I and II, respectively Experimental determination of critical 
strain energy release rates m modes I and II are also reported In Section 5 4 
determination of strength of the double lap carbon fibre reinforced plastic 
(CFRP) joints These results (bond strength) are compared with those obtamed 
from the Imear elastic finite element analysis considenng various fracture cntena 
in Section 5 5 

Concluding remarks are presented in Chapter 6 Section 6 1 contains the 
results and discussion of the present work Also comparison of elastic, elasto- 
plastic and elasto-visco-plastic results and the computation (execution) times 
are earned out and listed in this Section Section 6 2 presents the hmitations 
of the present analyses Finally, in Section 6 3 few suggestions are listed for 
future extension of this work. 
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(a) (b) (c) (d) 


Figure 1 1 Joints of equivalent load capacity with various fastenings, (a) overlap 
5 in, resorcinol-formaldehyde adhesive joint (b) overlap 9 in, four | m bolts (c) 
overlap 11 m, two 2^ in spht ring connectors and one | in bolt and (d) overlap 
12 ^ in, twenty-four 4 in 7 gauge nails dnven in double shear into drilled holes 


(a) 

(b) 

(e) 

(d) 

(e) 

(f) 



Figure 1 2 Schematics of some common adhesively bonded joints (a) single 

lap (b) double lap (c) butt with straps (d) stepped (e) scarf and (f) tubular lap 
joints 






Chapter 2 


Elastic analysis 


2.1 Introduction 


From a critical review of literature, it is observed that the finite element 
technique is found to be a suitable tool for the analysis of complex structures 
such as bonded joints In addition, fracture mechanics approach provides an 
accurate estimation of strength of these joints However, results, reported by 
various researchers m the literature, are found to be distributed over a wide 
range (Johnson, 1987 and Davies, et ol , 1992) Therefore, it is important to 
estimate the strength of these joints accurately 

Though, It is useful to include the material non-lineanty of the adhesive, 
as a first step, in this Chapter, both adhesive and adherend are considered 
to be linear elastic in behaviour In addition, elastic analysis of the joints 
confirms the importance of vanous parameters and also forms the basis for 
non-hnear analysis in the later Chapters Finite element code with fracture 
mechanics approach is developed and validated with existing results from the 
literature Elastic analyses of a smgle lap, a stepped lap and a double lap joints 
are carried out Various metals and fibre reinforced plastic (FF-P) composites 
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are employed as adherends in the analyses In addition, seven fracture cnteria 
are also considered Effect of mesh refinement on the computed strain energy 
release rates is investigated Influence of various parameters, such as (i) crack 
location, (ii) thickness of the adherend, (iii) Young’s modulus of the adherend, 
(iv) stacking sequence in composite adherends, (v) Young’s modulus of the 
adhesive, (vi) overlap length of the adhesively bonded joint, and (vii) thickness 
of the adhesive on the strength of the adhesively bonded joints, are also 
investigated 


2.2 Formulation 


I Crack closure technique 


Wang and Crossman (1982) discussed Griffith’s strain energy release rate 
concept, Gnffith’s cntenon, Irwin’s crack closure integral and influence of shape 
of the body, loading, location, size and orientation of crack on the stress intensity 
factors They also emphasised the crack initiation and crack growth criteria 

Gnfl!ith’s criterion is employed to predict the crack growth (Broek, 1984) 
The stability of the crack growth is governed by the following conditions 


G(a + Aa) < Gc 
G(a + Aa) > Gc 
G(u + Aa) ’ = Gq 


stable crack growth 
unstable crack growth 
neutral crack grwvth 


( 21 ) 
(2 2 ) 
( 23 ) 
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where, Gc is the critical strain energy release rate and Aa represents 
an infinitesimal crack extension from the initial crack size ’a’ as shown in 
Figure 2 lb 

For isotropic raatenals the distributions of stresses near the crack-tip could 
be expressed in terms of stress intensity factors Kj, Kjj and Km (Broek, 1984) 
Strain energy release rates m mode I loading for plane stress and plane strain 
cases are expressed in terms of the crack-tip stress field and functional of the 
SIF’s on the crack size 


ri„-\ - ^ 

(plane stress) 

(24) 

Gi(a) = 

(plane stram) 

(2 5) 


a, E and p represent the crack length, Young’s modulus and Poisson’s 
ratio, respectively For a general case of mixed mode fracture the total stram 
energy release rate is given by 


G = Gj + Gji + Gju 


(2 6 ) 


where G/, Gn and Gm denote the components of mode I, II and III 
fracture, respectively Ii|^eneral in fibre remforced composite materials, fracture 
IS a combination of all the three modes Irwin showed that the elastic strain 
energy released dunng an incremental crack extension is equal to the work 
done m closing the incremental crack (Wang and Crossman, 1982) This crack 
closure method provides an approadi for the direct determmation of strmn 
energy release rates from only crack-tp stresses and displacements^ When the 
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crstck IS opened as depicted in Figure 2 lb, the work done to close the crack 
opening could be written as 


1 “ 

AW = ^J^At^da (2 7) 

0 

aO' is the relative displacements between the mating crack surfaces along 
Aa ^ IS the surface stress distribution along infinitesimal crack extension Aa 
when the crack is closed as shown in Figure 2 Ic According to Irwin’s crack 
closure equivalence, the available energy release rate G for a crack size a is 
expressed by 


G{a) 


lim ^r-r- 

Aa —>0 2,^(1 


Aa 

j9A0da 


(2 8 ) 


Substituting the components of the surface stresses ^ and the corresponding 
relative displacements AU into Equation (2 8) yields components of G The 
crack closure representation is convenient for adaptation in numerical computation 
One among other methods to evaluate G is J integral technique, which is 
discussed by Johnson (1987) 


I Finite element technique 

A numencal technique to calculate the strain energy release rate involves a finite 
element solution of Irwin’s crack closure integral as given by Equation (2 8). In 
the finite idement represpntation the ctmtlmiouis stress and displacement fields 
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of the solid are approximated by the nodal forces and nodal displacements, 
respectively Figure 2 2a represents the finite element modelling near a crack-tip 
region A crack of length a is shown with the crack tip at node c The finite 
element solution determines the displacement components (u, v, w) of the crack 
tip node c An incremental crack extension Aa is introduced replacing the crack 
tip node c with two separate nodes / and g as depicted in Figure 2 2b For 
the new crack geometry the finite element solution for the nodal displacements 
(u,v,w)f and (u,v,w)g are found for nodes / and g respectively The crack 
opening is then collapsed by applying equal and opposite forces at nodes / and 
g such that their common displacements match the displacements found earlier 
at c These forces correspond to the internal nodal forces which exist at node 
c before it opens The work required to close the crack opening (Figure 2^Lc) 
is approximated by 


AW ~ ^[Fx(uf - Ug) + Fy(Vf ~ Vg) + F,(wf - Wg)] (2 9) 

where Fx, Fy and Fg are normal components of nodal forces required to 
bring nodes / and g together The energy release rates for the three crack 
extension modes are approximated by 


Gj = 

2Aa 

( 210 ) 

Gij == 


( 211 ) 

Gin = 

»,)) 

( 2 12 ) 


This method does not require the calculation of the stresses and strains 
as stress and strain fields are approximated by nodal forces and displacements. 
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respectively Therefore, this method has the advantage of less computation and 
more accuracy as compared to other methods, where one needs to compute 
the stresses and strains 


I Prediction of failure load 

A computer code is developed to predict the failure load, which is detailed 
m the following main algonthm and the flow chart is depicted in Figure 2 3 
Though the algonthm and the flow chart are meant for linear elastic material, 
these are extended later to mclude elasto-plastic and visco-pIastic matenals 
by adding corresponding modules Details of these additional modules are 
explained in Chapters 3 and 4 The loads are updated in each iteration, strain 
energy release rate is computed for a pre-asigned crack path along the interface 
of the adhesive and adherend The computed strain energy release rate value is 
compared with its cntical value within 1 percent tolerance Once this condition 
IS satisfied paired node is opened (Wang, 1982) and the corresponding failure 
load and crack length are noted 

The failure load is observed to be strongly influenced by the fracture 
criten^^ and the critical strain energy release rates Seven possible fracture 
cnteria^are considered and the corresponding failure load vs crack length 
variations are compared in Figure 2 4 As will be shown later, in Section 5 4, 
the fracture criterion Gj + Gn > Gjc + Gnc is the most suitable one, among 

the seven critena considered Therefore, m the present analysis this failure 
criterion is employed 


I Mam algonthm 

i 

'Step 01 read input dafa such as geometrical, material and mesh 
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parameters and various options 

step 02 calculate stiffness matrix {K] 

step 03 update applied load AP 

step 04 solve for displacement 

A6 = [Ky^AP 

step 05 update displacement, strain and stress 
8 = 8 + AS 
e = e + Ae 
a = a + Aa 

step 06 if plastic flow occurred 

(i e a major portion of adhesive has become plastc) 
or, maximum allowed load is achieved 
or, all prescnbed paired nodes are opened 

(i e critical crack length is reached) then go to step 09 

step 07 * compute the strain energy release rate (G) 
and compare with its cntical value (Gc) 

if 1 OlGc > G > Gc then open the paired node and go to step 
03 

else go to step 03 

step 08 if elastic analysis else go to step 07 
if elastO'plastic analysis go to 
step 01 of elasto-plastic olgonthm m Section 3 2 
if visco-plastic analysis go to 
step 01 of visco-plasttc algorithm m Section 4 2 

step 09 • stop 
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I Automatic mesh generation 

A finite element mesh generation computer code is developed, which automatica- 
lly generates refined meshes for problems considered in the present work 
vu (i) a cantilever beam, (ii) a perforated plate, (iii) a cracked lap shear 
(CLS) specimen, (iv) single edge notch (SEN) specimen, (v) adhesively bonded 
single lap, double lap and stepped lap joints This makes use of Chebyshev’s 
function to generate necessary mesh refinements at crack-tips and other stress 
concentration zones Thus, with less number of elements one could achieve 
reasonably accurate results A typical finite element mesh of adhesively bonded 
joint with more refinement at edges is shown in Figure 2 5 of next Section 
The code also generates data for the paired nodes at the interfaces and nodal 
forces and displacement boundary conditions The paired nodes, are opened 
in a secjuence along a pre-defined (crack) path, simulating the crack growth, 
which facilitates calculation of the strain energy release rates 


2.3 Validation of finite element code 


To validate the computer code of finite element developed a single lap joint 
and a stepped lap joint are considered Stress and strain distributions along 
the overlap lengths are compared with those available m open literature 


23.1 Single lap joint 

To check the accuracy of computer code of finite element analysis, a single 
lap joint (Figjure 2 5) is employed Details of the joint are listed in Table 2 2 
In Figure 2 6 the computed peel stress and shear stress distributions along the 
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overlap are compared with those of Hiregoudar (1993) Results of the present 
analysis (with 1816 CST elements) are in good agreement with Hiregoudar 
(1993) The differences are of the order of 6 and 12 percent, for the peel and 
shear stresses along the overlap length, respectively (Figure 2 6) Maximum 
difference is observed near the end region of the overlap This difference is 
found to be minimum near the central portion of the overlap Because of the 
effective refinement in the finite element mesh near the edges and near the 
interfaces of adhesive and adherends, m the present analysis, predicted stress 
distributions are expected to be more accurate 


2.3 2 Stepped lap joint 

Further, vahdation of the computer code of finite element developed is earned 
out on a stepped lap jomt Figure 2 7 illustrates the joint with a typical finite 
element mesh In this joint, thick adherends are employed and thus it is also 
called as thick adherend shear test (TAST) specimen Specifications of the joint 
are listed in Tkble 2 3 In Figures 2 8 and 2 9 the stress and strain distributions 
along the overlap at the interface of adherend and adhesive are compared with 
those of Su and Mackie (1993), respectively Similarly, in Figures 2 10 and 2 11 
stress distnbutions along the overlap at the midline of adhesive are compared 
Results of the present analysis (with 1816 CST elements) are observed to be 
in good agreement with Su and Mackie (1993) Near both ends (edges) of 

the overlap there are differences m these distributions Because of the 

effective refinement in the finite element mesh near the edges, th© distributions 
i^are expected to be more accurate The stress and strain distnbutions along the 
overlap and across the thickness of the adhesive are shown m Figures 2 12, 2 13 
and 2 14, respectively Because of the antisymmetric joint configuration of the 
joint about the midhne o^ the adhesive, stress and stram distributions m layers 
1 and 5 are found to be antisymmetnc Similar distributions are also observed 
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in layers 2 and 4 


2.4 Validation of crack closure technique 


2 4.1 Single edge notch specimen 

A single edge notch (SEN) specimen (Figure 2 15) is considered as a sample 
problem to validate the proposed technique The physical dimensions of the 
specimen and matenal properties are listed m Table 2 4 and a schematic of 
finite element mesh used, is illustrated in Figure 2 15 

Using modified crack closure integral and constant strain finite element 
stress analysis (AXISOL) Rybicki and Kanninen (1977) evaluated stress intensity 
factors and applied to three configurations (i) a double cantilever beam test 
specimen, (u) a finite width plate containing an embedded crack and (iii) a 
bolt fastened double lap joint containing two radial cracks They observed 
that the results computed employing higher order elements are less accurate 
as compared to those obtained with constant strain triangular (CST) elements, 
which uses nodal forces and nodal displacements Therefore, a two-dimensional 
linear small displacement plane strain CST element (Zienkiewicz and Thylor, 
1989) IS employed for the present work 

Finite element mesh with fine grids near the crack tip and paired nodes on 
the line MN (Figure 2 15) are generated A uniform far field stress o- is applied 
transverse to the crack to allow it to fail in mode I fracture Matenal of the 

specimen is considered to be aluminium Stress intensity factor is expressed as 
(Broek, 1984) 
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( 211 ) 


where 


/ 



- 0 « + IS ’ - 38 48 53 85 (^)‘ (2 1 


where a and W represent the crack length and width of the specimen, 
respectively 

The exact solution of Gj is calculated with Equations (2 4) and (2 14) Gi 
IS also computed using finite element analysis employing meshes of 520, 1360 
and 2688 elements, respectively and the results are compared The error in 
the estimation of Gj is reduced rapidly as the grid is more and more refined 
The estimated error is observed to be less than 1 6 percent for a finite element 
mesh of 2688 elements (Figure 2 15) 


2 4 2 Cracked lap shear specimen 

A four-point bending cracked lap shear (CLS) specimen (Figure 2 16) is also 
considered for the analysis Adherends and adhesive are considered to be 
aluminium and epo^ resm, respectively and their material properties and 
physical dimensions are listed in Table 25 In the present analysis, a fine 
mesh, of 5786 CST elements, is generated with proper refinement near the 
interface lines and paired nodes along the bottom interface line Schematic 
representation of finite element mesh used is illustrated in Figure 2 16 This 
problem, is investigated treating it as plane stress and plane strain models and 
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the results (strain energy release rate) are compared with the corresponding 
experimental results of Fernlund and Spelt (1991c) 

Ishikawa, et ol (1991) investigated the mixed mode fracture criteria for 
adhesive joints and carried out fracture tests on acryl adherends and epoxy resm 
adhesive They concluded that total strain energy release rate is a dominant 
parameter, which characterise the crack propagation of the adhesive joints 
In the present work also, results based on the mixed mode fracture criterion 
(Gj + Gn > Gic + Gjjc) shows better agreement with the results obtained from 
the experiment This is listed in Table 5 4, in Chapter 5 

Therefore, the total energy release rate G(= Gj + Gjj) is considered as 
the cnterion for crack propagation Variations of Gi, Gu and G with crack 
length a, are plotted in Figure 2 17, for plane strain and plane stress cases and 
compared with the expenmental result of Fernlund and Spelt (1991c) The 
difference is found to be less than 2 7 percent for plane strain model, while 
this is of the order of 5 5 percent for the plane stress model Therefore, the 
plane strain model represents the situation relatively better than the plane stress 
model It is also verified that the ratio of Gjj and Gj remains constant for 
the same thickness of adherends as mentioned by Fernlund and Spelt (1991c) 
However, in the Figure 2 17, initially strain energy release rate increases with 
crack length (till a = 2 mm), which is mainly due to the edge effect For the 
fmlure load P; (6 33 N per mm width), G is observed to be very close to the 
cri(ical strain energy release rate (G^ = 0 366 kJm"^) 

The finite element model predicts a value of G^ ~ 0 35kJmr\ while the 
results predicted from Figure 2 17 is 0 366* The present results are likely 
to be more accurate as elasticity approach is employed m the present analysis, 
where as solid mechanics approach was used by Fernlund and Spelt (1991c) 
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2.5 Mesh refinement 


To investigate the influence of the effect of mesh refinement on the computation 
of strain energy release rate, a double lap joint (Figure 2 18) is analysed with 
three different levels of refinements consisting of 680, 1184 and 1816 CSX 
elements Because of symmetry only one half of the joint is considered for 
the analysis The problem is treated as one of plane strain model The joint 
specifications are listed on Thble 2 6 The crack is introduced and allowed 
to grow and the corresponding strain energy release rates are computed for 
the above refined meshes Variations of strain energy release rates with crack 
length are shown in Figure 2 19 

It IS observed that the nature of the curves do not change wth mesh 
refinement Refinement level does not have any significant change in the strain 
energy release rates values as well In the investigation of adhesively bonded 
joint, refined mesh with 1816 elements are considered However, more refined 
mesh IS useful to analyse the free-edge effect particularly, when the thickness 
of the adhesive is very small A case study is presented in the subsection 2 66, 
where effect of adhesive thickness on joint strength is discussed 


2.6 Parametric study 




Effects of crack location, thickness of adherend. Young’s modulus of adherend, 
stacking sequence of composite adherend, Young’s modulus of adhesive, overlap 
length and thickness of adhesive on the strength of the double lap joint 
(Figure 2 18) are investigated Because of symmetry only one half of the joint is 
considered for the analysis The structure is treated as a plane strain model A 
finite element mesh consists of 1816 CST elements (Figure 2 5) is employed for 
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the an alysis From the computation of strain energy release rate, it is observed 
that there exist a crack length at which Gi and Gu values are equal This 
IS indicated in Figure 2 19 Further, before crack length attains a value of a,, 
mode I IS found to be dominating and later mode II dominates 


For convenience, three non-dimensional parameters are defined viz 


e = 

1 = 
t = 


Ei 

Er 

Lo 

Lt 

h 

tr 


(2 15) 
(2 16) 
(2 17) 


where, E„ and JSy are the elastic moduli of adhesive and adherend Lo and Lt 
denote the overlap length and the total length of the bonded joint In addition 
tg and tr represent the thicknesses of adhesive and adherend, respectively 
Assuming constant adherend thickness overlap length Lo and Young's modulus 
of adherend Er, parametnc study is carried out 


The value of e is 1/20 for Epon VIII adhesive and aluminium adherend For 
convenience, this value is also mentioned to be 1/20 even for Epon VIII adhesive 
and various CFRP composite adherends For example, if e (= EjEr)is 1/40 
implies that, the adhesive considered has a Young’s modulus of 1/2 times that 
of Epon VIII 


2 6.1 Crack location 


Effect of crack location on the strength of the double lap joint (Figure 2 18) is 
also investigated The adherends and adhesive of the joint, are considered to 
be aluminium and Epon VIII (Table 27), respectively 
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Eight different crack locations at the interfaces (Figure 2 20c) and correspon- 
ding crack growths are considered, which are described m Table 2 8 In the 
present work the criterion for failure is based on Gj + Gu > Gq Finite 
element analysis is carried out for two different boundaiy conditions viz free 
side type and constrained side type (Figure 2 20a-b) Variation of strain energy 
release rate G, with crack length a, for the free side type boundary condition is 
plotted in Figure 2 21 For the crack location 3, it is found that the stresses are 
maximum as there is no constraint imposed on the adherend at this location 
and hence free to deform 

In the constrained side type boundary condition (Figure 2 20b) mode II is 
found to be dominating For this boundary condition, vanation of strain energy 
release rate with crack length is depicted in Figure 2 22 It is observed that 
crack growths at diagonally opposite locations 1 and 3 are equally prone to fail 
Similarly, crack growths at locations 2 and 4 show same vanation of G with a 
Further, crack growths at locations 5, 6, 7 and 8 exhibit similar variations of G 
with a 

For both the boundaiy conditions (free side and constrained side types) 
the crack propagation remains stable till the crack length reaches the critical 
value ttcr and beyond that the growth becomes unstable leading to failure The 
corresponding strain energy release rate is denoted as G^r (Figure 2 21) This 
figure also indicates the ao- and Gcr for this boundary condition For eight 
different crack locations considered for both the boundaiy conditions critical 
crack length acr’s and corresponding Gcr’s are listed m Table 2 9 

From Figures 2 21 and 2 22, it may further be observed that Gj has 
increasing and decreasing nature initially for the cases of crack locations 1, 2, 
3 and 4 and smooth increasing vanation for crack locations 5, 6, 7 and 8 This 
complex behaviour at locations 1, 2, 3 and 4 can be attnbuted to free-edge 
effect The strain energy release rate Gcr, for the crack growth at location 3, 
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for the free side type boundary condition, is computed to be 195 kJ/m This 
value IS greater than that for crack growth at location 3 for the constrained side 
boundary condition The free side boundary condition results in a relatively 
weaker joint Thus, for the parametric study the crack location 3 and the free 
side type boundary condition are considered Further, the finite element mesh 
IS refined near crack location 3 instead of refinement near all edge points viz 
1, 2, 3 and 4 


2.6 2 Thickness and Young’s modulus of the adherend 

Failure loads are computed for various crack locations in the double lap joint 
(Figure 2 20c and Table 2 8) and vanous combinations of adherends of different 
thicknesses and matenals (Tkble 212) Four different thickness ratios are 
considered vu tr(upper)/tr(lower) = 0 5,1,2 and 4 Thickness of the lower 
adherend was maintained constant Metal to metal, composite to metal and 
composite to composite joints are considered for the analysis In metals, 
aluminium, copper, brass and steel are considered In the case of composite 
matenals glass-, kevlar-, carbon- and graphite-, fibre reinforced plastics are 
employed The elastic properties of the metal and composite adherends are 
listed in Tables 2 10 and 211, respectively Failure loads for metal to metal, 
composite to metal and composite to composite joints are listed m Tkbles 2 13- 
2 14, 2 15 and 2 16, respectively All the minimum (critical) failure loads among 
each set of the four locations considered are indicated with 't* It is found 
that m most of the cases location ’3’ is observed to be critical In the case 
of aluminium to brass and carbon fibre reinforced plastic (CFRP) to graphite 
fibre reinforced plastic (GFRP) joints (Tables 2 13 and 2 16) for the thickness 
ratio 05 crack locations 1’ and ’4’ yield same failure load This result could 
be refined using reduced load steps m the analysis However, the difference 
may not be insignificant 
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For the thickness ratio 0 5 in aluminium to aluminium joint crack location 
’r IS found to be critical (Table 2 13) Location ’4’ yields minimum failure 
load for the thickness ratios 0 5 and 1 of aluminium to copper joint While, 
location ’4’ gives minimum failure load for the thickness ratios 0 5, 1 and 2 of 
aluminium to steel joint Further, location ’4’ yields minimum failure load for 
the thickness ratios 0 5, and, 0 5 and 1 of carbon composite to brass and carbon 
composite to steel, respectively (Table 2 15) Failure load is observed to be 
higher for thinner adherends For example the failure loads are in increasing 
order for thickness ratios, 0 5, 1,2 and 4 


2 6.3 Stacking sequence 

Effect of stacking (layup) sequence of the composite adherends on the strength 
of the double lap joint (Figure 2 23) is investigated The adherends and adhesive 
of the joint, are considered to be carbon fibre reinforced plastics (CFRP) and 
Epon VIII (Table 2 7), respectively Unidirectional 0° and 90° layups and all 
symmetnc combinations of 0°, 90°, 45°, and -45° layups are considered for 
the study The ordering of lamina is from top layer to bottom layer The 
plies corresponding to 0°, 90°, 45°, and -45° fibre orientations are identified 
as phes, 1, 2, 3 and 4, respectively, e g [1324]s implies [0°/45°/90°/ - 45°]s' 
laminate Further, the transformed stiffriess coefficients required for the analysis 
IS obtained from the relations given in Jones (1975) 

Variation of failure load with crack length for all the 24 different combinations 
of stacking sequences in composite adherend [1324]5 are shown in Figures 2 24 
and 2 25 Symmetric laminates of combmations 1324, 1342, 1234, 1243, 1423, 
and 1432 are found to have approximately equal and higher strength than 
rest of the combinations This study indicates that if one employs 0° plies 
in the immediate neighbourhood of adhesive one can achieve a stronger joint 
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(Figure 2 25 ) Out of total 24 possible combinations the maximum bond strength 
IS observed for the laminate with [ 1324 ]s as ply sequence This is termed as 
optimal layup in the present work The joint with unidirectional laminates of 
0° as the adherends is found to be strongest 

The vanations of failure load with crack length are discussed for the 
different crack locations (Figure 2 20 c) in [Olg, OI^J, [9645, 9O49], {O45, OOl^j, 
[ 90 ^, OI5] laminate pairs constituting the top and bottom adherends. eg 
[0^5, 90^5] represents [0°]8 and [ 90°]8 unidirectional layups constitute bottom 
and top adherends, respectively Failure load vs crack length curves are 
illustrated for [Ol^, OI5], [9OI5, 90 ^], [01^, 9OI5], [ 90 | 5 , 0|s] laminate pairs for 
crack locations at 1,3 and 4 in Figures 2 26 , 2 27 and 2 28 , respectively It is 
observed that the trend of load with crack length in [O45, OI5] and {904^, 9O45] 
layups are opposite to those m [O45, 90 ^] and [ 90 | 5 , O45] layups The failure 
loads for these combinations with for the most likely crack growth and the 
corresponding crack lengths are listed in Table 2 17 and shown in Figure 2 29 


2.6 4 Ybung’s modulus of adhesive 

Effect of adhesive elastic properties on the strength of the double lap joint 
(Figure 218 ) is also investigated Strain energy release rate is plotted in 
Figure 2 30 as a function of crack length for various adhesive properties of 
double lap joint with aluminium adherend and Epon VIII adhesive, (Table 2 7) 
The values of Oa,’s and corresponding Gx’s for different e (— JSg/jE'r) values are 
listed in Tkble 2 18 Failure loads are calculated in terms of fracture toughness 
Gc of the adhesive material, for various values of e and are listed in Table 2 18 
The dimension of Gc is kJ/rn^ As the adhesive Young’s modulus £?, increases e 
also decreases This results in, decrease in G (Figure 2 30 ) and as a consequence 
corresponding failure load P; increases For various e values, failure loads are 
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calculated in terms of fracture toughness Gc of the adhesive material and listed 
in Table 2 18 Critical crack length Oo-, is observed to remain constant for all 
values of e As the Young’s modulus of the adhesive increases, e increases and 
slope of the curves G vs a decrease This indicates unstable crack growth 
and higher failure load are associated with higher adhesive Young’s modulus 
Therefore, it is advantageous to use adhesive materials with higher Young’s 
modulus values which leads to high Young’s moduli ratio e resulting in a joint 
of high strength 

Failure loads and the crack length vanations are studied for the joint of 
carbon fibre reinforced plastics (CFRP) adherends and Epon VIII adhesive 
(Thble 2 7 ), respectively For vanous e values and layups the predicted failure 
loads P/’s are shown in Tkble 2 19 It is observed that as e increases P/ 
increases appreciably for [O45, O45] laminates, while it is reversed for [OI5, 9O45] 
and [9Q45, Olg] laminates and remains fairly constant for [9O45, 9045] laminates 
Also crack growth is more stable in case of [9O45, O45] layup In the case of 
optimal layup as e increases, strength increases Therefore, in order to obtain 
higher fracture strength, adhesive with higher elastic modulus (high value of e) 
for [O45, O45] and optimal layup laminate and lower elastic modulus (low value 
of e) for [O45, 9O45] and [9O45, O45] laminates a-e suggested Variations of 
failure load with the crack length for different e values in [O45, and optimal 
layup [1324g, 1324s] are shown m Figures 2 31 antf and 2 32, respectively 


2.6.5 Overlap length 

Effect of overlap length on the strength of the double lap joint (Figure 2 18 ) is 
studied Strain energy release rate is plotted m Figure 2 33 as a function of crack 
length for various overlap lengths in a double lap joint of aluramium adherend 
and Epon VIII adhesive, (Tkble 27 ) The values of a^s and corresponding 
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Ga’s for different I (= io/A) values are listed in Table 2 20 Failure loads are 
calculated in terms of fracture toughness (Gc) of the adhesive material, for 
various values of I and listed in Table 2 20 As the overlap length increases I 
also increases As I increases strain energy release rate decreases (Figure 2 33a 
and 2 33b) and thus failure load P/ increases 


With an increase m the overlap length I, a* initially increases and then 
decreases But, with an increase in the overlap length I, ag- initially decreases 
and later mcreases The failure load P/ is increased substantially as 1 is increased 
from 0 166 to 0 250 (Table 2 20) As I increases slope of the curves G vs a 
decrease This indicates faster and unstable crack growth and, low failure loads 
are associated with small overlap length Therefore, it is advantageous to use 
large overlap length, which leads to high overlap ratio I, resulting in a joint 
of more strength For various 1 ratios, failure loads are calculated in terms of 
fracture toughness Gc of the adhesive material and listed in Table 2 20 


Failure loads and the crack length variations are studied for the joint of 
carbon fibre reinforced plastics (CFRP) adherends and Epon VIII adhesive 
(Table 27), respectively For various 1 values and laminate combinations the 
predicted failure loads P/s are shown in Table 221 It is observed that as 1 
increases Pf increases appreciably for all combination of 01 ^, 9 O 45 plies, and 
optimal layup For [0^, 01 ^] layup P/ is almost double as compared to rest 
of the layups including optimal layup More stable growth is observed in the 
case of [ 90 ^ 5 , 0^] layup In particular for [01^, and optimal layup 1 
preferably be more than 0 250 Variations of failure load with the crack length 
for different 1 values in [O 45 , O 45 ] and optimal [1324sf» 1324s] laminates are 
shown in Figures 2 34 and 2 35, respectively As also revealed by these figures 
for higher fracture strength, large overlap length (high value of I) is suggested. 
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2.6 6 Adhesive thickness 

Effect of adhesive thickness on the strength of the double lap joint (Figure 2 18) 
IS investigated Strain energy release rate is plotted in Figure 2 36 as a function 
of crack length for various adhesive thickness ratios t (= *./*,) in double lap 
joint of aluminium adherend and Epon VIII adhesive, (Table 27) Failure 
loads and the crack length variations are studied for the joints of unidirectional 
(0°) carbon fibre reinforced plastics (CFRP) adherends and Epon VIII adhesive 
(Table 2 7) With increase in t, adhesive thickness increases and the mesh 
refinement decreases in the adhesive region As the adhesive thickness f* 
increases (t increases), G increases (Figure 2 36) and thus failure load F/ 
decreases Cntical crack length (aa) is observed to decrease with adhesive 
thickness It is also found that vanation of strain energy release rate with crack 
length for t = 0 0125 and t = 0 0100 is nearly same The slope of the curve 
(G vs a) increases as adhesive thickness increases 

It IS also observed that as the crack grows from the location 3 (Figure 2 20c) 
along the interface, strain energy release rate of mode I, initially increases 
to a maximum, then decreases and again increases It can be observed from 
the curves that this local free-edge effect is confined to a small crack length, 
which is of the order of size of the adhesive thickness (t*) For the study of 
effect of adhesive thickness on bond strength and the free-edge effect, a more 
refined mesh consisting of 3376 elements and paired nodes at the interfaces is 
employed From Figure 2 37, it can be observed that there is an increase and 
then decrease in the value of G close to the free edge As the adhesive thickness 
increases (t value increases) this local variation diminishes and is confined to a 
small zone near the free-edge Further, for very small adhesive thickness (t = 
0 0100) this local effect is found to be negligible Associated slope of G vs a 
curve monotonically decreases, as a consequence the crack growth is unstable 

Present analysis indicates faster and unstable crack growth and, lower 
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failure loads are associated with large adhesive thicknesses Therefore, it is 
advantageous to use small adhesive thickness, which leads to low thickness ratio 
t resulting in a joint of high strength However, the value of t below 0 0125 has 
no extra advantage For various t ratios, failure loads are calculated m terms 
of fracture toughness (Gc) of the adhesive matenal and listed in Table 2 22 


For various t values and laminate combinations the predicted failure loads 
P/s are shown m Table 2 23 As t increases the failure load decreases for 
[O 45 , 9 O 45 ] and [901^, O 45 ] layups In all the laminates, there is no significant 
change in the failure load for t < 0 0125 For t = 0 2000 and t = 0 1000 the 
crack growths are found to be stable for [90^5, 90^^], [ 01 ^, 90ls] and [90^^, OI^J 
layup sequence and unstable crack growth is observed in [O^^, OI 5 J layup and 
optimal layup One can therefore, conclude that for higher joint strength, small 
thickness of adhesive (low value of t) is preferable for [OI 5 , [OJ^, 9 OI 5 J, 

optimal layup Variations of failure load with the crack length for different t 

values in [ 01 ^, 0 |g] and [13245, 13245] (optimal layup) are shown in Figures 2 38 
and 2 39, respectively 


Though the nature of G vs a curves are not affected by mesh refinement 

for e and 1 vanations, it is observed that refined mesh is useful for the study of 

free-edge effect particularly for small adhesive thickness Further, it is observed 

that high adhesive stiffness e, high overlap length I and low adhesive thickness 

t are desirable for a strong double lap joint Critical crack length is found 

to be a function of geometrical parameters I. t and independent of material 
parameter e 
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2.7 Conclusions 


On the basis of linear elastic analysis carried out in this Chapter, the following 
conclusions are drawn 


• The finite element approach with paired nodes modelling the crack provides 
a simple computational tool to predict joint failure However, fine gnds 
near the crack-tip are necessary for accurate estimation of stram energy 
release rates 

• Use of Chebyshev’s function is found to be effective for automatic mesh 
generation with various refinements at critical zones 

• Stress distributions along the overlap length of the single lap joint are 
predicted by the present computer code and the results agree with the 
those reported by Hiregoudar (1993) 

• Stress distnbutions are also predicted along the overlap length and across 
the adhesive thickness of thick adherend shear test (TAST) specimen and 
the results are observed to agree with the those reported by Su and Mackie 
(1993) 

• In case of single edge notch (SEN) specimen, with 2688 elements and with 
fine mesh near the crack-tip, error estimated m the value of strain energy 
release rate is less than 1 6 percent (as compared to close form solution) 

• Strain energy release rate computed by the present finite element analysis 
with fracture mechanics approach agrees well with the experimental result 
for cracked lap shear (CLS) specimens (Fernlund and Spelt, 1991c) The 
diflference is less than 2 7 percent. 
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Cracks located at the edges are more critical than those at the central 
part of the overlap region Further, in symmetric boundary conditions the 
diagonally opposite cracks are equally prone to failure 

Introduction of 0® ply m the immediate neighbourhood of adhesive results 
in a stronger joint In addition, a proper choice of layup sequence results 
in an efficient bonded joint 

Use of adhesives of high Young’s moduli results in stronger bonded joint 

Employment of longer overlap length results in stronger bonded joint 

Use of small thickness of adhesive results in stronger bonded joint 

For bonded joints with very thin adhesives free-edge effects are found to 
be significant Moreover, for the analysis of bonded joints having very 
small adhesive thicknesses, more refined meshes are preferred to account 
for the free-edge effects 
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Seven fracture criteria considered in the finite element analysis 


Criterion type 

Expression 


1 

Gi + Gn 

> 

Gic + Gnc 

2 

Gi 

> 

Gic 

3 

Gjj 

> 

Gnc 

4 

Gi + Gn 

> 

Gic 

5 

G] + Gji 

> 

Gnc 

6 

Gi Gn 

Gic Gnc 

> 

1 

7 

\Gic) \Gnc) 

> 

1 
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Table 2 2 Specifications of single lap joint (Hiregoudar, 1993) 


Single lap joint 

overlap length (£„) 
total length (Lt) 
applied stress (a) 
adhesive thickness (i,) 
adherend thickness (tr) 

Adhesive (FM 73) 


Young's modulus (Es) 2 21 Gpa 

Poisson's ratio (i/,) 0 43 

yield stress ) 40 Mpa 

Adherend (Steel) 

Young's modulus (Er) 211 3 GPa 
Poisson's ratio (i/^) 033 

yield stress (cto ) oo 


16 mm 
80 mm 
200 MPa 
030 mm 
1 60 mm 
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Table 2 3 Specifications of stepped lap joint (thick adherend shear test specimen 
of Su and Mackie, 1993) 


Joint configuration 


overlap length (L^) 8 mm 

total length (Lt) 125 mm 

applied load (P) 3 5 kN 

total joint thickness (t) 13 35 mm 

width of the joint (W) 25 4 mm 


Adhesive (rubber-toughened epoxy) 

Young’s modulus (Es) 6 60 Gpa 

Poisson’s ratio (i^s) 0 38 

thickness of adhesive (**) 0 65 mm 


Adherend (Steel) 

Young’s modulus (Er) 210 Gpa 

Poisson’s ratio (ur) 03 

thickness of adherend (tr) 6 35 mm 


Table 2 4 Physical dimensions and matenal properties of SEN specimen 
(Figure 2 15) 


length (L) 

100 mm 

width (W) 

50 mm 

thickness 

1 mm 

crack length (a) 

2 mm 

Young’s modulus (E) 

70 GPa 

Poisson’s ratio (u) 

03 

stress applied (tr) 

01 GPa 
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Table 2 5 Physical dimensions and material properties of CLS specimen 
(Figure 2 16) 


Adhesive (Cybond 4523 GB) 


Young’s modulus 

Poisson’s ratio 

thickness 

width 

length 


8 GPa 
0 35 

0 21 mm 
10 mm 
180 mm 


Adherend (Aluminium 6061~T6) 


Young’s modulus 

Poisson’s ratio 

thickness 

width 

lengths 

load (Pf) 


70 GPa 
03 

3 06 mm 
10 mm 

180 mm and 240 mm 
6 33 N/mm 
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Table 2 6 Specifications of the adhesively bonded joint (Figure 2 18) 


Joint specifications 

total length (L^) 80 mm 

overlap length (Xo) 40 mm 

load applied (P) 1000 N/mm 

width (W^) 1 mm 

Adhesive (Epon VIII) 

Young’s modulus {Eg) 3 5 GPa 

Poisson’s ratio (i/,) 0 41 

thickness (t,) 2 mm 

Adherend (Aluminium) 

Young’s modulus {E^) 70 GPa 

Poisson’s ratio {vr) 0367 

thickness (f,.) 10 mm 
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Table 2 7 Specifications of double lap joint 


Double lap joint 

overlap length (La) 40 mm 

total length (Lt) 80 mm 

Adherend (Aluminium) 

Young’s modulus (Er) 70 GPa 

Poisson’s ratio (j/^) 0 3 

adherend thickness (tr) 10 mm 

Adherend (carbon fibre reinforced plastic) 


longitudinal modulus (£7,) 137 GPa 

transverse modulus (Et) 7 GPa 

shear modulus Gu 45 GPa 

Poisson’s ratio (uu) 0 3 

adherend thickness (t^) 10 mm 

Adhesive (Epon VIH) 

Young’s modulus (Es) 3 5 GPa 

Poisson’s ratio (i/,) 0 41 

thickness (t,) 2 mm 

fi-acture toughness (Gc) 0 366 kj/m^ 


fracture criterion Gj + Gu > Gc 
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Table 2 8 Description of eight different crack locations and their growth 
direction (Figure 2 20c) 


Crack Location 

Interface 

From 

To 

1 

bottom 

left edge 

center 

2 

bottom 

right edge 

center 

3 

top 

right edge 

center 

4 

top 

left edge 

center 

5 

bottom 

center 

left edge 

6 

bottom 

center 

right edge 

7 

top 

center 

left edge 

8 

top 

center 

right edge 


Table 2 9 Computed critical strain energy release rates for two different types 
of boundary conditions and eight different crack locations (Figure 2 20a-c) 


crack growth 
location 

free side 

acr(mm) Gcr(kJfm^) 

constrained side 
acr(mm) GcrikJ/ir?) 

1 

172 

75 

172 

14 2 

2 

098 

05 

098 

05 

3 

172 

19 5 

172 

14 2 

4 

098 

57 

098 

05 

5 

21 

62 

21 

65 

6 

21 

62 

21 

65 

7 

21 

58 

21 

65 

8 

21 

58 

21 
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Table 2 10 Material properties of metal adherends 



aluminium 

brass 

copper 

steel 

Young’s modulus (GPa) 

70 

90 

120 

210 

Poisson’s ratio 

03 

0 25 

0 25 

025 
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Table 2 11 Material properties of composite adherends 


glass fibre reinforced plastic 

longitudinal modulus (Ei) 39 GPa 

transverse modulus (£?<) 10 GPa 

shear modulus Gu 7 GPa 

Poisson’s ratio (u^) 0 25 

kevlar fibre reinforced plastic 

longitudinal modulus (Ei) 83 GPa 

transverse modulus (Et) 5 6 GPa 

shear modulus Gu 2 1 GPa 

Poisson’s ratio (i/^) 0 34 

carbon fibre reinforced plastic 

longitudinal modulus (Ei) 137 GPa 

transverse modulus (Et) 7 GPa 

shear modulus Gu 4 5 GPa 

Poisson’s ratio (u^) 0 3 

graphite fibre reinforced plastic 

longitudinal modulus (Ei) 159 GPa 

transverse modulus (Ei) 10 9 GPa 

shear modulus Gu 6 4 GPa 

Poisson’s ratio (ui) 0 38 



Table 2 12 Specifications of double lap joint 


Double lap joint 

overlap length (Lo) 120 mm 

total length (Lt) 240 mm 

lower adherend thickness (f^) 2 mm 

Adhesive (Epon VIII) 

Young’s modulus (JSg) 3 5 GPa 

Poisson’s ratio (i/*) 0 41 

thickness (t*) 0 1 mm 

fracture toughness (Gc) 1 11 kJ/m^ 

fracture cnterion Gj + Gn > Gc 
where Gc ~ Gjc + Guc 
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Table 2 13 Computed failure loads in N/mm for joints of aluminium and vanous 
metals 


tr(upper) 

tr(lower) 


crack locations 
1 2 3 


4 



aluminium to aluminium 


05 

144 99t 

26499 

212 49 

148 74 

1 

249 99 

212 49 

167 49t 

242 49 

2 

459 99 

18249 

141 24t 

429 99 

4 

849 99 

163 74 

129 99t 

774 99 


aluminium to copper 


05 

129 99 

42249 

339 99 

126 24t 

1 

212 49 

32499 

26124 

197 49t 

2 

369 99 

26874 

212 49t 

324 99 

4 

654 99 

23499 

182 49t 

564 99 


aluminium to brass 


05 

137 49t 

324 99 

264 99 

137 49t 

1 

234 99 

257 49 

204 99t 

219 99 

2 

414 99 

216 24 

171 24t 

377 49 

4 

744 99 

193 74 

152 49t 

669 99 


aluminium to steel 


05 

118 74 

699 99 

579 99 

113 12t 

1 

186 25 

534 99 

429 99 

167 49t 

2 

30249 

41499 

32499 

264 99t 

4 

519 99 

347 49 

272 49t 

429 99 


t 


minimum failure load among the four crack locations 
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Table 2 14 Computed failure loads in N/mm for joints of steel and various 
metals 


tr (upper) 
trQower) 


crack locations 
1 2 3 


4 


steel to aluminium 


05 

324 99 

197 49 

159 99t 

339 99 

1 

594 99 

174 99 

141 24t 

609 99 

2 

1134 99 

163 74 

129 99t 

1134 99 

4 

— 

156 25 

124 37t 

— 


steel to copper 


05 

272 49 

294 99 

242 49t 

272 49 

1 

474 99 

253 74 

204 99t 

467 49 

2 

879 99 

227 49 

182 49t 

849 99 

4 

— 

216 25 

171 24t 

— 


steel to brass 


05 

294 99 

234 99 

197 49t 

30249 

1 

53499 

20874 

167 49t 

53499 

2 

999 99 

189 99 

152 49t 

98499 

4 

— 

18249 

144 99t 

— 


steel to steel 


05 

24249 

474 99 

399 99 

234 99t 

1 

414 99 

399 99 

324 99t 

384 99 

2 

744 99 

347 49 

272 49t 

654 99 

4 

137499 

31749 

242 49t 

1194 99 


t minimum failure load among the four crack locations 
— failure load is > 1000 N/mm 
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Table 2 15 Computed failure loads m N/mm for joints of carbon fibre reinforced 
plastic (CFRP) composite and various metals 


tr(upper) 
tr (lower) 

1 


crack locations 
2 3 4 


CFRP to aluminium 


05 

234 99 

193 74 

156 24t 

208 74 

1 

399 99 

163 74 

133 74t 

354 99 

2 

74499 

148 74 

122 49t 

654 99 

4 

1419 99 

14499 

116 87t 

125499 


CFRP to copper 


05 

29499 

234 99 

197 49t 

302 49 

1 

534 99 

208 74 

167 49t 

534 99 

2 

999 99 

189 99 

152 49t 

984 99 

4 

-• 

18249 

144 99t 

— 


CFRP to brass 


05 

20874 

23499 

189 99 

186 24t 

1 

354 99 

197 49 

159 99t 

309 99 

2 

624 99 

17499 

152 49t 

984 99 

4 

— 

18249 

144 99t 

— 


CFRP to steel 


05 

159 99 

47499 

384 99 

144 99t 

1 

249 99 

369 99 

30249 

219 99t 

2 

40749 

30999 

257 49t 

354 99 

4 

729 99 

279 99 

234 99t 

624 99 


t 


minimum failure load among the four crack locations 
-- failure load is » 1000 N/mm 
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Table 216 Computed failure loads in N/mm for joints of carbon fibre reinforced 
plastic (CFRP) composite and various fibre reinforced composites 


tr (upper) 
tr(lower) 

1 


crack locations 
2 3 4 


CFRP to glass 


05 

362 49 

12812 

99 99t 

377 49 

1 

654 99 

11312 

88 74t 

684 99 

2 

1149 99 

105 62 

8312t 

1269 99 

4 

131499 

103 74 

8124t 

1449 99 


CFRP to kevlar 


05 

29124 

223 74 

167 49t 

339 99 

1 

489 99 

18249 

141 24t 

579 99 

2 

83499 

159 99 

124 37t 

999 99 

4 

1449 99 

14874 

114 99t 

1449 99 


CFRP to CFRP 


05 

27249 

208 74 

163 74t 

178 74 

1 

45249 

163 74 

129 99t 

302 49 

2 

77499 

13749 

111 24t 

527 49 

4 

1314 99 

12437 

99 99t 

924 99 


CFRP 

to graphite 


05 

249 99t 

369 99 

29499 

249 99t 

1 

429 99 

287 49 

231 24t 

414 99 

2 

729 99 

24249 

193 74t 

729 99 

4 

1299 99 

21624 

174 99^ 

1299 99 


t mmimum failure load among the four crack locations 
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Table 2 17 Computed failure loads (P/’s) for various combinations of layups 
in different crack growths 


Crack 

location 

OI5] 

Failure loads in N/mm 
[90S,, 90S,] [OS,, 90S,] (9QS,, OS,] 

[13245, 13245) 

1 

45 0 

26 0 

242 0 

15 0 

63 2 

3 

28 5 

15 0 

12 0 

48 0 

26 8 

4 

1QO/1 

1300 

fno iA CO /Oj 

420 

rko / >icoi 

162 0 

42 0 

56 8 


1324 [O74579OV - 45°] 


Table 2 18 Predicted failure loads of adhesively bonded joints with different 
Young’s moduli of the adhesives (L7Lt= 1/2, 1/5, no of elements 

employed = 1816) 


Es/Er 

a*(mm) 

Ga,(kJ/l 

0 0125 

172 

22 4 

00250 

191 

14 9 

00500 

21 

10 5 

01000 

2 42 

78 


.cr(mm) 

G„ikJ/m^) Pf/VGc(kN/m) 

172 

293 

184 74 

172 

28 2 

18831 

172 

19 2 

228 21 

172 

145 

262 61 
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Table 2 19 Computed failure loads (P/’s) for various adhesive stiffnesses and 
different combinations 


EJB, 

[O4S1 

Failure loads in N/mm 
[90;s,90Jsl 10Js,90SJ [905s, 0 ;,] 

[1324s, 1324s] 

00125 

215 

15 0 

12 0 


20 0 

0 0166 

23 5 

15 5 

13 0 

16 75 

22 3 

0 0250 

26 0 

15 5 

12 75 

16 25 

24 0 

00500 

28 5 

15 0 

12 0 


26 8 

01000 

29 5 

14 0 

no 

13 75 

28 0 


1324 [07457907-45°] 


Table 2 20 Predicted failure loads of adhesively bonded joints with different 
overlap lengths {tjty- 1/5, Es/Er= 1/20, no of elements employed = 1816) 


LojLt 

aj;(mm) 

G^(kJ/m^) 

acr(mm) 

Gcr(kJ/rr^) 

PfJVGcikN/m) 

0 083 

mm 

350t 

137 

610 

40 48 

0125 

26^ 

250t 

132 

452 

47 03 

0166 


200t 

111 

312 

56 61 

0 250 

235 

50 

132 

48 

144 33 

0 500 

215 

10 

172 

19 2 

228 21 


t represents actual value is greater than the specified one 






78 


CHAPTER 2 ELASTIC ANALYSIS 


Table 2 21 Computed failure loads (P/’s) for various overlap lengths and 
laminae combmations 


LolU 

[01^, 0|g] 

Failure loads in N/mm 
(90|s,90|sl [0;s,90Js] [90;s,0|sl 

(1324s, 1324sl 



6 25 

75 

9 62 

117 


12 5 

80 

8 87 

115 

140 



no 

105 

14 0 

180 


285 

15 0 

120 

15 0 

28 5 

1324 

[O7457907 - 45°] 




ble 2 22 Predicted failure loads of adhesively bonded joints with different 

icknesses of the adhesives {LoJLt- 

1/2, EsfEr = 

= 1/20, no of elements employed 

3376) 






tshr 

Oj(mm) 

G^{kJ/m^) 

acr(nim) 


Pil^G^kN/m) 

00100 

016 

41 

003 

6 85 

382 072 

00125 

019 

42 

013 

78 

358 049 

0 0250 

032 

45 

023 

8 65 

321 908 

0 0500 

055 

54 

036 

10 2 

313 110 

01000 

100 

70 

067 

12 7 

280 607 

0 2000 

193 

117 

1 1 

20 4 

221 396 
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Table 2 23 Computed failure loads (P/’s) for various adhesive thicknesses and 
laminae combinations 


ts/tr 

[®4S5 ^45] 

Failure loads in N/mm 
[90i”s.90Ssl [0«,90S,,1 r9OS.„0Jsl 

[1324s, 1324 s1 

0 0100 

39 0 

16 0 

1175 

13 25 

410 

00125 

380 

16 0 

1175 

13 25 

400 

00250 

36 5 

15 5 

1125 

13 25 

36 0 

0 0500 

34 0 

15 5 

1125 

14 0 

33 0 

01000 

32 0 

15 5 

115 

14 75 

310 

0 2000 

28 5 

15 0 

120 

15 0 

28 5 


1324 [ 07457907 - 45 °} 
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Figure 2 1 Finite element mesh at crack-tip illustrating the crack closure 
technique 
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Y V 



Figure 2 2 Schematic representation of crack closure technique 




Figure 2 3* Flow chart of linear elastic finite element analysis for crack growth 
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I 



Figure 2 4 Variation of failure load and crack length for different possible 
fracture criteria (Table 21) 
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Figure 2 5 Schematic of single lap joint considered by Hiregoudar (1993) and 
the present finite element mesh consists of 1816 elements with an enlarged view 
of the edge region 
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Figure 2 6 Comparison of stress distnbutions along the overlap length of the 
single lap joint considered by Hiregoudar (1993) 
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Figure 2 7 Schematic of the stepped lap joint (thick adherend shear test TAST 
specimen of Su and Mackie, 1993) and the finite element mesh with an enlarged 
view of the stepped portion 
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Figure 2 
interface 




Comparison of stress distribution along the overlap length at 
adherend and adhesive of the stepped lap joint 




87 


Figure 2 9 
interface of 



Overlap Length 



Comparison of strain distnbution along the overlap length at 
adherend and adhesive of the stepped lap joint 
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Figure 2 10. Companson of stress distnbution along the overlap length at 
midline of adhesive of the stepped lap joint 
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Figure 2 12 Normal stress and strain distributions along the overlap 
length and across the adhesive thickness of the stepped lap joint 















Error In G|C/ ) 
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Pf =633N/mm, Width =10mm 



I — 3 06 mm 
j— — I -]ri. 0 21 mm 
^3 QS mm 


Pf 



representation of cracked lap shear (CLS) specimen 
r point bending and the finite element mesh consists of 2180 elements 


STRAIN ENERGY RELEASE RATE(kJ/m 
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Figure 2 17 Variation of strain energy release rates with crack length for plane 
stress and plane strain models of the cracked lap shear (CLS) specimen 
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Figure 2 18. Schematic of adhesively bonded double lap joint and the finite 
element mesh showing the crack location at ’3’ 
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Figure 2 21* Variation of strain energy release rate with crack length for crack 
growths at eight different locations (Figure 2 20c) for free side type boundary 
condition 
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Figure 2 22 Variation of strain energy release rate with crack length for crack 
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Figure 2 23 Schematic of double lap joint of composite adherends [1324s, 
1324s] with a magnified section view 
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Crack Length (mm) 



Figure 2 24 Variation of failure load with crack length for symmetric combinati- 
ons of 0°,45°,90° and -45® layups in composite adherends, and for crack 
location at ’3’ 
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Figure 2 25 Variation of failure load with crack length for various symmetnc 
combinations of 0°,45°,90° and -45“ layups m composite adheren-ds, with 0 
layups being present m the immediate nei^bourhood of the adhesive and, for 
crack location at ’3’ 
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Figure 2 26 Vanation of failure load with crack length for [O45, Ol^}, [9O45, 9O45], 
[O45, 9O45], [90^5, 0|g] laminate pairs for crack location at T’ 



Crack length (mm) 

Figure 2 27 Variation of failure load with crack length for [O45, O45], [9O45, 90^^ ), 
[O45, 9O45], [9O45, O45] laminate pairs for crack location at ’3’ 
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Crack length (mm) 

Figure 2 28 Vanation of failure load with crack length for [O 45 , 0^], [ 9 O 45 , 904 ^], 
[O45., 9O45], [9O45, O45] laminate pairs for crack location at ’4’ 



Crack length (mrn) 

Figure 2 29 Vanation of failure load with crack length for most likely interface 
edge crack growths in CFRP unidirectional and 9015^) laminate pairs 
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Figure 2 30 Variation of total strain energy release rate with crack length for 
crack initiation at location ’3’ and different e (= values in aluminium to 

aluminium joints 
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Figure 2 31 Variation of failure load with crack length for crack initiation at 
location ’3’ and different e (= Eg/jEr) values m [O 45 , O^g] CFRP laminates 
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Figure 2 32 Variation of failure load with crack length for crack initiation at 
location ’3’ and different e (= values in [13245, 1324s] (optimal) CFRP 

laminates 
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Figure 2 33 Variation of total strain energy rele^e rate with 

crack initiation at location ’3’ and different 1 (= values in aluminium 

aluminium joints 
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Figure 2 34 Vanation of failure load with crack length for crack initiation at 
location ’3’ and different 1 (= Lo/LO values in [Ol^, O^g] CFRP laminates 




Ill 



CRACK LENGTH (mm) 


Figure 2 35 Variation of failure load with crack length for crack ^ 

location ’3’ and different 1 (= values in [1324s, 1324s] (optimal) CFRP 

laminates 
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Figure 2 36. Variation of total strain energy release rate with crack lencth for 
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Figure 2 37: Variation of total strain energy release rate with crack length for 
crack initiation at location ’3’ and different t (= values in aluminium 
aluminium joints with very fine mesh refinement 
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Figure 238 
location ’3' 


Variation of failure load with crack length for crack initiation at 
and different t (= values in [01^, Ol^] CFRP laminates 
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Figure 2 39 VariEtion of failure load with crack length for crack initiation at 
location ’3’ and different t (= values in [132%, 132%] (optimal) CFRP 
laminates 



Chapter 3 


Elasto-plastic analysis 


3.1 Introduction 


Adhesive materials play a major role in the load transfer mechanism of bonded 
jomts These are mostly polymenc materials and exhibit material non-lmearities 
em VIZ elasto-plastic and visco-plastic Therefore, it is important to consider 
the non-hnear behaviour of the adhesive for a more realistic analysis In 
this Chapter, adhesive is considered to be elasto-plastic in behaviour, and 
adherend is considered to be linear elastic Finite element code with material 
non-lmeanty (elasto-plastic behaviour), is developed and validated with existing 
results from the literature Elasto-plastic analyses of double lap joint are carried 
out Growths of plastic zones are studied In addition, influence of various 
parameters, such as (i) stacking sequence in composite adherends, (ii) fracture 
toughness of the adhesive, (m) overlap length of the joint, and (iv) thickness 
of adhesive on the strength of adhesively bonded joints, are investigated 
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3.2 Formulation 


Yamada., et al (1968) presented a classical solution for continuum elasto-plastic 
problems by means of finite element approach This method is based on a plastic 
stress-strain matrix, which is derived by inverting the Prandtl-Reuss equations 
of theory of plasticity Zienkiewicz, ct al (1969) suggested an imtial stress 
computational process, which showed rapid convergence than other alternative 
approaches This also permits large load increments without violating the yield 
criterion They presented a general formulation for the elasto-plastic matrix 
for evaluating stress increments with an associated flow rule Several solutions 
showing stress distribution, strain development and growth of plastic enclaves are 
given for both the von Mises and Coulomb (Drucker) type yield surface Owen 
and Hinton (1980) explained in detail the finite element approach to material 
non-hnearities (elasto-plastic and visco-plasticity) Whitcomb and Dattaguru 
(1984) developed a general computer code for the geometrical and material 
non-lmear analyses of structures (GAMNAS) 


3.2.1 Constitutive relations 

In non-linear elastic behaviour, the stresses can be expresses as a function of 
strain as <t = a(e) The main difference of plastic formulation from non-hnear 
elastic formulation is that such a unique relationship between stress and strain 
IS not available Although the stresses at any level of strain have to he on or, 
within the current yield surface, it is difiicult to determine the exact value of 
each component 

Theoretical as well as practical mterest in plasticity has stimulated much 
work m the numerical solution of elasto-plastic problems by finite element 
method and other procedures Small strain displacement elasto-plastic solution 
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using finite element method is discussed by Yamada, et al (1968) The 
difference m elasto-plastic analyses reported by other researchers lie between 
the apparent forms of non-lmear solution algorithms used, and the form of 
constitutive relations postulated Some common computational approachs used 
with incremental theory of plasticity m finite element method are (i) incremental 
tangent stiffness, (ii) initial strain and (iii) initial stress methods 

In incremental tangent stiffness method the change of material constitutive 
matrix forming element stiffness, assembling and inverting an overall stiffness 
matrix for each step of analysis is earned out In the initial strain and the 
initial stress approaches, in all the iterations a simple resolution of constant 
linear elastic problem is earned out The initial stress method is found suitable 
for solving problems with ideal plasticity or with strain hardening 


I Yield entena 

A yield criterion is a hypothesis concerning the limit of elasticity under any 
possible combination of stresses The suitability of any proposed yield criterion 
should be vahdated by conducting experiments Yielding occurs when some 
scalar function of the principal stress differences reaches a critical magnitude, 
mathematically it can be expressed as 


/((Ti — <T 2 , cr 2 — <73, <73 — (Ti) — Constant 


IWo different yielding criteria namely (i) Tresca and (ii) von-Mises criteria 
are described in brief 

Tresca yield criterion Perhaps the simplest function, which satisfies the 
Equation (3 1) is of the form 
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a^ — CTj = con stant 


(3 2) 


Above equation can be interpreted as, when largest of the three stresses 
- cTz), (az - crz), and (crs - aj) attains a critical constant value, (for a given 
material) yielding occurs 

von Mises yield criterion Another admissible function, which satisfies 
Equation (3 1) is of the form 


(<Ti — crz)^ + (<72 — crs)^ + (<73 — CTi)^ = constant 


(3 3) 


In this type of function, each of the principal stresses contributes to yieldmg 
rhis criterion can also be interpreted as, when the shear strain energy reaches 
1 critical constant value, yielding occurs 


y / ((<rx - <7j,) 2 + (cTy - <7^)2 + (<7* - <7*)2 + + 7^J)/2 0 - <70 


(3 4) 


I Stress-strain relation 

iooke’s law for an isotropic elastic matenal can be written as 


(3 5) 


{<7} = P]{e} 
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where [D] is the matenal matrix An equivalent relation m plasticity is 
dependent on present state of stress and deformation Yamada, et al (1968) 
have developed a plastic stress-strain relation and presented in a matrix, which 
IS termed as Elasto-plastic matrix This matrix is derived by inverting Prandtl- 
Reuss equations obeying von Mises yield cnterion For a two dimensional plane 
strain model this matrix is written as 


[D,\ 

E 

l-\-u 

l-y _ ZsA 
l-2u S 

/ / 

-Z 

l-2«/ S 

/ t 

^ 

l-2u S 

l~2z/ S 

$ 

t 

S 

(3 6) 



/ 

L ^ 

1 

.... 

S 

1 

1 


Pnmes represent the deviatonc 

stress components and 


S = 

2d^ Et 

3 ^ ^ 3i! 

L 2(l+«/)J 



(3 7) 


where, Er and 5- are the tangent modulus and effective stress, respectively 

According to von Mises criterion, effective stress for two dimensional problems 
is written as 


a- > <ro 
a 




where, <ro represents yield stress of the matenal 


(3 8) 
(3 9) 
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3.2.2 Finite element technique 


Similar to the computer codes developed by Owen and Hinton (1980) and 
Whitcomb and Dattaguru (1984), a new code is written to predict the failure 
load, which is explained in the flow diagram given in Figure 3 1 The module, 
which takes into account the non-lmear (plastic) behaviour of the material is 
indicated by the dashed portion m the Figure Details of the elastic and plastic 
analyses are explained in the main algorithm in Section 2 2 and plastic algonthm 
of this Section, respectively 


I Non-lmear plasbc algonthm 


step 01 


step 02 : 

step 03 : 


compute the followings 
Dep = D^(S-) 

B” = Bo + Sjvi(6") 

n 

calculate proportionality factor H 

O-n+l _ 

- 

CTq - a" 

compute residual stress to be compensated 


step 04 : update stress 

step 05 • calculate residual forces 

A'tpep = + Af 

S7 

where 

Af represents nodal forces 
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step 06 check for convergence 
< tolerance 

if converged go to step 07 of mom olgontKm in Section 2 2 
else continue 

step 07 solve for displacement 
AS = [K]'^AP 

step 08 calculate mcremental stresses and strains 
Ac = [B"f A6 
Aa = [Dep]Ae 

step 09 update displacements, strains and stresses 

^n+l = ^ ^ 

gn+l = gn + ^gn 

g^n+l - + /io-" 

Step 10 , go to step 01 


As mentioned instep 01, S” and isf” are functions of the displacement 
vector of iteration One is expected to update these values in eveiy iteration 
But there is a possibility of divergence in results using these values However, 
using constant values of these parameters and Newton Raphson’s modified 
approximation (Owen and Hinton, 1980 and Zienkiewicz and Taylor, 1991) 

convergence is assured In addition, this reduces substantially the computational 
time and memoiy 

K the effective stress a oi any element in the adhesive, exceeds the yield 
v^ue of the matenal o-„, then the material behaves elasto-plastic and clasto- 
plastic analysis is to be earned out Module, which takes care of plastic part 
of elasto-plastic analysis, is indicated by dashed box of Figure ^ 1 
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Elasto-plastic calculations are iterated till the elasto-plastic residual (out-of- 
balancc) force ^ij>rp becomes smaller than a prescribed tolerance value Strain 
energy release rates arc calculated and, as mentioned earlier in Chapter 2 and 
the fracture criterion Gj + Gji > Gic + Guc is employed to predict the failure 
load A comparison of these loads, for linear elastic model and plastic models 
with various tolerance values, and the effect of strain hardening are discussed 
below 


I Convergence entenon 


Elasto-plastic analysis predicts a relatively higher failure load m comparison 
with that of linear elastic analysis This might be due to the fact that elasto- 
plastic analysis results in a considerable reduction in peak stresses as compared 
with elastic analysis This agrees with the finite element solutions obtained 
by Thomsen (1992) The failure load increases with the decrease m tolerance 
value and saturates for a tolerance value ~ 0 15 percent (Figure 3 2) 


I Strain hardening 

Figure 3 3 represents the strain hardening behaviour of the adhesive matenal 
Predicted failure loads remains same for strain hardening parameter k is in 
the range 0 2 to 0 8 Hart-Sraith (1973) also observed that arbitrary bi-linear 
representations leads to the same result It can be concluded that, little accuracy 
could be be gained by more complex representation of adhesive characteristics 
However, in the elasto-plastic analysis size of plastic zone is reduced rapidly, 
with the increase in hardening parameter k 
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3.3 Validation of computer code 


In order to validate the code developed the two-dimensional elasto-plastic beam 
problem considered by Zienkiewicz, et ol (1969) is analysed using plane strain 
CST element undergoing small displacement Young's modulus, Poisson’s ratio 
and yield stress of the beam are considered to be 70 GPa, 02 and 243 0 
MPa, respectively Figure 3 4 shows the comparison of the results using 300 
isoparametnc quadrilateral four noded elements, 75 CST elements and 150 CST 
elements with that of 200 triangular elements employed by Zienkiewicz, et al 
(1969) It is observed that the computed load-displacement response with 150 
CST elements are in fairly good agreement with that reported by Zienkiewicz, 
ct al (1969) 


3.4 Growth of plastic zone 


The growth of plastic zone in the adhesive is investigated for three different 

displacement boundary conditions Further this growth is studied with crack 
propagation 


3.4.1 Boundary conditions 


Adhesively bonded joints with aluminium adherend and Epon VIII (Adams 
and Wake, 1984) adhesive (Table 3 1) are chosen for the study of the growth 
of plashc zone tor various boundaiy conditions Growth of plastic zones for 
fte boundary conditions of types- I, 2 and 3 are illustrated in Figures 3 5a, 
35b and 3 5c, respectively It is observed that the growths of plastic zones 
are sensitive to the displacement boundary conditions For type-l boundaiy 
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condition the growths of plastic zones are along the diagonal joining the left 
top and right bottom corner of the adhesive layer While for types- 2 and 3 
boundaiy conditions the growths of plastic zones are significantly different It 
IS along the left bottom corner to the midpoint of the top adherend-adhesive 
interface In other words for type-1 boundary condition failure is cohesive type 
(limited within the adhesive only) While for types- 2 and 3 boundary conditions 
failure is adhesive type (not limited to the adhesive only, interface between the 
adherend and adhesive is also involved) Further, m type-1 boundaiy condition 
the adhesive is under pure shear leading to stronger joint, as compared to other 
boundary conditions where, peel stresses are present 


3.4.2 Crack propagation 


Double lap joint (Figure 3 6) with carbon fibre reinforced plastic CFRP adherend 
and Epon VIII (Adams and Wake, 1984) adhesive (Tkble 3 2) is considered to 
study the growth of plastic zone associated with crack propagation As explamed 
earlier crack propagation is modelled by sequential opening of paired nodes 
at the interface of adherend and adhesive Figure 37 illustrates the growths 
of plastic zones with the crack propagation Crack location ’3’ is considered 
Size of plastic zone decreases faster as the crack propagates and plastic zone 
altogether disappears when the crack becomes large in size 


3.5 Parametric study 


Effects of stactang sequence, critical strain energy release rate of the adhesive, 
overlap length and thickness of the adhesive on the strength of the double lap 
joint (Figure 3 6) are investigated Because of symmetry only one half of the 
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joint IS considered for the analysis The problem is treated as one of plane 
strain model A finite element mesh consists of 1816 CST elements is employed 
for the analysis 

Three non-dimensional parameters are defined viz 


(3 10) 

(311) 

(312) 

where, C^c and Gc are the critical strain energy release rates of the adhesive 
used and the standard adhesive (epoxy), respectively and Lt denote the 
overlap length and the total length of the bonded joint In addition tg and L 
represent the thicknesses of adhesive and adherend, respectively C^q is assumed 
to be in multiples of Gc In addition, assuming constant critical strain energy 
release rate Gc of the standard (epoxy) adhesive, adherend thickness L and 
overlap length L„ parametnc study is carried out The adherends and adhesive 
of the joint, are considered to be carbon fibre reinforced plastics (CFRP) and 
Epon VIII (Table 3 3), respectively Crack located at ’3’ is found to be critical m 

subsection 2 6 2 Therefore, only crack location at ’3’ is considered for further 
discussion m this Chapter 



3 5.1 Stacking sequence 


A symmetnc quaswsotropic combination is considered for the analysis to 
investigate the effect of stacking sequence on the bond strength Most 
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widely used fibre orientations 0°, 45°, 90° and -45° are considered Among 
all the possible symmetric combinations of these plies, it is observed in 
subsection 2 6 3 that [0°/45°/90°/ - 45°]s (optimal) layup, has the maximum 
bond strength among its combinations Also the presence of 0°, ply m 
the immediate neighbourhood of the adhesive results in stronger joint than 
the plies of any other orientation Elasto-plastic analysis is earned out 
for [0°/0°/0°/0°]s, [OV45°/90°/ - 45°]s and [0°/0°/90°/90°]s layups of CFRP 
adherends It is found that employing [0°/0°/90°/90°]g layup results in a 
stronger joint is obtained than [0°/0°/0°/0°]s layup This combination is found 
to be a stronger joint than [0°/45°/90°/ - 45°]^ layup The failure loads are 
computed for various adhesive critical strain energy release rate g (= G'c/Gc) , 
overlap length 1 (= Lo/Lt)&nA adhesive thickness t (= tg/tr) ratios for these 
three layups Detail of these results are discussed in the respective subsections 


3.5.2 Fracture toughness of adhesive 

For various g (= C?^/C?c) values and layups the computed failure loads are 
listed in Table 3 4 It is observed that as g increases failure load increases for 
all the three layup sequences considered Adhesive with higher value of critical 
strain energy release rate ( high g value) leads to a stronger joint Therefore, in 
order to achieve a higher fracture strength of the joints, adhesives with higher 
critical strain energy release rates are suggested Variation of failure load with 
crack length for various g values in [0°/0°/0°/0°]s, [0°/45°/90°/ ”45°]5 (optim^) 
and [0°/0°/90°/90°]s CFRP laminates are listed in Figures 3^, 3 9 and 340, 
respectively 

Though fracture toughness of adhesive is generally denoted by Kw, the 
stress intensity factor for ideal plane strain adhesive material, in the present 
study, for convenience critical strain energy release rate Gc is also termed as 
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fracture toughness 


3 5.3 Overlap length 

For various values of 1 (= LofLt) ^.nd the three layups chosen the computed 
failure loads are listed in Table 3 5 It is observed that as I increases failure 
load increases for all the three layup sequences For higher fracture strength of 
the bonds longer overlap lengths (higher l^^ar^su^ested Variation of failure 
load with crack length for vanous I values^in^j^^5“/90°/ - 45"]*, (optimal) 
and [0‘’/0°/90V90“]5 CFRP laminates are shown m Figures 3 11, 3 12 and 3 13, 
respectively 


3.5 4 Adhesive thickness 

For various t (= tg/t,) values and layups the computed failure loads are listed 
in Table 3 6 It is observed that as t increases failure load decreases for all 
the three lajmp sequences For higher fracture strength of the bonded joints, 
small thicknesses of the adhesives ®u88®®ted Variation of failure 

load with crack length for various t values m^[0°/45V90“/- 45“ ]s (optimal) and 

[0“/0“ /90°/90“]s CFRP laminates are illustrated in Figures 3 14, 3 15 and 3 16, 
respectively 


3.6 Conclusions 


On the basis of elasto-plastic analysis earned out in this chapter, the following 
conclusions are drawn 
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• The cldstic analysis predicts lower failure loads as compared to the elasto- 
plastic analysis Therefore, for a realistic analysis elasto-plastic behaviour 
IS suggested to be considered 

• Strain hardening ( 20-80 percent) has no significant effect on the strain 
energy release rate 

• Growth of plastic zone is sensitive to the displacement boundary conditions 
The plastic zone is initiated at the free-edges and gradually spreads over 
the core of the adhesive Growth of plastic zone is also sensitive to the 
crack propagation (size of the crack) 

• Strain energy release rate is sensitive to layup sequences in composite 
adherends A proper choice of layup sequence in the composite adherend 
results in an efficient bonded joint 

• Use of adhesive of high fracture toughness and large overlap lengths 
results m stronger bonded joint, whereas, employment of small thickness 
of adhesive results in stronger bonded joint 
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Table 3 1 Specifications of double lap joint 


Adherend (aluminium) 

Young’s modulus (£>) 70 GPa 

Poisson’s ratio (ur) 0 33 

thickness (tr) 10 mm 


Adhesive (Epon VIII) 


Young’s modulus (E,) 3 5 GPa 

Poisson’s ratio (i/,) 0 41 

thickness (t,) 2 mm 

yield stress (cTo) 515 MPa 


Table 3 2 Specifications of double lap joint 


Adherend (CFRP) 


longitudinal modulus (Ei) 137 GPa 

transverse modulus (Et) 7 GPa 

Poisson’s ratio (i/ft) 0 3 

thickness (tr) 10 mm 


Adhesive (Epon VIIl) 

Young’s modulus (Es) 3 5 GPa 

Poisson’s ratio (i/J) q 41 

thickness (i,) 2 mm 

yield stress ((T,) 51 5 MPa 

fracture toughness (Gc) 0 36 kJ/m^ 
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Table 3 3 Specifications of double lap joint 
Double lap joint 


overlap length (Lo) 40 mm 

total length (!»<) 80 mm 

Adherend (carbon fibre reinforced plastic) 

longitudinal modulus (Ei) 137 GPa 

transverse modulus (Et) 7 GPa 

shear modulus Gu 4 5 GPa 

Poisson’s ratio (vu) 0 3 

adherend thickness (i^) 10 mm 

Adhesive (Epon VIII) 

Young’s modulus (Eg) 3 5 GPa 

Poisson’s ratio (y*) 0 41 

thickness (Q 2 mm 

fracture toughness (Gc) 0 36 kJ/m^ 


fracture cntenon Gj + Gjj > Gc 
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Table 3 4 Computed failure loads (P/’s) for various g (= ) values and 

different laminae combinations 



ttcr (mm) 

Failure loads in N/mm for Lo/Lt= 050 & 0 20^ 

[070V070“]s [07457907-45“]s I070790790“k 

0 80 

0117 

22 000 

22 750 

21 250 

0 85 

0117 

23 125 

23 500 

22 375 

0 90 

0117 

26 125 

24 250 

22 750 

0 95 

0117 

25 000 

25 000 

23 125 

100 

0117 

25 750 

26 125 

24125 

105 

0117 

26 500 

26 500 

24 625 

1 10 

0117 

27 250 

27 250 

25 000 

1 15 

0117 

28000 

29500 

26 875 

120 

0117 

28 750 

29 500 

26 500 


Table 3 5 Computed failure loads (P/’s) for various 1 (« 
different laminae combinations 


values and 




^cr 


, Failure loads in N/mm for Ge/Gc= 1 00 & tJt 
(mm) [0707070°]^ [07457907 -45°].9 fO 


0 20 


/0°/90°/90°lg 


0 250 
0277 
0 312 
0357 
0416 
0500 


0 087 
0 091 
0 095 
0100 
0107 
0117 


16 375 
15 250 
15 250 
17 125 
19000 
25 750 


16 375 

14 875 

15 625 

16 937 
19375 
26125 


16375 
15 250 
15 625 
17500 
19000 
24 250 
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Table 3 6 Computed failure loads (P/’s) for various t (= values and 
different laminae combinations 


ta/tr 

ttcr (ram) 

Failure loads m N/mm for (?'e/Gc= 
[OVOVOVO'-Js [O 7457907 - 45°]s 

100 050 

[070790790“]s 

Blllig 

0117 

47 50 

41 500 

5650 


0117 

3175 

29 500 

3100 


0117 

29 50 

30250 

26 50 


0117 

28 75 

28 000 

25 75 


0117 

25 75 

26125 

2425 














35 



Crack Length (mm) 


Figure 3 2 Schematic of convergence of failure loads for various tolerance 
values 



Figure 3.3* Schematic of hardemng parameter k 
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Figure 3 4 Companson of deflections at location ’C’ of cantilever beam 
(Zienkiewicz, et al , 1969) 
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(a) Type - 1 Boundary Condition 



(b) Type - 2 Boundary Condition 



(c) Type - 3 Boundary Condition 


Figure 3 5: Boundary conditions and the corresponding growth of plastic zones 
in the adhesive matenal 
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Figure 3 6 Schematic of adhesively bonded double lap joint and the finite 
element mesh showing the crack location ’3’ 
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Figure 3 7 Growth of plastic zone with crack propagation 



Crack Length (mm) 


Figure 3 8- Variation of failure load with tomates 

location ’3’ and different g (= Go/C?c)values m [0 /O CFRP lamina 
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Figure 3 9 Variation of failure load with crack length for crack initiation at 
location ’3’ and different g (= G'^/Gc) values in [0°/45“/90'’/ - 45®ls (optimal) 
CFRP laminates 
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Figure 3 10. Variation of failure load with crack length for xrack initiation 
at location and different g (= G^/Gc) values in [O°/O°/^/90^]s CFRP 
laminates 
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Figure 3 11* Vanation of failure load with crack length for crack initiation at 
location 3 and different 1 (= values in [0°/0°/0°/0°]s CFRP laminates 




350 



Figure 3 12 Vanation of failure load with crack length for crack imtiafton at 
location 3 and different 1 (= I./i,)values in [0“/45V90V - 45°]s (optunal) 
CFRP laminates 
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Figure 3 13* Vanation of failure load with crack length for crack initiation at 
location 3 and different 1 (= values in [070790V90°ls CFRP laminates 



Figure 3 14 Vanation of failure load with crack length for crack initiation at 
location ’3’ and different t (= values in [QP /0° /QP /Q^s CFRP laminates 




^5 



Crack Length (mm) 


Figure 3 IS- Vamtion of failure load with crack length for crack inihahon at 
loition ’3’ and different t (= f./*.)value» in [ff/dS'/SO”/ - dS-]* (optimal) 


CFRP laminates 
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Figure 3 16 Vanation of failure load with crack length for crack initiation at 
location ’3’ and different t (- t./tr) values in [0“/0®/90“/90°]s CFRP laminates 




Chapter 4 


Visco-plastic analysis 


4.1 Introduction 


Adhesives are mostly polymenc materials and their responses are time dependent 
Groth (1986) concluded that events like moisture diffusion, delayed failure are 
necessary and should be taken into account for an accurate analysis of adhesively 
bonded joints At high stress levels, the non-linear visco-plastic effects could 
include creep strains These would lead to significantly different results from 
the linear predictions and therefore, these effects are important to be mcluded 
With the increase in load level, generally, adhesive layers become visco-plastic 
earlier than the adherends Therefore, in this analysis, adhesive is considered as 
visco-plastic in behaviour, while adherends are considered to be linear elastic 

Finite element code, with elasto-visco-plashc behaviour, is developed and 
validated with the available results from the hterature Visco-plastic analyses 
of single lap, stepped lap, and double lap joints are earned out Influences 
of vanous parameters, such as (i) stactang sequence in composite adherends, 
( 11 ) critical strain energy release rate of the adhesive, (ni) overlap length and 
(IV) thickness of the adhesive on the strength of the adhesively bonded joints. 
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are also investigated 


4.2 Formulation 


Biot (1954) has derived equations to represent visco-elastic and relaxation 
responses by employing potential and dissipation functions A large number of 
phenomena involving interaction of diffusion, chemical reaction, heat transfer, 
mechamcal deformation etc are included in this theory Zienkiewicz and 
Cormeau (1974) developed an elasto-visco-plastic algorithm, which is found to 
be an efficient tool and is applicable to a wide range of non-linear (matenal) 
problems The solution technique developed for visco-elastic behaviour and 
could be coupled with visco-plastic phenomenon with series However, most 
efficient way of dealing with large deformation and dynamic effects is not 
resolved Zienkiewicz, ct al (1974) discussed the finite element method of 
solution of elasto-visco-plastic engineering problems For the computational 
convenience, bulk material is considered to be composed of several layers, where 
each layer might have different material charactenstics This model was limited 
to isotropic and includes few material laws Kanchi, et ol (1978) discussed 
a visco-plastic technique and employed this to study the problems of plastic 
and creep mstabihty Schapeiy (1984) discussed pertinent aspects of J integral 
and energy release rate theory for non-lmear elastic media He developed 
methods of quasi-static deformation and fracture analysis for non-linear visco- 
elastic media The corresponding principles, which provide the basis for the 
analysis are applicable to crack growth, crack closing and other problems such 
as ablation and interfacial contact and separation 
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I Strain rate 

Groth (1990) presented expressions for visco-elastic , visco-plastic and visco- 
elasto-plastic strain rate as written in Equations (4 1), (4 2) and (4 4), respectively, 
for the one dimensional models of adhesive raatenal under pure shear as shown 
in Figure 4 1 


€vp = 'Yxylr^y-ro] \fT^y>rQ (4 2) 

= 0 if < To (4 3) 

€vep = 7xy[Ta:y-Gzy€ve] + 7xy[Tzy-To] if > Tfl (4 4) 

” *V*y \y^xy if Tiy ^ Tq (4 5) 


From the above Equations (4 1), (4 2) and (4 4), it is observed that, the 
expressions for strain rates are different for different adhesive models Visco- 
plastic strain rate €vp for two dimensional problems is expressed as 


Cyp = 7{^)^ 


(4 6) 


where, 7, # and ^ are fluidity parameter, visco-plastic flow function and 
flow vector, respectively (^) indicates for positive values of ^ Fluidity 
parameter is a material constant Different choices are recommended (Owen 
and Hinton, 1980 and Zienkiewicz and Taylor, 1991) for the function # and 
three most common representations are 
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# = 


(4 7) 


\ CTq / 

<p = 

r-<-i 

(4 8) 


\ CTq / 

^ = 

F 

(4 9) 


< 7-0 


where, M and N are material constants and <ro is yield stress Flow rule 
jP, is assumed as 


^ + + (410) 

Flow vector (&) is given as 



(411) 


Deviatoric stress components are denoted by primes and these are obtained 
by 



(4 12) 


4 represents the second invariant of deviatoric stresses and written as 
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•4 


-(<ry^ + (TyCy) + T^y 


(4n) 


I Visco-pIastic consbtuUve relation 


The visco-pIastic stress-strain relation is written as 


= [d-‘ + C"]‘‘ 


(414) 


D IS the material matrix of linear elastic matenal and C" is given as 


(7‘ = 


(415) 


where, 6> is a constant. At,, is the incremental time of n‘* interval and is 
given as 

=: piTTii + pzmz 


where, pi, pz, mi and mz are expressed as 


P2 = 


I Z M \ 


(417) 

(418) 


mi 


2/3 -1/3 0 

-1/3 2/3 0 

0 0 2 


( 419 ) 


7 
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/ / 

/ / 
<=^x<^V 


/ / 



2o^a;'7a;y 




(4 20) 


4.2.1 Finite element technique 

A computer code, similar to one developed by Owen and Hinton (1980), is 
wntten to predict the failure load, which is explained m the Figure 4 2 The 
module, which takes into account the non-linear visco-plastic behaviour of the 
matenal is indicated by the dashed portion in the flow chart Detail of the 
elastic and visco-plastic analyses are explained in the mam algorithm of Chapter 
2 and the following visco-plastic algorithm, respectively Visco-plastic algorithm 
IS explained as follows 


I Non-linear (visco-plashc) algorithm 

step 01 at t = tn find the following 
B- = Bo + Br,L(S-) 

+ C"]-i 

= J[B^fD!^B^da 
^ = t(^)^ 

step 02 compute incremental displacements and stresses 

(if i = to calc + A^) 
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Af>- = 

2l(r” = Ap(B’‘4\6” - e"p41t") 


step 03 

update displacements and stresses 

6”+’ = 6" + Co2l6” 

^n+l _ ^ 

step 04 

calculate visco-plastic strain 

step 05 

calculate proportionality factor H 

OrH+l _ 

Jtl = 

CTO — <r” 

step 06 

compute residual stress to be compensated 

(Tr = 

step 07 

update stress 

= (1 - 

step 08 • 

apply equilibrium correction 


/<? 

where 

^n+l ^ /[j3«+i]r^j»+idl2 + /«+! 

/? 

where 

j^n+i represents nodal forces 

step 09 • check if visco-plastic strain rate is close to steady state value 
< Css 1 e tolerance 

if steady state is reached go to step 07 of mom olgonthm in 

Section 2 2 

else go to step 01 

S”, or«, and represent strain displacement matrix, stress tensor, 

visco-plastic stress-strain relation and stiflE^ness matrix, respectively of n^^ iteration 
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B", <r”, C", and if" are functions of the displacement vector of n"* iteration 
One is supposed to update these values in every iteration But there is 
a possibility of divergence in the results by using the iiT” values Using 
constant values of B” and and Newton Raphson’s modified approximation 
(Owen and Hinton, 1980 and Zienkiewicz and Taylor, 1991) one can achieve 
convergence Further, this reduces substantially the computational time and 
memory However, and are updated in each time step 

If the effective stress of any element in the adhesive, exceeds the yield 
value of the material <to, then the material state becomes visco-pIastic and 
visco-plastic analysis is earned out Visco-plastic module of elasto-visco-plastic 
analysis is indicated by dashed box of Figure 4 2 

Visco-plastic calculations are repeated until the strain rate has reached 
a steady state value Then stram energy release rates are calculated and, as 
mentioned earlier in Chapter 2 The fracture criterion Gi Gji > Gic 4* Guc 
IS employed to predict the failure loads 

The developed code is quite general, and could be used for visco-elastic, 
visco-plastic and elasto-visco-plastic models of adhesive materials One needs 
to incorporate the corresponding behaviour of the models for adhesives m the 
expression for stram rate 


4.3 Validation of computer code 

4.3 1 Perforated plate 

In order to validate the developed computer code the two dimensional visco- 
plastic perforated plate problem considered by Zienkiewicz, et al (1974) is 
considered Schematic of this plate is shown in Figure 43 The physical 
dimensions and the material properties of the plate are listed in Tkble 4 1 The 
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stress level is plotted in Figure 4 4 against the time for an element P, which is 
in plastic range The result is in fairly good agreement with that by Zienkiewicz, 
et al (1974) 


Further, to check the accuracy of the computer code developed, a single lap 
joint and a stepped lap joint are investigated Stress and strain distributons 
along the overlap lengths of these joints are compared with those available in 
open literature 


4.3.2 Single lap joint 

To check the accuracy of computer code of visco-plastic finite element analysis 
a single lap joint (Figure 45) is employed Details of the joint are listed in 
Table 4 2 Stresses, along the overlap length of single lap joint (Hiregoudar, 
1993), are computed from elastic and visco-plastic analyses and these results 
are compared in Figure 4 6 The maximum difference in the computed peel 
stress along the overlap length of the single lap joint is less than 3 6 percent 
only While, this difference for shear stress is less than 11 4 percent Stresses 
predicted by the present visco-plastic analysis and visco-plastic analysis done 
by Hiregoudar (1993) are also compared in Figure 4 7 The difference in the 
computed peel stress along the overlap length of the single lap joint is of the 
order of 11 percent It is also observed that this difference is maximum near 
the edges and gradually reduces to a minimum towards the central part of the 
overlap Because of the proper refinement in the finite element mesh near the 
edges, stress distnbutions are expected to be more accurate 
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4 3 3 Stepped lap joint 

Further, validation of the computer code of finite element developed is carried 
out on a stepped lap joint (Figure 4 8) Joint specifications are listed in 
Table 43 In Figures 4 9 the shear stress and shear strain distributions, m the 
adhesive along the overlap and near the interface of adherend and adhesive^ are 
compared in Figure 49 with those of Su and Mackie (1993) The difference is 
found to be of the order of 11 7 percent It is also observed that this difference 
IS tnaYimiim near the edges and gradually reduces to a minimum towards the 
central part of the overlap length Because of the more refined mesh near the 
edges as compared to Su and Mackie (1993), present results are expected to 
be accurate 


4.4 Parametric study 


Effects of crack location, stacking sequence, critical strain energy release rate 
of adhesive material, overlap length and adhesive thickness on the strength of 
the double lap joint (Figure 4 10) are investigated Because of symmetry only 
one half of the joint is considered for the analysis The problem is treated as 
one of plane strain model and a finite element mesh consisting of 1816 CST 
elements is employed for the analysis 

Three non-dimensional parameters are defined viz 


g = 


Gc 

Lo 

Lt 


1 


(4 21) 
(422) 
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t 



(4 23) 


where, G^c and Gc are the fracture toughnesses of the adhesive used and 
the standard adhesive (epoxy^), respectively Lo and Lt denote the overlap length 
and the total length of the bonded joint In addition t, and tr represent the 
thicknesses of adhesive and adherend, respectively is considered to be in 
multiples of Gc Assuming constant fracture toughness of the adhesive (Gc), 
adherend thickness tr and overlap length Lo parametric study is earned out 
The adherends and adhesive of the joint, are considered to be carbon fibre 
reinforced plastics (CFRP) and Epon VIII (Table 44), respectively Cntical 
crack location ’3’ is considered for further discussion m this Chapter 


4.4.1 Stacking sequence 

Symmetric quasi-isotropic combinations are considered to mvestigate the effect 
of stacking sequence on the bond strength Four fibre orientations are 
considered They are 0% 45°, 90° and -45° Among all the possible symmetnc 
combinations, it was observed in Chapter 2 (elastic analysis) that [0°/45°/90 / — 
45°]5 has the maximum bond strength among its combmations Also the presence 
of 0°, ply in the immediate neighbourhood of the adhesive results m stronger 
joint Visco-pIastic analysis is earned out for [0°/0°/0°/0°]s, [0°/45°/90°/ - 
45°]s and [0° /0° /90° /90°]s layups of CFRP adherends It is found that employing 
[0° /0° /90° /90°]s layup results in a stronger joint is obtamed than [0°/0 /O /O ]s 
layup This combination is found to be a stronger joint than [0“/45°/90°/-45°]5 
layup The failures loads are computed for vanous adhesive fracture toughnesses 
ratio g (= dc/Gc), overlap lengths I (- Wii)and adhesive thickness ratio 
t (= t,/tr)foT these three layups Detail of these results are discussed in the 

respective sections 
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4 4.2 Fracture toughness of adhesive 


For various g (= C'c/Cc) values and layups the computed failure loads arc 
listed m Table 4 5 It is observed that as g increases failure load increases for 
all the three layup sequences considered Adhesives with higher critical strain 
energy release rate ( high g value) lead to a stronger joint Therefore, in order 
to achieve a higher fracture strength of the joint, adhesives with higher critical 
strain energy release rate are suggested Variation of failure load with crack 
length for vanous g values in [0®/0“/0“/0®]s’, [0®/45‘’/90“/ — 45"|5 (optimal) 
and [0°/0°/90°/90°]s CFRP laminates are listed in Figures 4 11, 4 13 and 4 12, 
respectively 

Though fracture toughness of adhesive is generally denoted by Kicy the 
stress intensity factor for ideal plane strain adhesive material, in the present 
study, for convenience cntical strain energy release rate Gc is also termed as 
fracture toughness 


4.4J Overlap length 


For various values of I (= Z-o/I-t) and the three layups chosen the computed 
failure loads are listed in Table 4 6 It is observed that as 1 increases failure 
load increases for all the three layup sequences For higher fracture strength 
of the bond larger overlap length Variation of failure 

load with crack length for various I values m^f0®/45°/90°/ - 45®]5 (optimal) 

and [070790°/90°]5 CFRP laminates are shown in Figures 4 14, 4 16 and 4 15, 
respectively 
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4 4 4 Adhesive thickness 

For various t (= values and layups the computed failure loads are listed m 
Tabic 4 7 It is observed that as t increases failure load decreases for the three 
layup sequences For higher fracture strength of the bond small thickness of 
the suggested Variation of failure load with crack length for 

various t values inJ0745°/90®/ - 45°]s (optimal) and [0‘’/0°/90“/90°]5 CFRP 
laminates are listed in Figures 4 17, 4 19 and 4 18, respectively 


4.5 Conclusions 


On the basis of elasto-visco-plastic analysis earned out in this chapter, the 
following conclusions are drawn 


• The clastic analysis predicts lower failure loads as compared to the elasto- 
visc 0 “plastic analysis Therefore, for a realistic analysis elasto-visco-plastic 
behaviour is suggested to be considered 

• Strain energy release rate is highly influenced by the adhesive properties 
such as fluidity parameter and flow rule 

• Stress distnbutions, along the overlap length of the single lap joint, are 
predicted by the present elasto-visco-plastic analysis and the results are 
found to agree with the those reported by Hiregoudar (1993) 

• Stress distributions along the overlap length and across the adhesive 
thickness of the stepped lap joint, are predicted by the elasto-visco-plastic 
and the results agree with the those reported for thick adherend shear 
test (TAST) specimen by Su and Mackie (1993) 
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• Strain energy release rate is sensitive to layup sequences in composite 
adherends A proper choice of layup sequence m the composite adherend 
results in an efiScient bonded joint 

• Use of adhesive of high fracture toughness and large overlap lengths 
results in stronger bonded joint, whereas, employment of small thickness 
of adhesive results in stronger bonded joint 
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Table 4 1 

Specifications of the perforated plate 

problem (Zienkiewicz, et al , 

1974) 

Perforated plate 



dimension of plate 

36 mm X 20 mm 


diameter of the circular cut out 

10 mm 


Young’s modulus (E) 

70 Gpa 


Poisson’s ration (i/) 

02 


stress applied (a) 

95 MPa 


yield stress (ao) 

243 Mpa 


fluidity parameter (7) 

001 Sec-i 



162 


CHAPTER 4 VISCO-PLASTIC ANALYSIS 


Table 4 2 Specifications of single lap joint (Hiregouciar, 1993) 


Single lap joint 


overlap length (L^) 16 ram 

total length (Lt) 80 mm 

applied stress (a) 200 MPa 

adhesive thickness (t,) 0 30 mm 

adherend thickness 1 60 mm 


Adhesive (FM 73) 


Young’s modulus (Eg) 
Poisson’s ratio (i/s) 
yield stress (cr^ ) 
fluidity parameter (7 ) 
constant N 
constant A 


2 21 Gpa 
043 
40 Mpa 
4495 10 
1426 
14 


flow function ( 0 ) 




Adherend (Steel) 


211 3 GPa 
033 
00 


Young’s modulus (Er) 
Poisson’s ratio (i/^) 
yield stress (o-o ) 
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Table 4 3. Specifications of stepped lap joint (thick adherend shear test 
specimen) of Su and Mackie (1993) 


Joint configuration 


overlap length (!<,) 8 mm 

total length (Z-t) 125 mm 

applied load (P) 3 5 kN 

total joint thickness (t) 13 35 mm 

width of the joint (W) 25 4 mm 


Adhesive (rubber-toughened epoigr) 

Young’s modulus (Es) 6 60 Gpa 

Poisson’s ratio (ug) 0 38 

thickness of adhesive (*,) 0 65 mm 


Adherend (Steel) 

Young’s modulus (Er) 210 Gpa 

Poisson’s ratio (vr) 0 3 

thickness of adherend (fr) 6 35 mm 
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Table 4 4 Specifications of double lap joint 


Double lap joint 


overlap length (Lo) 40 mm 

total length (Lt) 80 mm 

Adherend (carbon fibre reinforced plastic) 

longitudinal modulus (Ei) 137 GPa 

transverse modulus (Et) 7 GPa 

shear modulus Gu 4 5 GPa 

Poisson’s ratio (i/n) 0 3 

adherend thickness (^r) 10 mm 

Adhesive (Epon VIII) 

Young’s modulus (E^) 3 5 GPa 

Poisson’s ratio (y,) 0 41 

thickness (t,) 2 mm 

fracture toughness (Gc) 0 36 kJ/m^ 


fracture cntenon G/ + G// > Gc 
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Table 4 5 Computed failure loads (P/’s) for various g (= G^/Gc) values and 
different laminae combinations 


GdGc 

a„ (mm) 

Failure loads 
(0 / 070 V 0 ° 1 . 

in N/mm for LofLi- 
? [070V90V90“ls 

= 050 & 0 20 

[O 7457907 - 45°l5 

01 

0164 

5 760 

5 984 

5 696 

02 

0164 

8 096 

8416 

8 032 

0"^ 

0164 

9 920 

10304 

9 856 

04 

0164 

11520 

11904 

11392 

05 

0164 

12 800 

13 312 

12 672 

06 

0164 

14 080 

14 592 

14 016 

07 

0164 

15 232 

15 872 

15 168 

08 

0164 

16 384 

17 024 

25 472 

09 

0164 

17472 

18304 

16 192 

10 

0164 

18560 

19 456 

17408 

1 1 

0164 

19 840 

20 736 

18 944 

12 

0164 

20992 

22016 

20 096 

13 

0164 

22400 

23 296 

21 632 

14 

0164 

23 552 

24576 

23 040 

15 

0,164 

24704 

25 984 

24 320 


Table 4.6: Computed failure loads (P/’s) for vanous 
different laminae combinations 

1 (= LoZ-^t) values and 

lotU 

Ocr (mm) 

Failure loads in N/ram for Gfc/Gc 
[0707070“}5 {070790790"]g 

= 040 & 0 20 

[07457907-45°]5 

0307 



— 

— 

0333 



4607 

4 271 

0363 

0059 

5112 

4 968 

4 607 

0400 

0065 

5 412 

5 448 

5 040 

0444 

0 073 

5 832 

5 880 

5 520 

0500 

0164 

11616 

12048 

11472 


plastic flow 
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Table 4 7. Computed failure loads (P/’s) for various t (= 4, /L) values and 
different laminae combinations 



ttcr (mm) 

Failure loads in N/mm for C^c/^c 
[07070V0°]s (070790790°js 

= 040 & LJLt= 050 

[07457907 - 450 ].,; 

0012 

0164 

16 384 

17 792 

!Fi28 ^ ■ 

0 025 

0164 

15 488 

16 000 

14 848 

0 033 

0164 

15 040 

15 296 

14 528 

0 050 

0164 

14336 

14 592 

13 696 

0100 

0164 

13 056 

13 440 

12 672 

0 200 

0164 

11616 

12 048 

11472 



Figure 41 Adhesive material models (a) the three 
(b) the Bingham visco-plastic model, and (c) the five 
viscoplastic model '■ ^ 


parameter solid model, 
parameter viscoelastic- 
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Figure 4 2: Flow chart of visco-plastic finite element analysis 
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Figure 43 Schematic of perforated plate under uniform tension (Zienkiewicz, 
et al , 1974) showing the plastic element V and the finite element mesh of the 
plate 
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Tina (Saconds) 


Figure 4 4 Vanation of stress with tune for plastic element V of the perforated 
plate (Figure 4 



Figure 4 5» Schematic of single lap joint 


considered by Hiregoudar (1993) 
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Overlap Length 


Figure 4 6 Comparison of elastic and visco-plastic stress distribution along the 
overlap length of single lap joint (Hiregoudar, 1993) 
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Figure 4.7: Comparison of peel stress distribution along the overlap length of 
single kp joint (Hircgoudar, 1993) 



Figure 4.8: Schematic of the stepped lap joint (thick adherend shear test TAST 
specimen considered by Su and Mackie, 1993) 
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Figure 4 9 Companson of shear stress and strain distribution along the overlap 
length at interface of adherend and adhesive of the stepped lap joint 
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Figure 4 10. Schematic of adhesively bonded double lap joint and the finite 
element mesh showing the crack location ’3’ 
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Figure 4 11 Variation of failure load with crack length for crack initiation at 
location 3’ and different g (= values in [0707070“ I5 CFRP laminates 
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laminates 


50 

initiation 
’]s CFRP 
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Figure 4 13 Variation of failure load with crack length for crack initiation at 
location ’3’ and different g (= Crc’/^o) values in [0®/45°/90“/ ■“ 45°]s (optimal) 
CFRP laminates 



177 



Crack length (mm) 

Figure 4 14 Variation of failure load with crack length for crack initiation at 
location 3’ and different I (= Xo/i*) values in [0°/0°/0°/0°]s CFRP laminates 



Crack length (mm) 

Figure 4 IS. Vanation of failure load with crack length for crack imtaton at 
loLon ’3- and different I (= Vi.)valucs m [0“/0“/9(F/90»lr CTRP laminates 
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Figure 4 16 Vanation of failure load with crack length for crack initiation at 
location ’3’ and different 1 (= values in I0“/457907 ~ 45‘’]5 (optimal) 

CFRP laminates 
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Figure 4 17. Variation of failure load with crack length for crack initiation at 
location ’3’ and different t (= ig/tr) values in [0°/0®/0°/0°]5 CFRP laminates 



180 



Crack length (miu) 


IoS^n‘‘-3*an’?d?ff? "“h «ack length for crack inittahon at 

tion 3 and different t (= values m [0»/0»/90“/90»Is CFRP laminates 
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Figure 4.19: Variation of failure load with crack length for crack initiation at 
location ’3’ and different t (= f,/tr)values in [0745790°/ - 45°]s (optimal) 
CFRP laminates 



Chapter 5 


Experiments 


5.1 Introduction 


The fracture cntena employed m structural analyses are based on pure modes 
I and II and various combinations of them Davies, et ol (1992) discussed 
interlaminar fracture testing of carbon fibre reinforced epoxy and PEEK 
composites They found large amount of scatter m the results of strain 
energy release rates in modes I and II, respectively Therefore, these fracture 
parameters are required to be determined experimentally Also, accuracy in 

X. 

the testing would reflect in the estimation of the strength of a bonded joint 


In this Chapter expenmental determination of critical strain energy release 
rates in modes I and II, and the necessary theory is explained Also, the 
strengths of the double lap bonded joints are determined experimentally and 
compared with those obtained from the finite element analysis considering 
various fracture criteria With the limited observations, the most critical fracture 
criterion is suggested 
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5.2 Strain energy release rate in mode I 


The energy approach discussed by Carlsson and Pipes (1987) is based on a 
thermodynamic criterion for fracture by considenng the energy available for 
crack growth of the system, and the surface energy required to extend an 
existing crack Energy balance of a cracked body can be represented as 


H ^ tJ-'W ( 51 ) 

where U and W represents the potential ener®i, the elastic strain 
energy stored in the body and the work done by the movement of external 
forces, respectively The strain energy release rate in defined as 


G, = -g (5 2) 

where, dA is the newly created cracked area Figure 5 1 shows a typical 
load versus displacement response for a cracked body Assuming ^rack growth 

IS under constant load (fixed load) P and the vertical displacement is u 

f fir. 


I For contant load 


dW * PAv 
dtf = 


(5 3) 


area area ABC 


( 54 ) 
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Av)-^Pv 

(5 5) 


(5 6) 

-\pAv 

(5 7) 


For linear elastic body, the reactions between the load and displacement 
may be expressed as 


V = CP 


(5 8 ) 


where, C is the compliance of the specimen Using Equations (5 7) and 
(5 8) m Equation (5 2), the strain energy release rate under constant load is 
written as 


Gi = 


2 dA 


(5 9) 


I For contant displacement 


Similarly, for crack groth under constant displacement condition (fixed grip) 
work done W" vanishes and rraer 


(5 10) 


dff = area AFC - area ABC 
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= ^(P-AP)v-ipv 


1 

r^p 


(511) 

(512) 


and 


Gj 


1 dP 
2^dA 


(513) 


Substituting Equation (5 8) in Equations (5 9) and (5 13), respectively one 
can obtain 


= 

dc 

20^0 A 

(514) 

Gi = 

P^dC 

2 dA 

(515) 


From Equations (5 9) and (5 15), crack growths under constant load and 
constant displacement yields the same expression for mode I strain energy release 
rate Therefore, depending on the convenience one can conduct the experiment 
either under constant load or constant displacement Several methods have 
been adapted to experimentally determine the mode I cntical strain energy 
release rate These include 


• Compliance method using DCB specimen 


• area method 
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• modified area method 

• method using theoretical formulae 

• method using non-linear beam theory 

Compliance method using DCB specimen is widely used in the case of 
metals and also in composites because of the following features 

- crack growth is stable under constant displacement loading (faxed grip) 

- error introduced in analysing data is small 

- each DCB specimen will provide a Gjc value 


Wilkins, et al (1982) discussed the characterisation of dclamination growth 
m graphite epoxy Further, descnption on various factors that might affect the 
load displacement response of slender DCB specimen under mode I testing of 
compliance method is given by Keary, et ol (1985) Following precautions are 
to be taken while conducting the experiment 

- the crack remains m a plane 

- crack-tip plane is perpendicular to the length of the DCB specimen 

- crack plane is perpendicular to the load line 

- the cantilevers of the specimens are thin plates and the contribution 
of shear stress developed due to bending moment is negligible in 
companson to bending loads 

- large deflections in beam bending, which could generate geometrical 
non-lmearity is avoided 


From elastic beam theory the displacement can be written as 
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V 


2Pa^ 

~'^EJ 


( 516 ) 


where, E, I and a represent the Young’s modulus, moment of inertia of 
cross section and the ciack length of the cantilever, respectively Using the 
compliance as expressed m Equation (5 8) and substituting for v as expressed 
in Equation (5 16) 


^ P 

(517) 

r - 2 “^ 

^ 3BI 

(518) 

C * Aic? 

(519) 

^ 3EJ 

(5 20) 

u4i IS a constant for a given specimen 

Also, the above equation could be 

rewritten in logarithm form as 



InC ~ In Ai + 3 In a 


Gi 


P^dC 
2 dA 


(5 22) 
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dC 
2W da 

EIW 


(5 23) 
(5 24) 


where, W represents the width of the DCB specimen Since, critical load 
Pc corresponds to critical strain energy release rate Gic, above equation is 
rewritten as 


EIW 

(5 25) 

Pc = \/GicEIWa-^ 

(5 26) 

1 

II 

(5 27) 

^2 = \/GicEIW 

(5 28) 


Ai is a constant for a given specimen Also Equation (5 27) could be 
rewritten in logarithm form as 


In Pc = 


(5 29) 


In Az— In a 
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Eliminating El from Equations (5 20) and (5 28), the critical strain energy 
release rate in mode I G;<- is expressed as 


2W (5 30) 

The unidirectional composite laminates are prepared by hand layup techniq- 
ue on flat steel plates with carbon fibre (G-808) with a warp to weft ratio 
of 91 9 and area density 220 gm/m^ (manufacturer Brochier SA, France) and 
the resin Resin is prepared by mixing Araldite LY 556 and Hardener HY 
951 (manufacturer Ciba-Geigy Ltd , India) with a weight ratio 10 1 These 
composite laminates are cut to 30 mm width strips Each pair of strip is glued 
with epoxy resm with bi-axially oriented polypropylene (BOPP) filler of thickness 
20 fxm IS embedded at one end This BOPP membrane provides openmg space 
to attach the fixture with hinge facilities for mode I loading condition and a 
sharp through width crack After cunng BOPP stnp is removed, holes are 
made and the fixture is fitted with screws To measure the crack length at each 
load increment, on the thickness sides thin strips of graph paper are glued 
Specimen is loaded m mode I loading as shown in Figure 5 2 Experunent 
IS earned out in displacement control condition (fixed grip) with a cross head 
speed of 1 mm per minute Initially it was loaded till the crack just begins to 
grow' and unloaded till the load level comes down to zero value Loading and 
unloading cycles were repeated for 5 to 8 times till the crack grows to half of the 
length of the specimen At the end of each loading cycle, the crack sizes were 
noted from either sides of the DCB specimen with the help of a magnifying 
lens and a torch Average value of these crack lengths were considered for the 
calculation of Gjc 

Load versus displacement curves for all 8 specimens are obtained A typical 
load versus displacement curves for loading and unloadmg cycles is shown m 
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Figure 5 3 For crack length corresponding compliance C are calculated from 
loading cycle of the graph and the corresponding critical loads Pc arc noted 
From Equation (5 30) it is clear that knowing Ax and Az one can hnd G/c for 
a DCB specimen The values of Ai and Az can be found out graphically from 
the logarithm plots of Pc vs a and G vs a, respectively Also, these values could 
be determined by linear regression of data, which is given as 


1 " 

InAi — — G, + 3 Zn a,} 
1 " 

InAz = — Y]{ln Pcx — In a,} 
n ^ 


(5 31) 
(5 32) 


where, n represents the number of load cycles Thus, using the above 
linear regression method, critical strain energy release rate for mode I (Gw) is 
calculated for the SEN specimens and listed m Table 5 1 The average value 
of Gw (0 446 kJ/m?') is considered to predict the strength of bonded joints 

These Gw values are not close a single value The reasons could be 
non-umform thicknesses of the adherends and adhesives, defects in adherend 
and adhesive, uneven surface preparation, improper gluing etc 


5.3 Strain energy release rate in mode II 


For the determination of critical strain energy release rate in mode II loading, 
standard experimental procedure using ENF specimen and the necessary theory 
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IS explained by Carlsson, el el (1986) In the present work, an ENF specimen 
IS employed for the experimental determination of critical strain energy release 
rate in mode II condition This is a three point bending flexure specimen 
(Figure 5 4) with an embedded through width crack at the interface of two 
laminates The specimen has been found to produce shear loading at the 
crack-tip without much friction between the surfaces (fretting) From the beam 
theory an expression for the strain energy release rate in mode II is explained 
by Carlsson, et al (1986) The expression is written as 


Guc 


9P^Ca^ 

2W(2L'^ + 3a?) 


(5 33) 


where, C, P, a, W, and L represents the the compliance, the applied load, 
the crack length, the width and the semi-span length of the specimen, respectively 
(Figure 5 4) This expression does not include the action of shear surfaces and 
the surface of friction between crack surfaces (fretting) They also derived the 
condition for stable crack growth, which is expressed as 


a > OIL 


(5 34) 


But usually the crack length a is of the order of 0 5i leading to unstable 
crack growth A load displacement curve is obtained for each specimen The 
specimen under three point bending is shown in Figure 5 4 A cross head 
speed of 1 mm per minute was maintained during the loading Seven such 
specimens were tested and the critical strain energy release rates for mode II are 
determined and listed on Tkble 5 2 The average of these values (0 mkJ/rn?) 
IS considered for the prediction of the jomt strength 
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5.4 Strength of adhesively bonded joints 


Composite laminates are prepared as mentioned earlier These arc cut to 30 
mm, 20 mm and 15 mm width strips With proper surface preparation end 
tabs are glued Further surfaces are prepared to make the double lap joints 
(Figure 5 5) After curing tensile strength of these bonded joints are determined 
expenmentally and compared with those obtained from the elastic finite element 
analysis with various fracture cnteria Details of the failure criteria are listed 
m Chapter 2 Also a companson of failure load vs crack length curves based 
on these failure critena for the same specimen are illustrated 

The physical dimensions of the double lap bonded joints are listed m 
Table 5 3 Failure loads {Pf) of these joints from testing and finite element 
analysis with various fracture critena are presented in Table 5 4 In some 
joints failure occurred at the interface (adhesive), while m some it occurred m 
the adherend and for others the failure was a combination of both interface 
(adhesive) failure and adherend failure 'iQ'pes of failure occurred in the 
specimens are given in the Table 5 4 


5.5 Conclusions 


On the basis of experimental study and comparison of results of finite element 
analysis following conclusions are drawn 

• It IS observed that the experimental results are close to those results based 
on following fracture criterion 


(5 35) 


^7 "I" Gji > Gjc + Gjjc 
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where Gic and Guc are cntical strain energy release rates in modes I and 
II, respectively 

• The difference in critical strain energy release rate results, might be due 
to non-uniform thickness of the composite adherends and adhesive during 
specimen (joint) preparation 


Table 5 1 Experimental results for mode I critical strain energy release rate 
Gic 


Specimen 

No 

fracture toughness 
(kJ/rr?) 

average value 

1 

0726 


2 

0638 


3 

0642 


4 

0148 

0446 

5 

0246 


6 

0529 


7 

0233 


8 

0411 



Table 5 2 Experimental results for mode II cntical stram energy release rate 
Gijc 


Specimen 

No 

fracture toughness 
(kJ/m^) 

average value 
(kJ/m^) 

1 

0623 


2 

0778 


3 

0 666 

0664 

4 

0722 


5 

0587 


6 

0736 


7 

0536 
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Table 5 3 Various dimensions of double lap joint specimens (Figure 5 5) 


Specimen 

Number 

Lo 

(mm) 

U 

(mm) 

Width 

(mm) 

tr(Lower) 

(mm) 

(mm) 

tr(Upper) 

(mm) 

i 

103 

353 

19 13 

166 

009 

2 35 

2 

121 

339 

17 62 

126 

010 

175 

3 

102 

353 

17 65 

1 12 

010 

180 

4 

134 

325 

17 96 

1 19 

009 

173 

5 

121 

341 

17 66 

166 

010 

191 

6 

102 

217 

295 

172 

010 

166 

7 

101 

215 

290 

167 

009 

180 

8 

101 

216 

290 

163 

010 

176 

9 

98 

216 

285 

183 

010 

182 

10 

101 

217 

290 

185 

009 

182 

11 

97 

220 

285 

164 

010 

179 

12 

98 

220 

290 

167 

009 

168 

13 

98 

220 

290 

166 

010 

166 

14 

102 

220 

285 

171 

010 

169 

15 

90 

295 

18 9 

135 

010 

244 

16 

89 

296 

191 

133 

010 

2 53 

17 

100 

332 

13 8 

168 

010 

175 

18 

98 

336 

13 6 

170 

010 

165 
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Table 5 4 Comparison of experimental and finite element technique (elastic 
analysis) results for various fracture criteria 


experiment 


finite element analysis 


specimen 

number 

2 

4 

5 

6 

7 

8 

9 

10 
11 
12 
n 

14 

15 

16 

17 

18 


failure loads jP/s) 




criterion 
type 1 

entenon 
type 4 

criterion 
type 5 

entenon 
type 6 

entenon 
type 7 

failure 

type 

(10-"N) 

(10®N) 

(lO^N) 

(lO^N) 

(lO^N) 

(10®N) 

adhesive 

01569 

0 1250 

00791 

0 0978 

00875 

01043 

both 

0 2550 

01482 

0 0919 

01143 

01025 

01208 

both 

01467 

01501 

00958 

01153 

01050 

01237 

adhesive 

01414 

01506 

00958 

0 0895 

01066 

01249 

adhesive 

01159 

01155 

00856 

01040 

00927 

01084 

both 

0 2538 

01182 

00881 

01089 

00968 

01131 

both 

0 1397 

01370 

00868 

01069 

0 0939 

01111 

both 

01571 

01183 

00880 

01074 

0 0952 

01123 

both 

01675 

01232 

0 0784 

0 9564 

0 8483 

00999 

both 

0 2561 

01215 

0 0773 

00943 

00836 

0 0985 

both 

01369 

01369 

0 0874 

01063 

00942 

01111 

adherend 

01116 

01400 

00887 

01073 

00955 

01124 

adherend 

01402 

01408 

00893 

01080 

00961 

01131 

both 

01351 

01390 

00874 

01071 

00957 

01121 

adhesive 

0 2089 

01670 

01063 

01297 

01150 

01355 

adhesive 

01962 

0.1695 

01079 

01316 

01167 

01376 

adhesive 

0.1569 

01364 

00859 

01054 

00929 

01108 

adhesive 

01628 

01374 

00862 

01067 

00945 

01115 


whore, P, failure load, both both mtorface (adhesive) mdadherend 

entenon type 1 G, + Gn > G,c + Gnc , entenon type 4 

- criterion type 6 ^ ^ ^ 1 


criterion type 5 (7/4- Gn > Gjjc > 

entenon type 7 > 1 


\Cw/ \Oiic/ 
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Figure 5 1 Load-displacement behaviour for a cracked body under constant 
load and constant displacement 





Figure 5.2 Double cantilever beam (DCB) specimen under mode I loading 
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Figure 5.3: 
responses of 



schematic of load-displacement (crack opening displacement) 
double cantilever beam (DCB) specimen under mode I loading 



Figure 5 4: End notch flexure (ENF) specimen under mode II loading 
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adherend 

adhesive 

adherend 


Figure 5 5 Double lap joint under tensile loading 



Chapter 6 

General Conclusions 


In this Chapter, general conclusions, from linear elastic, elasto-plastic, elasto- 
visco-plastic analyses and experimental study, are listed A companson of 
computation time for different analyses are made Limitations of the present 
work and the future scope of this work are also discussed 


6.1 Results and discussion 

6.1.1 Comparison of results 

Execution time taken for linear elastic, elasto-plastic and elasto-visco-plastic 
analyses to open three pairs of nodes m the computer CONVEX C 220, are 
compared m Thble 6 1 The starting load level, material properties, physical 
dimensions, mesh refinement and crack location are maintained identical for all 
the three analyses Adhesive and adherends are considered to be Epon VIII 
and CFRP laminate, respectively The non-dimensional ratios g (= Gc/Gc), 
1 (« L(,/Lt)a.nd t (= f,/fr)are assumed to be 04, 1/2 and 1/5, respectively 
The convergence tolerance value, fracture toughness tolerance value, fracture 
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toughness, fluidity parameter and the flow rate function employed are 1 percent, 
1 percent, 0 15 kJra~^, 4 495 lO'VSec and Equation (6 2), respectively From the 
comparison of the results it is observed that the linear elastic analysis predicts a 
lower failure load as compared to both elasto-plastic and the elasto-visco-plastic 
analyses This might be due to the peak stresses calculated from elasto-plastic 
and elasto-visco-plastic analyses are reduced 




( 


ctq ) 


N 


-1 


E. 


(6 1 ) 
(6 2 ) 


where, JV is 1 426, <tq is the yield stress and the flow rule F is expressed 
as 


^ + 0r2 - + 2>rl^ 


(6 3) 


On the basis of elastic, elasto-plastic and elasto-visco-plastic analyses carried 
out, the following general conclusions are made 


• The finite element approach with paired nodes modelling the crack provides 
a simple computational tool to predict joint failure However, fine grids 
near the crack-tip are necessary for accurate estimation of strain energy 
release rates 

• The elastic analysis predicts a lower failure load as compared to the 
elasto-plastic and elasto-visco-plastic analyses Therefore, for a realistic 

analysis elasto-plastic and elasto-visco-plastic behaviour arc necessary to 
be included 
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• The general computer code developed, can model fairly accurately the 
elastic, clasto-plastic and clasto-visco-plastic behaviour of adhesively bonded 
joint 

• Stress distributions, along the overlap length of the single lap joint, are 
predicted by the present clastic and elasto-visco-plastic analysis and the 
results agree with those reported by Hiregoudar (1993) 

• Stress distributions, along the overlap length and across the adhesive 
thickness of thick adherend shear test (TAST) specimen (stepped lap joint) 
arc predicted by the clastic and elasto-visco-plastic analysis and the results 
agree with the those reported by Su and Mackie (1993) 

• In case of single edge notch (SEN) specimen, with 2688 elements and 
proper mesh refinement, error estimated in the value of fracture parameter 
is less than 1 6 percent (as compared with close form solution) 

• Strain energy release rate predicted by finite element analysis with fracture 
mechanics approach agreed well with the expenmental result reported by 
Fcrnlund and Spelt (1991c) for cracked lap shear (CLS) specimen The 
difference is less than 2 7 percent 

• Growth of plastic zone is influenced by the displacement boundary conditions 
and the size of the crack 

• Cracks occurring at the edges are more critical than those at the central 
part of the overlap region Further, m symmetne boundary conditions the 
diagonally opposite cracks are equally prone to failure 

. In case of composite laminates, introduction of 0“ ply m the immediate 
neighbourhood of adhesive results m a stronger joint Further, a proper 
choice of layup sequence results in an efficient bonded jomt 

• The results obtained from parametric study are similar for linear elastic, 
elasto-plastic and elasto-wsco-plastic analyses 
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• Use of adhesive of high Young’s modulus, high fracture toughness and 
large overlap length results in a stronger bonded joint 

• Use of small thickness of adhesive results in stronger bonded joint Free- 
edge effects are significant for joints with very thin adhesive For the 
analysis of bonded joints having very small adhesive thickness, more refined 
mesh IS preferable to account for the free-edge effect 

• Following fracture cntenon is observed to predict bond strength better 
than other basic cnteria (as compared with the experimental results) 


Gi + Gji > Gic + Giic (6 4) 

where Gjc and Guc are cntical strain energy release rates in modes I and 
II, respectively 

• The computed bond strengths employing elastic behaviour for the adhesive 
are m good comparison with those obtained from the experiments The 
difference in fracture toughness results, might be due to non-uniform 
thickness of the composite adherends and adhesive during specimen (joint) 
preparation 


6.2 Limitations of the present work 


In the present work, the bonded joints are modelled as two-dimensional plane 
stress and/or plane strain models It is assumed that the adhesive in bulk 
adhesive sample specimen and very thin layer of adhesive in the bonded joint, 
have same material properties In the non-linear analyses only adhesive material 
IS considered to be either elasto-plastic and elasto-visco-plastic Linear elastic 
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fracture mechanics (LEFM) is employed The cracks are assumed to be quasi- 
static and through width type and the crack growth is assumed to occur along 
the interface of adhesive and adherend (adhesive type failure) 


6.3 Suggestions for future work 


Present two-dimensional finite element analysis could be extended to three- 
dimensional one Thumb nail cracks could be employed Cohesive and mixed 
(both cohesive and adhesive) type of failure could be included Similar to 
overlap length, area of adhesive in contact with the adherends could be a more 
realistic parameter to study for the prediction of strength of bonded joints 
Geometrical non-lineanty could be included More number of expenments 
could help to find fracture toughnesses m pure modes (I, II and III) and mixed 
modes Prediction of failure load under adverse environmental conditions needs 
to be looked into Incorporation of non-linear fracture mechanics might be 
rational, particularly when the adhesive exhibit non-lmear behaviour Studies 
on use of primers in bonded joints could help to analyse veiy strong adhesive 

joints 
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Table 6 1 Comparison of results obtained from linear elastic, elasto-plastic and 
elasto-visco-plastic analyses 


crack length (mm) 

elastic 

failure loads in kN/mra 
elasto-plastic elasto-visco-plastic 

0 055 

13 824 

14496 

14 400 

0 219 

11612 

11379 

11592 

0 492 

10 886 

10867 

10 896 

execution time (Seconds) in Convex C 220 


2 432 

3 574 

4513 



References 


Adams, R D (1989) Strength prediction for lap joints, especially with composite 
adherends, a review Joumol of Adhesion, 30 219-242 

Adams, R D (1992) Strength predicted in bonded joints current thoughts and 
research at the university of Bnstol, U K Joumol of Adhesion, 37 47-50 

Adams, R D and Hams, J A (1987) The influence of local geometiy on 
the strength of adhesive joints Intemotionol Joumol of Adhesion and Adhesives, 
7(2) 69-80 

Adams, R D and Mallick, V (1993) The effect of temperature on the 
strength of adhesively-bonded composite-aluminium jomts Journal of Adhesion, 
43 17-33 

Adams, R D and Wake, W C (1984) Structural odhesive joints m engineenng 
Elsevier Applied Science Publishers London and New York 

Adams, R D, Coppendale, J, Mallick, V and H, Al-Halman (1992) The 
effect of temperature on the strength of adhesive joints Intemotionol Joumol 
of Adhesion and Adhesives, 12(3) 185-190 

Aivazzadeh, S and Verchery, G (1986) Stress analysis at the interface adhesive 
joints by special finite elements Intemotionol Joumol of Adhesion and Adhesives, 
6(4) 185-188 

Aivazzadeh, S , Bichara, M , Ghazal, A and Verchery, G (1988) Special 
mixed finite elements for mterfacial stress analysis of adhesively bonded joints 
Adhesively Bonded Joints Testing, Anolysis ond Design, WS Johnson, Ed , ASTM 
STP, 981 133-144 

Allen, K W and Shanahan, M E R (1975) The creep behaviour of structural 
adhesive jomts I Journal of Adhesion, 7 161-174 

Allen, K W and Shanahan, M E R (1976) The creep behaviour of structural 
adhesive joints II Journal of Adhesion, 8 4^56 


205 



206 


REFERENCES 


Alwar, R S and Nagaraja, Y R (1976) Elastic analysis of adhesive butt joints 
Journal of Adhesion, 6 279-287 

Bassani, J L and Qu, J (1989) Finite crack on bimatenal and bicrystal 
interfaces Journal of Mechanical Physics and Solids, 37(4) 435-453 

Bigwood, D A and Crocombe, A D (1989a) Bonded joint design analysis 
Adhesion - 13 Ed K W Allen, Chapter - 11, Elsevier Applied Science Publishers 
New York 163-187 

Bigwood, D A and Crocombe, A D (1989b) Elastic analysis and engineering 
design formulae for bonded joints International Journal of Adhesion and 
Adhesives, 9(4) 229-242 

Bigwood, D A and Crocombe, A D (1990) Non-linear adhesive bonded joint 
design analyses International Journal of Adhesion and Adhesives, 10(1) 31-41 

Biot, M A (1954) Theory of stress-strain relations in anisotropic viscoelasticity 
and relaxation phenomena Journal of Applied Physics, 25(11) 1385-1391 

Bishopp, J A (1992) The chemistry and properties of a new generation 
of toughened epoj^ matrices International Journal of Adhesion and Adhesives, 
12(3) 178-184 

Broek, D (1984) Elementaiy Engineering Fracture Mechanics Martinus Nijhoff 
Publishers Boston 

Brust, F W, Nakagami, M and Gilles, P (1990) Companson of elastic-plastic 
fracture mechanics techniques Fracture Mechanics Tenty-first Symposium, J P 
Gudas, J^ Joyce and E M Hackett, Eds ASTM STP, 1074 448-469 

Carlsson, L A and Pipes, R B (1987) Experimental characterisation of advanced 
composite materials Prentice-Hall Inc 

Carlsson, L A , Gillespie, J W (Jr) and Pipes, R B (1986) On the design and 
analysis of end notched specimen (ENF) specimen for mode II testing Journal 
of Composite Materials, 20 594-604 

Chai, H (1986) Bond thickness effect m adhesive joints and its significance 
for mode I interlaminar fracture of composites Composite Materials Testing 
and Design (Seventh Conference), J M Whitney, Ed , ASTM STP, 893 209-231 

Chai, H (1992) Experimental evaluation of mixed-mode fracture in adhesive 
bonds Experimental Mechanics, 32(2) 296-303 



REFERENCES 


207 


Chai, H (1993) Observation of deformation and damage at the tip of cracks 
in adhesive bonds loaded in shear and assessment of a criterion for fracture 
Intemotional Journal oj Fracture, 60(4) 311-326 

Chalkley, P D and Chiu, W K (1993) An improved method for testing the 
shear stress/strain behaviour of adhesives Intemotional Journal of Adhesion and 
Adhesives, 13(4) 237-242 

Chen, D and Cheng, S (1983) An analysis of adhesive-bonded single-lap 
joints Trons ASME Journal of Applied Mechanics, 50 109-115 

Chen, D and Cheng, S (1990) Stress distnbution in plane scarf and butt joints 
Trans ASME Joumol of Applied Mechanics, 57 78-83 

Chen, D and Cheng, S (1992) Torsional stress in tubular lap joints Intemotinol 
Journal Solids Structures, 29(7) 845-853 

Cheng, S , Chen, D and Shi, Y (1991) Analysis of adhesive-bonded jomts 
with nonidentical adherends Joumol of Engineering Mechonics, 117(3) 605-623 

Cheriy, B W and Ye, Y Q (1992) The behaviour of high-temperature 
anaerobic adhesives International Joumol of Adhesion and Adhesives, 12(3) 206- 
210 

Choi, H J and Thangjitham, S (1993) Thermally-induced interlaminar crack- 
tip singularities in laminated anisotropic composites Intemationol Journal of 
Fracture, 60(4) 349-361 

Conley, K M , Gu, W, Ritter, J E and Lardner, T J (1992) Observations 
on finger-like crack growths at a urethane acrylate/glass interface Journal of 
Adhesion, 39 173-184 

Conley, K M , Ritter, J E and Lardner, T J (1992) Subcritical crack growth 
along epoxy/glass interfaces Joumol of motenal Research, 7(9) 2621-2629 

Crocombe, A D and Bigwood, D A (1992) Development of full elasto-plastic 
adhesive joint design analysis Journal of Strom Analysis, 27(4) 211-218 

Crocombe, A D and Evan, I E J (1988) The interaction of adhesive jomt 
strength of adherend cladding Journal of Adhesion, 26 199-213 

Crocombe, A D , Bigwood, D A and Richardson, G (1990) Analysing 
structural adhesive joints for failure Intemotionol Journal of Adhesion ond 
Adhesives, 10(3) 167-178 



208 


REFERENCES 


Crompton, J S and Clark, J D (1989) Finite element modelling applied to 
crack propagation studies in bonded components Adhesion - 13 Ed K W Allen, 
Chapter-12, Elsevier Applied Science Publishers New York 188-200 

Czamocki, P and Piekarski, K (1986) Non-linear numerical stress analysis of 
a symmetnc adhesive bonded lap joint International Journal of Adhesion and 
Adhesives, 6(3) 157-160 

Davies, P, Kausch, H H , Williams, J G and Kinloch, A J (1992) Round- 
robin interlaminar fracture testing of carbon-fibre-reinforced epoxy and PEEK 
composites Composite Structures, 43(2) 129-136 

Davies, R and Khalil, A A (1990) Design and analysis of bonded double 
containment corner joints Intemationol Journal of Adhesion and Adhesives, 
1(1) 25-30 


Devitt, D F, Schapery, R A and Bradley, W L (1980) A method 
for determining the mode I delamination fracture toughness of elastic and 
viscoelastic composite matenals Journal of Composite Materials, 14 270-285 

Dom, L and Weiping, L (1993) The stress state and failure properties 
of adhesive-bonded plastic/metal joints International Journal of Adhesion and 
Adhesives, 13(1) 21-31 

Dung, N L (1992) Three-dimensional void growth in plastic materials 
Mechanics Research Communications, 19(3) 227-235 

Edde, F and Verreman, Y (1992) On the fracture parameters in a clamped 
cracked lap shear adhesive joint International Journal of Adhesion and Adhesives, 


IMund, U and Warbnng, A (1990) Analysis of elastic and elastic-plastic 
adhesive joints using a mathematical programming approach Computer Methods 
in Applied Mechanics and Engineenng, 78 19-47 


Felippa, C A (1975) Solution of linear equations with sltyline stored symmetnc 
matrix Computers and Structers, 5(1) 13-28 


Fernlund, G and Spelt, J K (1991a) Analytical method for calculating adhesive 
joint tracture parameters Engineering Fracture Mechanics, 40(1) 119-132 


Femlund, G and Spelt, J K (1991b) 
adhesive joints part 1, Analytical method 
Adhesives, 11(4) 213-220 


Failure load prediction of structural 
International Journal of Adhesion ond 



REFERENCES 


209 


Fernlund, G and Spelt, J K (1991c) Failure load prediction of structural 
adhesive joints part 2 experimental study International Journal of Adhesion 
and Adhesives, 11(4) 221-227 

Fmdlater, D (1987) A preliminaiy design approach for adhesively bonded 
joints Intemotionol Journal of Adhesion ond Adhesives, 7(3) 129-134 

Fujii, T (1993), Dynamic response of sandwich beams with an adhesive damping 
layer (general Maxwell method for a viscoelastic adhesive layer) International 
Journal of Adhesion and Adhesives, 13(3) 201-209 

Gilchrist, M D and Smith, R A (1993) Fatigue growth of cohesive defects 
in T-peel joints Journal of Adhesion, 42 179-190 

Goland, M and Reissner, E (1944) The stresses m cemented joints Trans 
ASME Journal of Applied Mechanics, 11(1) A17-A27 

Groth, H L (1986) Calculation of stresses in bonded joints using the 
substructuring technique International Journal of Adhesion and Adhesives, 
6(1) 31-35 

Groth, H L (1988a) Prediction of failure loads of adhesive joints using the 
singular intensity factor Fracture Mechanics Eighteenth Symposium, D T Read 
and R P Reed, Eds , ASTM STP, 945 278-284 

Groth, H L (1988b) A test specimen with constant stress intensity factor 
Adhesively Bonded Joints Testing, Analysis and Design, W S Johnson, Ed , ASTM 
STP, 981 98-104 

Groth, H L (1990) Viscoelastic and viscoplastic stress analysis of adhesive 
joints International Journal of Adhesion ond Adhesives, 10(3) 207-213 

Groth, H L and Brottare, I (1988) Evaluation of stress intensity factors in 
elastic-plastic matenals Journal of Testing and Evaluation, 16(3) 291-297 

Groth, H L and Brottare, I (1989) Apparent stiffness of a butt joints with 
a thick adhesive layer and elastic-plastic adhesives Journal of Testmg and 
Evaluation, 17(2) 131-134 

Groth, H L and Nordlund, P (1991) Shape optimization of bonded joints 
International Journal of Adhesion and Adhesives, 11(4) 204-211 

Habib, M, Aivazzadeh, S and Verchery, G (1988) Edge effect analysis 
in adhesively bonded composite tubes using axisymmetnc interface elemen^ 
Proceedings International Conference Southampton, on Computer Aided Desi 
Composite Matenal Technology, Eds C A Brebbia, W P deWilde ond W R 
415-428 



210 


REFERENCES 


Hanneman, S E and Kinra, V K (1992a) A new technique for ultrasonic 
nondestructive evaluation of adhesive joints Part I theory Experimental 
Mechanics, 32(2) 323—331 

Hanneman, S E and Kinra, V K (1992b) A new technique for ultrasonic 
nondestructive evaluation of adhesive joints Part II expenment Experimental 
Mechanics, 32(2) 332-339 

Harj-Rabab, Y , Aivazzadah, S and Verchery, G (1990) Edge effect analysis 
and influence of defects in adhesively bonded composite joints using interface 
finite elements Development and Design with Advanced Materials, Eds G C 
Sih, S V Hoa and J J Pindera, Elsevier Science Publishers B V, The Netherlands 
129-136 

Harris, J A and Adams, R D (1985) An assessment of the impact performance 
of bonded joints for use m high energy absorbing structures Proceedings of the 
Institution of Mechanical Engineers, 199(C2) 121-131 

Hart-Smith, L J (1973) Adhesive-bonded double-lap joints NASA Technical 
Report, CR 112235 

Hart-Smith, L J (1981) Stress analysis a continuum mechanics approach 
Developments in Adhesives-2, Ed A J Kmloch, Applied Science Publishers London 

Hashim, S A , Cowling, M J and Winkle, I E (1990) Design and assessment 
methodologies for adhesively bonded structural connections International 
Journal of Adhesion and Adhesives, 10(3) 139-145 

Hatton, T (1991) A stress-singulanty-parameter approach for evaluating the 
adhesive strength of single-lap joints JSME International Journal Senes I, 
34(3) 326-331 

Hatton, T, Sakata, S and Murakami, G (1989) A stress singulanty parameter 
approach for evaluating the mterfacial reliability of plastic encapsulated LSI 
devices Journal of Electronic Packaging, 111 243-248 

Hennksen, M (1984) Nonlinear viscoelastic stress analysis — a finite element 
approach Computers and Structers, 18(1) 133-139 

Hiregoudar, S (1993) Viscoplastic and geometric nonlinear finite element 
analysis of adhesively bonded lap joint M Sc Thesis Report, Dept of Civil 
Engineenng, Indian Institute of Science, Bangalore 

Ikegami, K, Takeshita, T, Mastruo, K and Sugibayashi, T (1990) Strength 
of adhesively bonded scarf joints between glass-reinforced plastics and metals 
IntemaUonal Journal of Adhesion and Adhesives, 10(3) 199-206 



REFERENCES 


211 


Imanaka, M , Kishimoto, W, Okita, K , Nakayama, H and Nagai, H (1989) 
Fatigue life estimation of adhesive bonded shaft joints Intemotionol Journal of 
Fracture, 41 223-234 

Ishikawa, H , Yuuki, R , Chung, N Y and Nakano, S (1991) Mixed mode 
fracture criteria on adhesive joints Proceedings of Intemotional Symposium on 
Mixed-Mode Fracture and Fatigue at Vienna AUSTRIA, during July 15-19, 1991 

Janardhana, M N, Brown, K C, Jones, R and Paul, J (1986) Effect 
of debonding adhesively bonded composite to metal joints in compression 
Composite Structures, 5(1) 1-14 

Jeandrau, J P (1986) Intrinsic mechanical characterization of structural 
adhesives International Journal of Adhesion and Adhesives, 6(4) 229-231 

Jeandrau, J P (1991) Analysis and design data for adhesively bonded joints 
International Journal of Adhesion and Adhesives, 11(2) 71-79 

Jiao, D and Rose, J L (1991) An ultrasonic interface layer model for bond 
evaluation Journal of Adhesion Science and Technology, 5(8) 631-646 

Johnson, W S (1987) Stress analysis of the cracked-lap-shear specimen An 
ASTM round-robin Journal of Testing ond Evaluation, 15(6) 303-324 

Jones, R M (1975) Mechanics of Composite Motenols Scnpta Book Comapany 
Washington D C 

Kairouz, k C and Matthews, F L (1990) Mechamsms of failure in bonded 
CFRP single lap joints with different stacking sequences Proceeding? of the 
Institution of Mechanical Engineers, FRC’ (1990), Fourth International Conference, 
C400/014 47-59 

Kanchi, M B , Zienkiewicz, O C and Owen, D R J (1978) The visco-plastic 
approach to problems of plasticity and creep mvolvmg geometncal nonhnear 
effects International Journal for Numerical Methods in Engineering, 12 169-181 

Keary, P E, Ilcewicz, L B, Shaar, C and Tlostle, J (1985) Mode I 
interlaminar fracture toughness of composites usmg slender double cantilever 
beam specimens Joumol of Composite Materials, 19 154—176 

Khalil, A A and Bayoumi, M R (1991) Effect of loading rate on fracture 
toughness of bonded joints International Joumol of Adhesion and Adhesives, 
11(1) 25-29 

Kim, K S, Kim, W T, Lee, D G and Jun, E J (1992) Opbm^ tubular 
adhesive-bonded lap joint of the carbon fiber epoxy composite shaft Composite 
Structures, 21(3) 163-176 



212 


REFERENCES 


Krenk, S (1992) Energy release rate of symmetric adhesive joints Engineering 
Fracture MecHomcs, 43(4) 549-559 

Lee, A W and Selby, J S (1992) Adhesive selection by means of the pal 
program Computerisation and Networking of Materials Databases Third Volume, 
T I Barry and K W Reynard, Eds , ASTM STP, 1140 

Lee, R J, Davidson, R and McCarthy, J C (1989) Composite to metal 
joining for transport application Adhesion - 13 Ed K W Allen, Chapter-10, 
Elsevier Applied Science Publishers New York 142-162 

Liechti, K M and Chai, Y S (1991) Biaxial loading experiments for determining 
mterfacial fracture toughness Trons ASME Journal of Applied Mechanics, 58 680- 
687 

Long, R S (1991) Static strength of adhesively bonded ARALL-1 joints 
Journal of Composite Materials, 25 391-415 

Mackie, R I and Vardy, A E (1990) Applying the coin-tap test to adhesives 
in civil engineering a numerical study International Journal of Adhesion and 
Adhesives, 10(3) 215-220 

Mall, S and Johnson, W S (1986) Characterization of mode I and mixed— mode 
failure of adhesive bonds between composite adherends Composite Materials 
Testing and Design (Seventh Conference), J M Whitney, Ed , ASTM STP, 893 322- 
334 

Mall, S and Ramamurthy, G (1989) Effect of bond thickness on fracture and 
fatigue strength of adhesively bonded composite joints International Journal of 
Adhesion and Adhesives, 9(1) 33-37 

Marcolefas, S , Kostopoulos, V and paipetis, S A (1991) Non— linear analysis 
of a metal— to-coraposite scarf joint Intemationol Journal of Mechanical Science, 
33(12) 961-973 

Mecklenburg, M F, Arab, C O , McNamara, D , Hand, H and Joyce, J A 
(1990) Adhesive fracture testmg Fracture Mechanics Tenty-first Symposium, 

J P Gudas, J A Joyce ond E M Hackett, Eds ASTM STP, 1074 307-321 

Miskioglu, I, Villmann, CR, Pawloski, J S and Panseau, D (1991) A 
photoelastic and FEM analysis of mterfacial crack propagation Experimental 
Mechanics, 31 135-139 

Munz, D and Yang, Y Y (1992) Stress singularities at the interface in bonded 
dissimilar materials under mechanical and thermal loading Trans ASME Journal 
of Applied Mechanics, 59 857-861 



REFERENCES 


213 


Nagaraja, Y N and Alwar, R S (1980) Viscoelastic analysis of an adhesive- 
bonded plane lap joint Computers and Structers, 11 621-627 

Needleman, A (1990) Analysis of interfacial failure Applied Mechanics Review, 
43(5) S274-S275 

Ojalvo, I U and Eidinoff, H L (1978) Bond thickness eifects upon stresses 
in single-lap adhesive joints AIAA Journal, 16(3) 204-211 

Owen, D R J and Hinton, E (1980) Finite elements in plasticity Pineridge 
Press, Swansea U K 

Parker, B M (1990) The strength of bonded carbon fibre composite joints 
exposed to high humidity International Journal of Adhesion ond Adhesives, 
10(3) 187-191 

Parker, B M (1993) Environmental durabihty of aluminium joints with different 
pretreatments International Journal of Adhesion and Adhesives, 13(1) 47-51 

Pilarski, A and Rose, J L (1988) A transverse-wave ultrasonic oblique- 
mcidence technique for interfacial weakness detection in adhesive bonds Joumol 
of Applied Physics, 63(2) 300-307 

Pmdera, J T and Wang, G (1990) Cnticality of local three dimensional stresses 
on a sy mm etric adhesive bonded joint Development and Desigfi with Advonced 
Materials, Eds G C Sih, S V Hoo ond J J Pindero, Elsevier Science Publishers 
B V 119-128 

Qu, J and Bassani, J L (1989) Crack on bimatenal and biciystal interfaces 
Journal of Mechonical Physics arid Solids, 37(4) 417-433 

Qu, J and Bassani, J L (1993) Interface fracture mechanics for anisotropic 
matenals TVans ASME Joumol of Applied Mechonics, 60 422-431 

Qu, J and Li, Q (1991) Interfacial dislocation and its apphcations to interface 
cracks in anisotropic bimatenals Journal of Elastunty, 26 169-195 

Ramamurthy, T S and Rao, A K (1984) Analysis of bonded joints 
with arbitraiy adherend shapes and adhesive non-hneanties Joumol of the 
Aeronouticol Society of Indio, 36(1) 29-36 

Rao M D and Crocker, M J (1990) Analytical and experimental study of the 
vibration of bonded beams with a lap joint Trans ASME Journal of Vibrations 
and Acoustics, 112 444-451 



214 


REFERENCES 


Ratwani, M M and Kan, H (1986) Analysis of delammation propagation in 
composite layer of metal-to-composite stepped-Iap joints Composite Materials 
Testing and Design (Seventh Conference), J M Whitney, Ed , ASTM STP, 893 308- 
321 

Ratwani, M M , Kan, H P and Liu, D D (1982) Time-dependent adhesive 
behavior effects in a stepped lap joint AIAA Journal, 20(5) 734-736 

Reddy, J N and Roy, S (1988) Non-lmear analysis of adhesively bonded 
joints International Journal of Non-Linear Mechanics, 23(2) 97-112 

Richardson, G , Crocombe, A D and Smith, P A (1993) A comparison of two- 
and three-dimensional finite element analyses of adhesive joints International 
Journal of Adhesion and Adhesives, 13(3) 193-200 

Ritter, J E and Conley, K (1992) Moisture-assisted crack propagation 
at polymer/glass interfaces Intemotional Journal of Adhesion and Adhesives, 
12(4) 245-250 

Rossettos, J N and Zang, E (1993) On the peak shear stress in adhesive joint 
with voids Trons ASME Journal of Applied Mechanics, 60 559-560 

Roy, S and Reddy, J N (1988a) Fimte-element models of viscoelasticity 
and diffusion in adhesively bonded joints Intemationol Journal for Numerical 
Methods in Engineenng, 26 2531-2546 

Roy, S and Reddy, J N (1988b) A finite element analysis of adhesively bonded 
composite joints with moisture diffusion and delayed failure Computers and 
Structers, 29(6) 1011-1031 

Rybicki, E F and Kannmen, M F (1977) A finite element calculation of 
stress intensity factors a modified crack closure integral Engineenng Fracture 
Mechanics, 9 931-938 

Sawa, T and Kobayashi, T (1988) The strength of joints combining an adhesive 
with a bolt Journal of Adhesion, 25 269-280 

Sawa, T, Ishikawa, H and Muto, K (1992) The strength of joints combining 
adhesives with bolts JSME International Journal Senes I, 35(1) 38-44 

Sawa, T, Temma, K and Ishikawa, H (1989) Three-dimensional stress analysis 
of adhesive butt joints of solid cylinders subjected to external tensile loads 
Joumol of Adhesum, 31 33-43 

Sawa, T, Temma, K and Tfeunoda, Y (1989) Axisymmetric stress analysis of 
adhesive butt joints of dissimilar solid cylinders subjected to external tensile 
loads Intematumol Journal of Adhesion ond Adhesives, 9^) 161-169 



REFERENCES 


215 


Schapery, R A (1984) Correspondence principles and a generalized J integral 
for large deformation and fracture analysis of viscoelastic media Intemationol 
Journal of Fracture, 25 195-223 

Sen, J K and Jones, R M (1980a) Stresses in double-lap joints bonded with 
a viscoelastic adhesive part I theory and experimental corroboration AIAA 
Journal, 18(10) 1237-1244 

Sen, J K and Jones, R M (1980b) Stresses in double-lap joints bonded 
with a viscoelastic adhesive part II parametric study and joint design AIAA 
Journal, 18(11) 1376-1382 

Sethuraman, R and Maiti, S K (1989) Finite element analysis of doubly 
bonded crack-stiffened panels under mode I or mode II loading Engineering 
Fracture Mechanics, 34(2) 465-475 

Su, N and Mackie, R I (1993) IWo-dimensional creep analysis of structural 
adhesive joints International Journal of Adhesion and Adhesives, 13(1) 33-40 

Su, N, Mackie, R I and Harvey, W J (1992) The effects of aging and 
environment on the fatigue life of adhesive joints Intemationol Journal of 
Adhesion and Adhesives, 12(2) 85-93 

Suo, Z (1990) Failure of bnttle adhesive joints Applied Mechanics Review, 
43(5) S276-S279 

Temma, K, Sawa, T and Iwata, A (1990) Two-dimensional stress analysis 
of adhesive butt joints subjected to cleavage loads International Joumol of 
Adhesion and Adhesives, 10(4) 285-293 

Temma, K , Sawa, T and Tsunoda, Y (1990) Three-dimensional stress analysis 
of adhesive butt joints with disbonded areas and spew fillets Intemotional 
Joumol of Adhesion and Adhesives, 10(4) 294-300 

Temma, K, Sawa, T, Uchida, H and Nakano, Y (1991) A two-dimensional 
stress analysis of butt adhesive joints having a circular hole defect in the adhesive 
subjected to external bending moments Joumol of Adhesion, 33 133-147 

Thomsen, O T (1992) Elasto-static and elasto-plastic stress analysis of adhesive 
bonded tubular lap joints Composite Structures, 21(4) 249-259 

Thomson, R B and Thomson, D O (1991) Past expenences in the development 
of tests for adhesive bond strength Joumol of Adhesion Science ond Technology, 
5(8) 683-599 



216 


REFERENCES 


Tiwan, A , Hienneke, I , Edmund, G and Duke, J C (1991) Acousto- 
ultrasonic(AU) technique for assuring adhesive bond quality Journal of Adhesion, 
34 1-15 

Volkerson, O (1938) Rivet strength distnbution m tensile stressed rivet joints 
with butt strap cross-section Luftohforschung, 15 41-47 

Wang, A S D and Crossman, F W (1982) Fracture mechanics of sub-Iaminate 
cracks Technicol Report F 49620-79-C-0206, AFOSR, Drexel Univ , Philadelphia 

Wassell, G C, Clark, J D, Crompton, J S and Dickson, R F (1991) 
Fatigue within adhesive bonds International Journal of Adhesion and Adhesives, 
11(2) 117-120 

Whitcomb, J D and Dattaguru, B (1984) User's manual for GAMNAS- 
Geometrical and Material Nonlinear Analysis of Structures NASA Technical 
Memorandum, 85734 

Whitcomb, J D and Woo, K (1993) Analysis of debond growth in tubular joints 
subjected to tension of flexural loads Computers ond Structers, 46(2) 323-329 

Wilkms, D J , Eisenmann, J R , Camin, R A , Margolis, W S and Benson, R A 
(1982) Characterizing delamination growth in graphite epoxy Characterizing 
Delamination Growth m Graphite Epoxy, K L Reifsnider, Ed , ASTM STP, 775 168- 
183 

Wingfield, J R J (1993) Treatment of composite surfaces for adhesive bonding 
International Journal of Adhesion and Adhesives, 13(3) 151-156 

Yadagin, S , Reddy, C P and Reddy, T S (1987) Viscoelastic analysis of 
adhesively bonded joints Computers and Structers, 27(4) 445-454 

Yamada, Y , Yoshimura, N and Sakuri, T (1968) Plastic stress— strain matrix 
and Its application for the solution of elastic-plastic problems by the finite 
element method Intemationol Joumol of Science, 10 343-354 

Zenkert, D (1991) Strength of sandwich beams with interface debondings 
Composite Structures, 17(4) 331-350 

Zienkiewicz, O C and Cormeau, I C (1974) Visco-plasticity-plasticily and 
creep m elastic solids - A unified numerical solution approach International 
Joumol for Numerical Methods in Engineering, 8 821-845 

^enkiewicz, O C and Taylor, R L (1989) The finite element method, volume 1 
McGraw-Hiil Book Co New York, fourth edition 



REFERENCES 


217 


Zicnkicwicz, O C and Taylor, R L (1991) The jinite element method, volume 2 
McGraw-Hill Book Co New York, fourth edition 

Zicnkicwicz, O C, Owen, D R J and Cormeau, I C (1974) Analysis of 
viscoplastic effects in pressure vessels by the finite element method Nuclear 
Engineering and Design, 28 278-288 

Zienkiewicz, O C , Valliappan, S and King, I S (1969) Elasto-plastic solutions 
of engineering problems ’initial stress’, finite element approach International 
Joumol for Numerical Methods in Engineenng, 1 75-100 



